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Abstract 



This review article is devoted to analyze the main properties characterizing the cosmo- 
logical singularity associated to the homogeneous and inhomogeneous Mixmaster model. 
After the introduction of the main tools required to treat the cosmological issue, we re- 
view in details the main results got along the last forty years on the Mixmaster topic. We 
firstly assess the classical picture of the homogeneous chaotic cosmologies and, after a 
presentation of the canonical method for the quantization, we develop the quantum Mix- 
master behavior. Finally we extend both the classical and quantum features to the fully 
inhomogeneous case. Our survey analyzes the fundamental framework of the Mixmaster 
picture and completes it by accounting for recent and peculiar outstanding results. 
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1 Preface 



The formulation by Belinski, Khalatnikov and Lifshitz (BKL) at the end of the Sixties 
about the general character of the cosmological singularity represents one of the most 
important contributions of the Landau school to the development of Theoretical Cosmol- 
ogy after the work by Lifshitz on the isotropic Universe stability The relevance of the 
BKL work relies on two main points: on one hand, this study provided a piece wise an- 
alytical solution of the Einstein equations, improving and implementing the topological 
issues of the Hawking-Penrose theorems; on the other hand, the dynamical features of 
the corresponding generic cosmological solution exhibit a chaotic profile. 
Indeed the idea that the behavior of the actual Universe is some how related to an inho- 
mogeneous chaotic cosmology appears very far from what physical insight could suggest. 
In fact, the Standard Cosmological Model is based on the highly symmetric Friedman- 
Robertson-Walker (FRW) model, and observations agree with these assumptions, from 
the nucleosynthesis of light elements up to the age when structures formation takes place 
in the non-linear regime. 

Furthermore a reliable (but to some extent, model-dependent) correspondence between 
theory and data exists even for an inflationary scenario and its implications. Several ques- 
tions about the interpretation of surveys from which inferring the present status of the 
matter distribution and dynamics of the Universe remain open. However, in the light of 
an inflationary scenario, a puzzling plan appears unless the observation of early cosmo- 
logical gravitational waves become viable, since such dynamics cancels out much of the 
pre-existent information thus realizing local homogeneity and isotropy. 
The validity of the BKL regime must be settled down just in this pre-inflationary Universe 
evolution, although unaccessible to present observations. These considerations partly 
determined the persistent idea that this field of investigation could have a prevalent aca- 
demical nature, according to the spirit of the original derivation of the BKL oscillatory 
regime. 

In this paper, we review in details and under an updated point of view the main fea- 
tures of the Mixmaster model. Our aim is to provide a precise picture of this research line 
and how it developed over the last four decades in the framework of General Relativity 
and canonical quantization of the gravitational field. 

We start with a pedagogical formulation of the fundamental tools required to under- 
stand the treated topics and then we pursue a consistent description which, passing 
through the classical achievements, touches the review of very recent understandings 
in this field. The main aim of this work is to attract interest, especially from researchers 
working on fundamental physics, to the fascinating features of the classical and quantum 
inhomogeneous Mixmaster model. Our goal would be to convince the reader that this 
dynamical regime is physically well motivated and is not only mathematical cosmology 
effectively concerning the very early stages of the Universe evolution, which matches the 
Planckian era with the inflationary behavior. 

The main reason underlying the relevance of the Mixmaster dynamics can be identified 
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as follows: the pre-inflationary evolution fixes the initial conditions for all the subsequent 
dynamical stages as they come out from the quantum era of the Universe, and hence the 
oscillatory regime concerns the mechanism of transition to a classical cosmology provid- 
ing information on the origin of the actual notion of space-time. 

Moreover, we emphasize three points strongly supporting that the very early Universe 
dynamics was described by more general paradigms than the isotropic model, indeed 
concerning the generic cosmological solution. 

- The FRW model is backward unstable with respect to tensor perturbations, which 
increase as the inverse of the scale factor. Recent studies have shown how such 
instability holds for scalar perturbations too, as far as the presence of bulk viscosity 
is taken into account. 

- The inflationary scenario offers an efficient isotropization mechanism, able to rec- 
oncile the primordial inhomogeneous Mixmaster with the local high isotropy of the 
sky sphere at the recombination age. 

- In the Planckian era, the quantum fluctuations can be correlated at most on the 
causal scale, thus we should regard global symmetries as approximated toy models. 
In this respect, the assumption that the Universe was born in the homogeneous 
and isotropic configuration does not appear well grounded. On a more realistic 
point of view, the quantum dynamics has to be described in absence of any special 
symmetry, i.e. by the inhomogeneous Mixmaster Universe. 

This article is organized as follows: 
in Section 1 we review the fundamental tools for geometrodynamics, with particular rele- 
vance to the Hamiltonian formulation of General Relativity and to Singularity Theorems. 
In Section 2 we widely discuss the classifications of the homogeneous spaces and in par- 
ticular the dynamics of the types I, II and VII of the Bianchi classification. 
In Section 3 the classical features of the Mixmaster model are analyzed in details, both in 
the field equations formalism and in the Hamiltonian one. Particular attention is paid to 
the chaotic properties and to the cosmological implications of the homogeneous dynam- 
ics, also pointing out the effects on chaos of matter fields and of the number of dimen- 
sions. 

In Section 4 we face the quantization of the Mixmaster model in different framework: af- 
ter reviewing the standard Wheeler-DeWitt approach, we derive the full energy spectrum 
in the Dirac quantization scheme, and finally we discuss the more recent Loop Quantum 
Gravity approach and the Generalized Uncertainty Principle. 

In Section 5 we generalize the previous results to the inhomogeneous case, both on the 
classical and on the quantum level. We focus the attention also on several important top- 
ics of the Generic Cosmological Solution of the Einstein equations and its identification 
with the Mixmaster model. 



10 



2 Fundamental Tools 



In this Section we will analyze and discuss some of the fundamental features of the clas- 
sical theory of gravity, i.e. the Einstein theory of General Relativity (GR). In particular, 
after presenting the Einstein field equations, the way how the macroscopic matters fields 
are treated in such theory and a discussion about the tetradic formalism, we will analyze 
in some details the Hamiltonian formulation of the dynamics, which will be fundamental 
in the following. Then we will discuss the synchronous reference and finally review the 
fundamental singularity theorems although without entering in the rigorous proofs. 

Let us fix our convention, following Landau and Lifschitz B19|) : the space-time indices 
are given by the Latin letters in the middle of the alphabet, i.e. i,j,k... and run from 
to 3 (the Lorentzian indices are labeled as a, b, c). The spatial indices are the Greek 
letters a,/3,7... and run as 1,2,3. We adopt the signature (+, — , — , — ), unless differently 
specified. 



2.1 Einstein Equations 

The main issue of the Einstein theory of gravity is the dynamical character of the space- 
time metric, described within a fully covariant scheme. Assigned a four-dimensional 
manifold M, endowed with space-time coordinates x l and a metric tensor gij(x l ), its line 
element reads as 

ds 2 = gijfatdxL (2.1.1) 

This quantity fixes the Lorentzian notion of distances. The motion of a free test particle 
on M corresponds to the geodesic equation 

^ + I>/V = 0, (2.1.2) 

where u 1 = dx'/ds is its four-velocity, defined as the vector tangent to the curve, and 
= g m Tji m are the Christoffel symbols given by 

1 

r,' h„ = r, j, n = - (Bjgim + digj m ~ dmglj) ■ (2.1.3) 

The geodesic character of a curve requires to deal with a parallelly transported tangent 
vector u l . However, for a Riemannian manifold this curve extremizes the distance func- 
tional, i.e. it is provided by the variational principle 

5 [ ds = 5 [ dsjgijuhii = 0. (2.1.4) 
Jm Jm v 

If a test particle has zero rest mass, its motion is given by ds = and therefore an affine 
parameter must be introduced to describe the corresponding trajectory. 
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The equivalence principle is here recognized as the possibility to have vanishing Christof- 
fel symbols at a given point of M (or along a whole geodesic curve). The space-time 
curvature is ensured by a non-vanishing Riemann tensor 

R 1 ju = - 9/Ih + r)"r^ /c - Tf k T l mP (2.1.5) 

which has the physical meaning of tidal forces acting on two free-falling observers and 
their effect is expressed by the geodesic deviation equation 

u'V/^VfcS'') = R' ji m v)u l s m , (2.1.6) 

where s' is the connecting vector between two nearby geodesies. 
The Riemann tensor obeys the algebraic cyclic sum 

Rijki + Rnjk + Rikij = o (2.1.7) 

and the first order equations 

V w %, + VlRijmk + VfeRy/m = 0. (2.1.8) 

The constraint in Eq. (|2.1.8[) is called the Bianchi identity being identically satisfied by 
the metric tensor, meanwhile it is an equation for the Riemann tensor where covariant 
derivatives are expressed in terms of the metric. Contracting the Bianchi identity with 
g lk gfl, we get the equation 

VjG{ = 0, (2.1.9) 



where the Einstein tensor reads as 



dj^Rij-^Rgij, (2.1.10) 



in terms of the Ricci tensor Rjj = R l m and of the scalar of curvature R = g'^Rij- 
In order to get the Einstein equations in the presence of matter fields described by a 
Lagrangian density A m , we must fix a proper action for the gravitational field, i.e. for 
the metric tensor of the manifold M. A gravity-matter action, satisfying the fundamen- 
tal requirements of a covariant geometrical physical theory of the space-time, takes the 
Einstein-Hilbert matter form 



S = S £ _ H + S m = -— / d"x^{R-2KA m ) , (2.1.11) 
2k Jm 

where g is the determinant of the metric tensor gu and k is the Einstein constant. The 
variation of action (|2.1.11[) with respect to g' ] can be expressed in terms of the vector 5WJ 

as ma 

SW' = g' k 6T' lk -g l 'ST k lk , (2.1.12) 

for which the relation 

^'SRij = Vj5W (2.1.13) 
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holds (for detailed calculation see 1 5051 12421 13191 1349 \ ). providing the field equations in 
presence of matter 

Gij = KTij. (2.1.14) 
Here T» denotes the energy-momentum tensor of the matter field and reads as 

As a consequence of the Bianchi identity we find the conservation law V,T- = 0, which 
describes the motion of matter and arises from the Einstein equations. The gravitational 
equations thus imply the equations of motion for the matter itself. 

Finally, by comparing the static weak field limit of the Einstein equations (|2.1.14|) with the 
Poisson equation of the Newton theory of gravity, we easily get the form of the Einstein 
constant in terms of the Newton constant G as k = 8nG. 

The whole analysis discussed so far regards the Einstein geometrodynamics formula- 
tion of gravity. In the gravity-matter action, the cosmological constant term is considered 
as vanishing, although it would be allowed by the paradigm of General Relativity. This 
choice is based on the idea that such term should come out from the matter contribution 
itself, eventually on a quantum level. 



2.2 Matter Fields 



In General Relativity, continuous matter fields are described by the energy-momentum 
tensor Ty 15051 13191 . We will focus our attention only on the tensor fields and will not 
discuss more general entities as the spinor fields for which we recommend the exhaustive 
discussion in |505| and for the related bi-spinor calculus j|318|| . In particular, we briefly 
review the cases of the perfect fluid, the scalar and the electromagnetic fields. 

The definition of the symmetric tensor T;y as in Eq. (|2.1.15|) differs from the energy- 
momentum tensor obtained from the Eulero-Lagrange equations. In fact, in the Minkowski 
space-time (IR 4 , rjij) we have an alternative way to construct the energy-momentum ten- 
sor for a generic field Lagrangian. However, such tensor in general has not the symmetry 
requirement, thus we adopt the definition given as in Eq. (|2.1.15|) . Let us stress that, from 
another point of view, the equations of motion of the matter field 

V'Ty = (2.2.1) 

are obtained requiring the action for the matter S m introduced in Eq. (|2.1.11|) to be invari- 
ant under diffeomorphisms. For the Einstein equations, the Bianchi identities V'G« = 
may be viewed as a consequence of the invariance of the Hilbert action under diffeomor- 
phisms. Let us now list the principal aspects of the three fields under consideration. 

• The energy-momentum tensor of a perfect fluid is given by 

T i j = (p + p)u i U j -pg ij/ (2.2.2) 

where m, is a unit time-like vector field representing the four- velocity of the fluid. 
The scalar functions p and p are the energy density and the pressure, respectively, as 



13 



2 Fundamental Tools 



measured by an observer in a locally inertial frame moving with the fluid, and are 
related by an equation of state p = p(p). Since no terms describing heat conduction 
or viscosity are introduced the fluid is considered as perfect. For the early Universe 
we have p = (7 — l)p, where 7 is the polytropic index. In the particular case where 
p = we deal with dust, whose fluid elements follow geodesic trajectories. 



The Lagrangian density for the linear, relativistic, scalar field theory reads as 

Am = \ - m 2 f) (2-2-3) 

and in a curved space-time it is obtained by the minimal substitution rule — * gjj 
and 9, — > V;. The corresponding energy-momentum tensor is expressed as 



Tij = Vi<pVj<p - - gij (V k cp\/ k <p -m 2 cp) , (2.2.4) 



and the Klein-Gordon equation becomes 

Vv,t + m 2 ) <j> = 0. (2.2.5) 



• The Maxwell field is described by the Lagrangian density 

A m = (2.2.6) 
and the electromagnetic energy-momentum tensor reads as 

Tij = ^ (-FikF) + l §lj F kl F k ^ , (2.2.7) 

where F = dA is the curvature 2-form associated to the connection 1-form A — 
A{dx l , i.e. Fjj = VjAj — VyA,. We note that the trace of the energy-momentum 
tensor defined as in Eq. (|2.2.7|) is identically vanishing. Thus, the Maxwell equations 
in a curved space-time become 

V'F lk = -4nj k (2.2.8a) 
V[;F ;7c] =0, (2.2.8b) 

where j k is the current density four-vector of electric charge and [ ] denote the anti- 
symmetric sum, or in the forms language equations (|2.2.8|l are written as 

d*F = 4n*i (2.2.9a) 
dF = (2.2.9b) 

where d is the exterior derivative and * is the Hodge star operator |505l t3701. 
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2.3 Tetradic Formalism 



Let M be a four-dimensional manifold and e a one-to-one correspondence on it, e : M — > 
TM, mapping tensor fields on M to tensor fields in the Minkowski tangent space TM. 
The four linearly independent fields e\ (tetrads or vierbein) (a is the Lorentzian index) are 
an orthonormal basis for the local Minkowskian space-time and satisfy the condition 

e'atbi = nab (23.1) 

only. In Eq. <|2.3.1|) rj ab is a symmetric, constant matrix with signature (+, — , — , — ) and 
rj ah is the inverse matrix of rj ab . Let us define the reciprocal (dual) vectors e m , such that 
e \ e \ ~ ^ definition of e ai and by (|2.3.1|) , the condition e a ^ a = 5\ is also verified. 
In this way, the Lorentzian index is lowered and raised with the matrix \] ab , and the metric 
tensor gjj assumes the form 

gii = Vab<$Jj. (2.3.2) 

We denote the projections of the vector A' along the four e a as "vierbein components", 
A a = e' a Aj and A a = e"A' = t] ah A b . In particular, for the partial differential operator, we 
have d a = e l a di, below denoted as a comma ( ) /H = d a . The generalization to a tensor of 
any number of covariant or contra variant indices is straightforward. 

Let us introduce a couple of quantities which will directly appear in the Einstein equa- 
tions, the Ricci coefficients 'Yabc, an d their linear combinations K a \, c by the following equa- 
tions 

labc = V fc e a ;4e*, (2.3.3) 

^-abc = Ifabc Ifacb- (2.3.4) 

With some algebra, the Riemann and the Ricci tensors can be expressed in terms of ^y a \, c 

and of A„i, c as 

Rabcd = 7abc,d ~ labd,c + labf (j f cd ~ 7"^) + Ta/cT^M ~ lafd7 f bc , (2-3.5) 



Kb = ~\ (hb'c + X ba, c + A °ca,b + AC C fy 



+ h cd b ^cda + A cd b A cdfl - ^A b cd A acd + A c cd A fli , d + A c cd A b /^ 



(2.3.6) 



With the use of the tetradic fields e a -, we are able to rewrite the Lagrangian formulation of 
the Einstein theory in a more elegant and compact form. Let us introduce the connection 
1-forms col anc ^ * ne curvature 2-forms associated to it by 

R a b = dcv a b + oo a c A cv c b , (2.3.7) 

which is the so-called I Cartan structure equation. In this formalism, the action for GR in 
absence of matter field reads as 

S{e,co) = ~ f e abcd e a A e b A R cd , (2.3.8) 
4k Jm 



15 



2 Fundamental Tools 



where e a \, c & is the totally antisymmetric tensor on TM such that e^i = £ a bcd e< i el j eC k e i- We 
stress that here we are working a la Palatini, i.e. in a first order formalism where the 
tetrads e" and the spin connections co a b can be considered as independent variables. The 
variation of the action (|2.3.8[) with respect to the connections leads to the II Cartan struc- 
ture equation, i.e. to the equations of motion for <x>\ 

de a + co a b A e h = => w = w{e) (2.3.9) 

and to the identity e b A R b = 0, while the variation with respect to the tetrad e a leads to 
the Einstein vacuum equations 

e abcd e b AR cd = 0. (2.3.10) 

The solution to the II Cartan structure equation (|2.3.9|) can be written as cof = e s 'Vtfi£, 

which corresponds to the realtion cof 3 = r y ab e\ for the Ricci coefficients. The two different 
notations are adopted accordingly as 7 in ||319|| , and co in [505 J. 



2.4 Hamiltonian Formulation of the Dynamics 

The aim of the first part of this paragraph is the Hamiltonian formulation of GR in the 
metric formalism |505l 12801 14851 1310 11, and of the corresponding Hamilton- Jacoby theory 
in the following 12591 1297| . For the theory of constrained Hamiltonian systems we refer 
to the standard textbooks lfT8Tll25T1l . 



2.4.1 Canonical General Relativity 

The generally covariant system par excellence is the gravitational field, i.e. it is invari- 
ant under arbitrary changes of the space-time coordinates (four-dimensional diffeomor- 
phisms). To perform a canonical formulation of General Relativity one has to assume 
that the topology of M (the physical space-time) is M = R x Z, where E is a compact 
three-dimensional manifold (the three-space). As follows from standard theorems, all 
physically realistic space-times (globally hyperbolic) posses such topology; thus, M can 
be foliated by a one-parameter family of embeddings X f : E — > M, t £ R, of E in M. As a 
consequence, the mapping X:RxE-> M, defined by (x, t) — > X f , is a diffeomorphism 
of R x E to M. A useful parametrization of the embedding is given by the deformation 
vector field 

?)X'(r f) 

r(X) = dt = N(X)n'(X) + AT'(X) , (2.4.1) 

where N'(X) = N a X' a . The vector field T 1 , satisfying T'V,T = 1, can be interpreted as the 
"flow of time" throughout the space-time. In Eq. (|2.4.1)l . n 1 is the unit vector field normal 
to E f , i.e. the relations 

U'^ = 1 (2.4.2) 

\gijn l d a Xl = 

hold. The quantities N and N a are the lapse function and the shift vector, respectively. 
This way the space-time metric gu induces a spatial metric, i.e. a three-dimensional Rie- 
mannian metric h a p on each Ef, by — h x p = gjjdaX'dpXl . The space-time line-element 
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adapted to this foliation can be written as 

ds 2 = N 2 dt 2 - h afi (dx a + N K dt) (dx? + N?dt) . (2.4.3) 

This formallism was introduced by Arnowitt, Deser and Misner in 1962 [19J and is known 
as the ADM procedure. 

The geometrical meaning of N and N a is now clear: the lapse function specifies the proper 
time separation between the hypersurfaces X t (E) and X t+ $ t (L) measured in the direction 
normal to the first hypersurface. On the other hand, the shift vector measures the dis- 
placement of the point X t+ st(x) from the intersection of the hypersurface Xf + «(L) with 
the normal geodesic drawn from the point X f (x) (see Fig l2.1|) . In order to have a future 
directed foliation of the space-time, we stress the requirement for the lapse function N to 
be positive everywhere in the domain of definition. 




Figure 2.1: The geometric interpretation of the lapse function and of the shift vector (from 
||310I0 . In the figure latin indices denote spatial components, while greek ones space-time 
components. 



In the canonical analysis of General Relativity, the Riemannian metric h a p on Ef plays 
the role of the fundamental configuration variable. The rate of change of h x p with respect 
to the time label f is related to the extrinsic curvature K K p = — 1 /2£ Jn h lx p of the hypersurface 
E{ as 

K a f}(x, t) = -— (dtKp ~ (&N*h) K(i ) , (2.4.4) 
where L a denotes the Lie derivative along the vector field a. 

Let us now pull-back the Einstein Lagrangian density by the adopted foliation X : Ex 
E — > M and express the result X* in terms of the extrinsic curvature, the three-metric h K p, 
N and N a . This gives 

- 2j t (VhKl) - dp - h^d a N) , (2.4.5) 

where WjR is the n-dimensional curvature scalar and — g = NVh, being h the deter- 
minant of the three-dimensional metric. We are able to re-cast the original Hilbert action 
into a 3 + 1 form by simply dropping the total differential in (|2.4.5|) , i.e. the last two terms 
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on the r.h.s. of (|2.4.5|) 

S(h,N,N«) = y dtJ^d 3 xNVh (k^K^ - (K^) 2 + ( 3 >r) . (2.4.6) 

Performing a Legendre transformation of the Lagrangian density appearing in equation 
<|2.4.6|) we obtain the corresponding Hamiltonian density. Let us note that the action (|2.4.6|) 
does not depend on the time derivatives of N and N a , and therefore using the definition 
(|2.4.4|) and the fact that does not contain time derivatives, we obtain for the conjugate 
momenta 



±U^(x t t) = f . , <N =1±(k^- h^K\) (2.4.7a) 



k ' Sh K p(x,t) k 

TI(x, t) = — = (2.4.7b) 

V ; SN(x,t) V ' 

n a (x, t) = . 3S - = . (2.4.7c) 

From equations (|2.4.7|l follows that not all conjugate momenta are independent, i.e. one 
cannot solve for all velocities as functions of coordinates and momenta: one can express 
h a a in terms of h a p, N, N K and LT^, but the same is not possible for N and N a . By other 
words, we have the so-called primary constraints 

c(x,t) = n(x,t) = 0, c a (x,t) = n K (x,t) = o, (2.4.8) 

where "primary" emphasizes that the equations of motion are not used to obtain relations 
d2A8l . 

According to the theory of constrained Hamiltonian systems, let us introduce the new 
fields A(x, t) and A a (x, t) as the Lagrange multipliers for the primary constraints, making 
the Legendre transformations invertible and the corresponding action reading as 

S = J dt J^d 3 x h^U^ + M + N x TI a - (AC + A a Q + N a H a + NH) , (2.4.9) 

where 



H a = -IKyVpW? (2.4.10a) 

h = g^n^n " 5 - Vh {3) r (2.4.10b) 

in which 

1 

9ap 7 6 = — ^ (h a7 hps + hp y h aS - h a ph 7 s) (2.4.10c) 

is the so-called super-metric on the space of the three-metrics, and the functionals H a and 
H are the super-momentum and super-Hamiltonian, respectively. 

The classical canonical algebra of the system is expressed in terms of the standard Pois- 
son brackets as 
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{h a p(x,t),h 7 s(x',t)} = (2.4.11a) 
{ll^(x r t),W s {x',t)} =0 (2.4.11b) 
{h 7S (x f t),U a ^x',t)} =K8l/ s) 8\x-x'), (2.4.11c) 

where the parentheses ( ) denote symmetrized indices. We can analyze the meaning of 
the Hamiltonian of the system, bracketed in equation (|2.4.9[) 

% = J^d 3 x (AC + A a C a + N X H X + NH) = C(A) + C (T) + H (tfj + H(N). (2.4.12) 

The variation of the action (|2.4.9|) with respect to the Lagrange multipliers A and A a re- 
produces the primary constraints (|2.4.8|l . The consistency of the dynamics is ensured by 
preserving them during the evolution of the system, i.e. requiring 

C(x,t) = {C(x,t),0i} = 0, C a {x,t) = {C a (x,t),%} = 0. (2.4.13) 

The Poisson brackets in equation (|2.4.13|) do not vanish, but equal to H(x, t) and H a (x, t), 
respectively, and therefore the consistency of the motion leads to the secondary constraints 

H{x, t) = 0, H a (x, 0=0, (2.4.14) 

by means of the equations of motion. Let us observe that the Hamiltonian in General 
Relativity is constrained as "K ~ being weakly zero, i.e. vanishing on the constraint 
surface. This is not surprising since we are dealing with a generally covariant system. 

A problem that can arise is that the constraint surface, i.e. the surface where the con- 
straints hold, could not be preserved under the motion generated by the constraints them- 
selves, but this is not the case. In fact, the Poisson algebra of the super-momentum and 
super-Hamiltonian, computed using (|2.4.11a|l - (|2.4.11c|) , is closed. In other words, the set 
of constraints is of first class, i.e. the Poisson brackets between the Hamiltonian and any 
constraint weakly vanish, as arising from the relations 

{X,H(f)} = - H(L f N) (2.4.15) 

{X,H(f)} = H(L^f) + H(N(N,f,h)), (2.4.16) 

where / is a smooth function and N a (N,f,h) = h a P(Ndpf — fdpN). These equations are 
equivalent to the Dirac algebra ||181|| . 

Moreover, varying the action (|2.4.9|) with respect to the two conjugate momenta n and 
IF, we obtain 

N(x,t) = X{x,t), N a (x,t) = A x {x,t) (2.4.17) 

ensuring that the trajectories of the lapse function and of the shift vector in the phase 
space are completely arbitrary. 

The Hamiltonian of the theory is not a true one but a linear combination of constraints. 
From relations (|2.4.15|) and (|2.4.16|) it is possible to show that rather than generating time 
translations, the Hamiltonian generates space-time diffeomorphisms, whose parameters 
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are the completely arbitrary functions N and N tt , and the corresponding motions on the 
phase space have to be regarded as gauge transformations. An observable is defined as a 
function on the constraint surface that is gauge invariant, and more precisely, in a system 
with first class constraints it can be described as a phase space function that has weakly 
vanishing Poisson brackets with the constraints. In our case, A is an observable if and 
only if 

{A, JC(A, A a , N tt ,N)} 0, (2.4.18) 

for generic A, A a , N a and N. By this definition, one treats on the same footing ordinary 
gauge invariant quantities and constants of motion with respect to evolution along the 
foliation associated to N and N K . The basic variables of the theory, h a p and Tl K P, are not 
observables as they are not gauge invariant. In fact, no observables for General Relativity 
are known, except for the particular situation with asymptotically flat boundary condi- 
tions. 

Let us remark that the equations of motion 

h a *(x,t) = tl*Hx,t) = — rt~~7~ x (2-4.19) 

together with the eight constraints (|2.4.8|) and (|2.4.14|) are completely equivalent to the 
vacuum Einstein equations, Ru = |199| . 

2.4.2 Hamilton-Jacobi Equations for Gravitational Field 

The formulation of the Hamilton-Jacobi theory for a covariant system is simpler than the 
conventional non-relativistic version. In fact, in such case the Hamilton-Jacoby equations 
are expressed as 

H(*.^)-0, (2.4.20) 

where H and S(q a ) denote the Hamiltonian and the Hamilton functions, respectively. 
For GR, the Hamilton-Jacobi equations arising from the super-Hamiltonian and super- 
momentum read as 

HJS = % h s f* ,f " ®R = (2.4.21) 
HJ K S = -2/i KT V«— - = . (2.4.22) 

These four equations, together with the primary constraints (|2.4.8[) , completely define the 
classical dynamics of the theory. 

Let us point out how through a change of variable we can define a time coordinate. 
In fact, writing h a p = f^'^U^, with rj = h 1 ^ and detw a ^ = 1, the super-Hamiltonian 
Hamilton-Jacobi equation (|2.4.21|) can be written as 



^ J 5 = -T6 K (I)' + % U « U P>^jL ~ h 1 '^^,^), (2.4.23) 



ZK 01 
H — jUuyUfiS— 

where the potential term V comes out from the spatial Ricci scalar and V refers to spatial 
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gradients only. As we can see from Eq. (|2.4.23|) , tj has the correct signature for an internal 
time variable candidate. We will show later how this variable is nothing but a power of 
the isotropic volume of the Universe. 



2.5 Synchronous Reference 

In this paragraph we will focus our attention on one of the most interesting reference 
system, i.e. the synchronous one. For a detailed discussion see [319 J. It is defined by the 
following choice for the metric tensor 

goo = 1, goa = 0, (2.5.1) 

thus in the 3 + 1 framework we have to require N = 1 and N a = in (|2.4.3|) . The first 
condition in (|2.5.1|) is allowed from the freedom to rescale the variable t with the transfor- 
mation y/goodt, in order to reduce goo to unity and setting the time coordinate x° = t as 
the proper time at each point of space. The second one is allowed by the non-vanishing of 
\g a p | and allows the synchronization of clocks at different points of space. The elementary 
line interval is given by the expression 

ds 2 = dt 2 - dl 2 , (2.5.2) 

where 

dl 2 = h Kfi (f, x" 1 ) dx a dx£ , (2.5.3) 

in which the three-dimensional tensor h a p defines the space metric. 
In such a reference system, lines of equal times are geodesies in the four-space, as implied 
by the splitting definition. Indeed the four- vector u 1 = dx'/ds which is tangent to the 
world line (x 7 = const.), has components U° = 1, u a = and automatically satisfies the 
geodesic equations 

du' 

— +I>V = r 00 = 0. (2.5.4) 

The choice of such a reference is always possible and moreover it is not unique. In fact, 
considering a generic infinitesimal displacement 

t' = t + Z(t,xP), x' a = x a + £ a (t,xP) (2.5.5) 

and the associated four-metric change g'^ = gjj — 2V(j£y), the conditions of preserving the 
synchronous reference can be written as 

d t Z = t' = t + £(x p ) (2.5.6a) 
9 K £ = x' a = x a + J h^dt + (p a (xP) , (2.5.6b) 

(p a being generic space functions. 

In the reference defined by the metric (I2.5.2D , the Einstein equations are written in mixed 
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components as 



R°o = -~< ~ \4*$ = BnG (V ° - It) (2.5.7a) 
< = \ (</3 - *?J = BttGIJ (2.5.7b) 
= -?£ - (Vfocf) = 8tiG f if - ^r) , (2.5.7c) 



where 



k^ = ^, h=\Kp\. (2-5.8) 

P^p is the three-dimensional Ricci tensor obtained through the metric h^p which is used 
to raise and lower indices within the spatial sections. 

The metric allows to construct the three-dimensional Ricci tensor pf = h^P^y as 

fy» = Mlp - ^ + - ( 2 - 5 - 9 ) 

in which appear the pure spatial Christoffel symbols 

AJ, = \w s {d a h 5fi + 3^^ - 9^) . (2.5.10) 

From ( |2.5.7aD it is straightforward to derive, even in the isotropic case, the Landau-Ray chaudhuri 
theorem, stating that the metric determinant h must monotonically vanish in a finite in- 
stant of time. However, we want to stress that the singularity in this reference system is 
not physical and of course can be removed passing to another one. 



2.6 Singularity Theorems 

In this Section we investigate space-time singularities. We will show how singularities 
are true, generic features of the Einstein theory of gravity and how they arise under cer- 
tain, quite general, assumptions. We only enter in some details and we will not give the 
rigorous proofs, for which we refer to the standard literature I247l l505| . 
In particular, in this Section only we follow the signature of [505J (—,+,+,+)• 
We will treat here the classical aspects, referring to [92, 217J and references therein for a 
discussion on quantum cosmological singularities. 

After defining what a space-time singularity means we will present some basic technol- 
ogy, and finally we will discuss the theorems. 



2.6.1 Definition of a Space-time Singularity 

Let us clarify the meaning of a singularity of the space-time. In analogy with field theory, 
we can represent such a singularity as a "place" of the space-time where the curvature 
diverges, or where some similar pathological behavior of the geometric invariants takes 
place. Therefore, a first problem arises when characterizing the singularity as a "place". 
In fact, in General Relativity the space-time consists of a manifold M and a metric gj; 
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defined everywhere on M: a singularity (as the Big-Bang singularity of the isotropic cos- 
mological solution or the r = singularity in the Schwarzschild space-time) can not be 
considered as a part of the manifold itself. We can speak of a physical event only when 
a manifold and a metric structure are defined around it. A priori, it is possible to add 
points to the manifold in order to describe the singularity as a real place (as a boundary 
of the manifold). But apart from very peculiar cases, no general notion or definition of a 
singularity boundary exists |505 , 2081 1209112851 - Another problem is that singularities are 
not always accompanied by unbounded curvature as in the best known cases. There are 
several examples |247)| of a singularity without diverging curvature. In fact, as we will 
see, this feature is not the basic mechanism behind singularity theorems. 

The best way to clarify what a singularity means is the geodesic incompleteness, i.e. 
the existence of geodesies which are inextensible at least in one direction, but have only a 
finite range for the affine parameter. We can thus define a singular space-time as the one 
possessing one incomplete geodesic curve. 



2.6.2 Preliminary notions 

We initially define the notions of expansion and contraction of the congruence of time- 
like geodesies. The treatment of null geodesic congruences is conceptually similar and 
not reviewed here. 

Let O be an open set of a manifold M. A congruence in O is defined as a family of curves 
such that only one curve of this family passes through each point p £ O. Let be the 
vector field tangent to the geodesies such that = — 1, and define the tensor By = V/f, 
to introduce the expansion 9 (as well as its trace part), the shear o~u (as its symmetric, trace- 
free part) and the twist tOy (as its anti-symmetric part) of the congruence. Thus, By can be 
written as 

Bij = -9hij + <nj + coij , (2.6.1) 

where the spatial metric is hy = gy + £j£y. Given a deviation vector orthogonal to the 
vector field the tensor measures its failure from being parallelly transported along 
a geodesic in the congruence for this subfamily. 

The rates of change of 9, o~y and COy follow from the geodesic equation 

Z k V k Bij = g*V k V^i = g*(V y V^j + R kji = 

= - (Vyg*) (V,t&) + % = ~BjB ik + % . (2.6.2) 

Taking the trace of Eq. H2.6.2D , we obtain the Raychaudhuri equation 

e = g*v t = -\&- ^ + *»a<ji - RwSV • (2-6-3) 

This result will be fundamental when proving the singularity theorems. Let us focus 
our attention on the right-hand side of it: using Einstein equations, the last term can be 
written as 

Ru?? = K (jij - l -T gi ^j ?tf = K (t$$ + ir) . (2.6.4) 
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Let us state the physical criterion preventing the stresses of matter from becoming too 
large so that the right-hand side of (|2.6.4|) is not negative, obtaining 

Tij?& > -~T, (2.6.5) 

which is known as the strong energy condition B247|l . It is commonly accepted that every 
reasonable kind of matter would satisfy the condition (|2.6.5[) . Therefore, from the Ray- 
chaudhuri equation <|2.6.3)> one can see that, if the congruence is non-rotating 0, which 
means cou = 0, and the strong energy condition holds, 9 always decreases along the 
geodesies. More precisely, we get 9 + ^9 2 < 0, whose integral implies 

e-\T)>e Q l + ^T, (2.6.6) 

where 9q is the initial value of 9. For negative values of 9q (i.e the congruence is initially 
converging), 9 will diverge after a proper time no larger than t < 3/ \9q |. By other words, 
the geodesies must intersect before such instant and form a caustic (a focal point). Of 
course, the singularity of 9 is nothing but a singularity in the congruence and not a space- 
time one, since the smooth manifold is well-defined on caustics. 

Let us briefly discuss the meaning of the strong energy condition in the simple case of 
a perfect fluid, for which it reads as 

P + X>«>0, p + p a >0, (2.6.7) 

CL 

and is satisfied for p > and for negative pressure components smaller in magnitude 
then p. 

We need to introduce some notions of differential geometry and topology to translate 
the occurrence of caustics into space-time singularities. 

Let 7 be a geodesic with tangent v± defined on a manifold M. A solution of the geodesic 
deviation equation is rf satisfying □ (|2.1.6|) that constitutes a Jacobi field on 7. If it is non- 
vanishing along 7, but is zero at both p, q G 7, then p, q are said to be conjugate. 
It is possible to show that a point q G 7 lying in the future of p G 7 is conjugate to p if 
and only if the expansion of all the time-like geodesies congruence passing through p ap- 
proaches —00 at q (i.e. a point is conjugate if and only if it is a caustic of such congruence). 
A necessary hypothesis in this statement is that the space-time manifold (M,gij) satisfy 
Rij£'& > 0, for all time-like cf . Moreover, the necessary and sufficient condition for a 
timelike curve 7, connecting p,q G M, that locally maximizes the proper time between p 
and q, is that 7 is a geodesic without any point conjugate to p between p and q. 

An analogous analysis can be made for time-like geodesies and a smooth space-like hy- 
persurface E. In particular, let 9 = K = WKu be the expansion of the geodesic congruence 
orthogonal to E, Kjj being the extrinsic curvature of E. Then, for K < at the point q G E, 
within a proper time T < 3/ |X| a point p conjugate to E along the geodesic orthogonal to 
E exists, for a space-time (M,gij) satisfying Rjj^'^ > 0. As above, a time-like curve that 
locally maximizes the proper time between p and E has to be a geodesic orthogonal to E 

1 It is possible show that the congruence is orthogonal to the hypersurface if and only if = 
2 In this case, a minus sign appears in the right-hand side of <2.1.6t due to the different signature of the 
metric. 
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without conjugate point to E. 

The last step toward the singularity theorems is to prove the existence of maximum 
length curves in globally hyperbolic space-times. We recall that this is the case if they 
posses Cauchy surfaces in accordance with the determinism of classical Physics. Without 
entering in the details, in that case a curve 7 for which T attains its maximum value exists, 
and a necessary condition is that 7 be a geodesic without conjugate point. 

2.6.3 Singularity Theorems 

In the previous subsection we have summarized the necessary concepts to analyze the 
singularity theorems in some details, although we will discuss them without any proofs. 

Let a space-time manifold be globally hyperbolic, with Ry^ 1 ^ > for all the time-like 
Suppose that the trace K of the extrinsic curvature of a Cauchy surface everywhere 
satisfies K < C < 0, for a constant C. Therefore, no past-directed time-like curves A from 
E can have a length greater than 3/|C|. In fact, if there were any, a maximum length 
curve would exist and would be a geodesic, thus contradicting the fact that no conjugate 
points exist between E and p £ A: therefore such curve cannot exist. In particular, all 
past-directed time-like geodesies are incomplete. 

This theorem is valid in a cosmological context and expresses that, if the Universe is ex- 
panding everywhere at a certain instant of time, then it must have begun with a singular 
state at a finite time in the past. 

It is also possible to show that the previous theorem remains valid also relaxing the 
hypothesis that the Universe is globally hyperbolic. The price to pay is the assumption 
that E be a compact manifold (dealing with a closed Universe) and especially that only 
one incomplete geodesic is predicted. 

For the theorems proving null geodesic incompleteness in a gravitational collapse con- 
text, i.e. the existence of a singularity in a black hole space-time, we refer to the literature 
E07t 

We conclude stating the most general theorem, which entirely eliminates the assump- 
tions of a Universe expanding everywhere and the global hyperbolicity of the manifold 
(M,gij). On the other hand, we loose any information about the nature of one incomplete 
geodesic at least, since it does not distinguish between a time-like and a null geodesic. 
The space-time is singular under the following hypotheses: 

i) the condition RijV l vi > holds for all time-like and null v l , 

ii) no closed time-like curves exist and 

Hi) at least one of the following properties holds: 

a) {M,gjj) is a closed Universe, 

b) {M,gjj) possesses a trapped surfac^E 

c) there exists a point p € M such that the expansion 6 of the future or past di- 
rected null geodesies emanating from it becomes negative along each geodesic 
in this congruence. 

3 A trapped surface is a compact smooth space-like manifold T, such that the expansion 6 of either outgoing 
either ingoing future directed null geodesies is everywhere negative. 
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In particular, our Universe must be singular. In fact, the conditions i)-ii) hold and 6 for 
the past-directed null geodesies emanating from us at the present time becomes negative 
before the decoupling time, i.e. the time up to when the Universe is well described by the 
Friedmann-Robertson-Walker (FRW) model. 

The occurrence of a space-time singularity undoubtedly represents a breakdown of the 
classical theory of gravity, i.e. the General Relativity. The removal of such singularities 
is a prerequisite for any fundamental theory, like the quantum theory of the gravitational 
field. Singularity theorems are very powerful instruments, although do not provide any 
information about the nature of the predicted singularity. Unfortunately, we do not have 
a general classification of singularities, i.e. many different types exist and the unbounded 
curvature is not the basic mechanism behind such theorems. 

Let us conclude remarking how singularity theorems concern properties of differen- 
tial geometry and topology. Einstein equations are used only with respect to the positive 
curvature case, being only used in the Raychaudhuri equation. Nevertheless, no gen- 
eral conditions for non-singular solutions are known, and therefore it is not possible to 
disregard the singularity in General Relativity. 
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In this Section we introduce the homogeneous cosmological models, in order to describe 
the dynamics of the Universe toward the cosmological singularity with a more realistic 
approach then the one offered by the homogeneous and isotropic Friedmann-Robertson- 
Walker one. We will summarize the derivation of the Bianchi classification, with partic- 
ular attention to the dynamics of the Bianchi types II and VII, while we will dedicate the 
entire Section U to types VIII and IX, the so-called Mixmaster Il54ll367ll368| . 

3.1 Homogeneous Spaces 

The study of homogeneous models arises from breaking the hypothesis of space isotropy. 
We will focus in particular on space-times spatially homogeneous, without treating the 
non physical case of a space-time homogeneous manifold, i.e. where the metric is the 
same at all points of space and time, because it represents a Universe not expanding at 
all. We will follow in particular the description as in Refs. |505H319ll446ll473| . 

When relaxing the space isotropy hypothesis we gain a significantly larger arbitrariness 
in the solution. In particular, the homogeneous and isotropic model possesses a single 
gravitational degree of freedom given by the scale factor leaving free only the curvature 
sign BIO) . This larger class of solutions has still a finite number of degrees of freedom 
and very general properties characterize the dynamics during the evolution toward the 
initial singularity [54J . 

A spatially homogeneous space-time is defined as a manifold with a group of isome- 
tries, i.e. a group of transformations leaving invariant the metric g. The Killing vectors 
£ are the corresponding infinitesimal generators, with vanishing Lie derivative L?g = 
and whose orbits are the space-like hypersurfaces which foliate the space-time. We will 
discuss in which sense the metric properties are the same in all space points under homo- 
geneity. 

3.1.1 Killing Vector Fields 

The Lie algebras of Killing vector fields generate the groups of motions via infinitesimal 
displacements, yielding conserved quantities and allowing a classification of homoge- 
neous spaces. Consider a group of transformations 

x f -> x* =f (x,a) (3.1.1) 

on a space M (eventually a manifold), where {a a } « are r independent variables which 
parametrize the group and let a.Q correspond to the identity 

F (x,ao) = X* . (3.1.2) 
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Take an infinitesimal transformation «o + $ a r i.e. one which is very close to the identity so 
that 



& = f (x,a + Sa) « f (x,a ) + [ ) (z,«o)<W (3.1.3) 



sS2W 

i.e. 

x f< ^ ^ « + ^ (x) <Sfl fl = (1 + ^ fl £ fl ) , (3.1.4) 

where the r first-order differential operators {£„} are defined by £ fl = jot and corre- 
spond to the r vector fields with components {£«}■ These are the "generating vector 
fields" and when the group is a group of motions, they are called Killing vector fields, sat- 
isfying Hgg = 0, so that, under infinitesimal transformations (|3.1.3|) all points of the space 
M are translated by a distance Sx^ = Sa a ^ m the coordinates {x^} and 

ta (1 + &* fl £ fl ) x* « e*'" f «jc' < . (3.1.5) 

In fact, the finite transformations of the group may be represented as 

Xf -» ^ = e 6 "^' (3.1.6) 

where {0"} are r new parameters on the group. 

These generators form a Lie algebra[446], i.e. a real r-dimensional vector space with basis 
{£„}, which is closed under commutation, i.e. the commutators of the basis elements can 
be expressed as constant linear combinations of themselves 

[£,, ft] = £ fl ft - ftft = ±C e a6 & (3.1.7) 

where C c nfc are the structure constants of the Lie algebra ((+) refers to left-invariant groups, 
while ( — ) to right-invariant ones). Suppose {e a } is a basis of the Lie algebra g of a group 
G 

[e a ,e b ] = C c ab e c (3.1.8) 

and define 

lab = C c ad C d bc = lba (3.1.9) 
which is symmetric by definition, providing a natural inner product on the Lie algebra 

lab = e a -e b = j(e a ,e b ) ; (3.1.10) 

when det (j ab ) ^ 0, Eq. (13.1.10ft is non-degenerate and the groups for which this is true 
are called semi-simple. The r vector fields {e a } may be used instead of the coordinate 

basis | g^r | as a basis in which to express an arbitrary vector field on G. This basis is a 

frame. 

3.1.2 Definition of Homogeneity 

Suppose that the group acts on a manifold M as a group of transformations 

x* ^ f" (x,a) = (x) (3.1.11) 
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and let us define the orbit of x 

fc (x) = {/« (x) | a G G} (3.1.12) 

as the set of all points that can be reached from x under the group of transformations. 
Thus, the group of isotropy at x is 

G x = {a G G | /„ (x) = x} (3.1.13) 

i.e., it is the subgroup of G which leaves x fixed. Suppose G x = {flo} an d fc { x ) = M, i.e. 
every transformation of G moves the point x, and every point in M can be reached from 
x by a unique transformation. Since G|G X = {aao \a G G} = G, G is diffeomorphic to M 
and one may identify the two spaces. 

If g is a metric on M invariant under G, it corresponds to a left-invariant one on G, spec- 
ified entirely by the inner products of the basis left-invariant vectors fields e a . For three 
dimensions one obtains the family of spatially homogeneous space sections of the spa- 
tially homogeneous space-times. 

Given a basis {e a } of the Lie algebra of a three dimensional Lie group G, with structure 
constants C c ab , the spatial metric at each moment of time is specified by the spatially con- 
stant inner products 

e a -eb = gab (0 , (3.1.14) 

which are six functions of time. The Einstein equations, as we will apply, become ordi- 
nary differential equations for these six functions, plus whatever functions of time are 
necessary to describe the matter of the Universe. 

For both homogeneous and spatially homogeneous space-times, one needs only to 
consider a representative group from each equivalence class of isomorphic Lie groups 
of dimension four and three, respectively. In three dimensions the classification of in- 
equivalent three-dimensional Lie groups is called the Bianchi classification [80J and de- 
termines the various symmetry types possible for homogeneous three spaces, just as 
(k = +1,0,-1) classify the possible symmetry types for homogeneous and isotropic 
three-spaces (FRW). 

After obtaining all the three-dimensional Lie groups according to the Bianchi classifica- 
tion we will write down and discuss the corresponding Einstein equations. 



3.1.3 Application to Cosmology 

A homogeneous space-time is defined by space-like hypersurfaces E such that for any 
points 6 E there is a unique element T G G such that r(p) = q (in this case the Lie 
group acts simply transitively on each E). Such uniqueness implies dimG = dimE = 3, 
and G and E can be identified (for example, in the simplest case of translations group we 
have G = R 3 ); thus, the action of the isometries on E is just the left multiplication on 
G and tensor fields invariant under isometries are the left-invariant ones on G. In four 
dimensions one obtains the homogeneous space-times and the foliation of M turns out to 
be M = R x G. In order to preserve metric properties at all points, let us consider the 
group of transformations of coordinates which transform the space into itself, leaving the 
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metric unchanged: if the line element before the transformation has the form 

dl 2 = h a p (V, x 2 , x 3 ) dx^dx? , (3.1.15) 

then it becomes 

dl 2 = h afi (x n , x a , x' 3 ) dx' a dx lfi , (3.1.16) 

where h a p has the same form in the new coordinates. 

In the general case of a non Euclidean homogeneous three-dimensional space, there are 
three independent differential forms which are invariant under the transformations of the 
group of motions, however they do not represent the total differential of any function of 
the coordinates. We shall write them as to" = e a a dx a . Hence the metric (|3.1.16|) is re- 
expressed as 

dl z = nah {^dx a )(e\dx^ (3.1.17) 



so that the metric tensor reads as 

where rj a j, is a function of time only, symmetric in ab and in contravariant components we 
have 

h a P = tj ab (ty a (x^ b (xi), (3.1.19) 

where r] ab should be viewed as the components of the inverse matrix. All considerations 
developed in the previous Sections here apply straightforwardly. 

The relationship between the covariant and contravariant expression for the three basis 
vectors is 



1 



v 



1 



ei = - [e 2 A e 3 ] , e 2 = - e 3 A e 1 , e 3 = - e 1 A e 2 , (3.1.20) 



V L 



where e" and e a are to be understood formally as Cartesian vectors with components e a a 
and e", while v represents 

v =\e a K \=e l ■ [e 2 Ae 3 ] . (3.1.21) 

The determinant of the metric tensor d3.1.18b is given by 7 = tjv 2 where tj is the determi- 
nant of the matrix 77^ . 

The invariance of the differential form (|3.1.16b means that 

e" a (x) dx a = e a tt (x r ) dx' tt (3.1.22) 

and e a K on both sides of (|3.1.22|) are the same functions expressed in terms of the old and 
the new coordinates, respectively. 

The algebra for the differential forms permits to rewrite (|3.1.22|) as 

^ = 4 (x') C (x) . (3.1.23) 

This is a system of differential equations which define the change of coordinates x'P(x) in 
terms of given basis vectors. 
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Integrability over the system (|3.1.23[) is rewritten in terms of the Schwartz condition 

d 2 x'P 



dx^dx'r dx^dx a 



(3.1.24) 



which, explicitly, leads to 



delix^s (xf] _ d 4 ( X ') C S (x n 



4 (*) < (x) 



£ (x'\ 



Multiplying both sides of l|3.1.25|l by e^{x)ej{x)er^{x') and differentiating, the left-hand 
side becomes 



de a 7 (x) de a K (*) 



dx tt dx^ 



(3.1.25) 



J ( v i 



<£ (x 1 ) 4 (X>) 



?>4 (*') H oo 



dx' s 



dx'y 



(3.1.26) 



and the right-hand side gives the same expression but in terms of x. 

Since x and x' are arbitrary, both sides must be constant, and (|3.1.26|) reduces to 



c 



dxP dx a 



eX = C c ab , (3.1.27) 



which gives the relations of the vectors with the group structure constants C c ab . Multiply- 
ing (|3.1.27|) by el, we finally have 

e ^~ e ^- C ^c (3-1.28) 

Similarly, such expression in the forms language is given by the left-invariant 1-form 
co" satisfying the Maurer-Cartan equation 

dco a = -Cl c co b A co c . (3.1.29) 

By construction, we have the antisymmetry property from (|3.1.26|) or (|3.1.27|) 

C ah = -C ba . (3.1.30) 

Defining 

y_ — p r 

' dx a ' 



*a = el— K , (3.1.31) 



equation (|3.1.28|) rewrites as 

[X a ,X b ] =C c ab X c . (3.1.32) 
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Homogeneity is expressed as the Jacobi identity 

[[X a ,X b ],X c ] + [[X fc ,X c ],X fl ] + [[X c ,X a ],X b ] = (3.1.33) 

and explicitly 

C f ab C d cf + C f bc C d af + d ca C% = . (3.1.34) 

With this formalism, the Einstein equations for a homogeneous Universe can be written as 
a system of ordinary differential equations which involve only functions of time, provided 
all three-dimensional vectors and tensors are projected on the tetradic basis, while the 
explicit coordinate dependence of the basis vectors is not necessary for the equations 
ruling the dynamics. In fact, such choice is not unique as e a = A ba e b yields again a set of 
basis vectors. 

Introducing the two-index structure constants as C c ab = £ a bdC dc , where e abc = e abc is the 
Levi-Civita tensor (£123 = +1), the Jacobi identity (|3.1.34|) becomes 

e bcd C cd C ba = 0. (3.1.35) 

The problem of classifying all homogeneous spaces reduces to finding all inequivalent 
sets of structure constants. 



3.2 Bianchi Classification and Line Element 

The list of all three-dimensional Lie algebras was first accomplished by Bianchi [81 J and 
each algebra uniquely determines the local proprieties of a three-dimensional group. If 
a space is homogeneous and its Lie group is the "Bianchi Type N" (N=I,...,IX), the sub- 
classification of the Bianchi groups agrees with the one made by Ellis and MacCallum 
ITT851 14461 13191 . 
Any structure constant can be written also as 

C a bc = e bcd m da + S a c a b -S a b a c , (3.2.1) 

where the matrix m ab = m ba . The subclassification as class A and class B models refers to 
the cases a c = or a c 7^ 0, respectively. 

The Jacobi identity (|3.1.34[) written for the structure constants like in (|3.2.1|) reduces to 
the condition 

m ab a b = 0. (3.2.2) 

Without loss of generality, we can put a c = (a, 0, 0) and the matrix m ab can be described 
by its principal values n\,n^,n^. The condition (|3.2.2|) becomes nothing but an\ = 0, i.e. 
either a or n\ has to vanish. Condition (|3.1.35|) rewrites explicitly as 

[Xx,X 2 ] = -aX 2 + n 3 X 3 

[X 2 ,X 3 ] = n x X x (3.2.3) 
[X 3 ,Xx] = n 2 X 2 + aX 3 

where a > 0, {ti\, n 2 , n 3 ) can be rescaled to unity, and we finally get the Bianchi classifica- 
tion as in the Table 13.1 1 
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Let us note that the Bianchi type I is isomorphic to the three-dimensional translation 
group R 3 , for which the flat FRW model is a particular case (once isotropy is restored) 
and analogously the Bianchi type V contains, as a particular case, the open FRW. 

Table 3.1: Inequivalent structure constants corresponding to the Bianchi classification. 



Type 


a 


ni 


n 2 


n 3 


I 














II 





1 








VII 





1 


1 





VI 





1 


-1 





IX 





1 


1 


1 


VIII 





1 


1 


-1 


V 


1 











IV 


1 








1 


VII 


a 





1 


1 


III (a = 1) ) 
VI (a ^ 1) j 


a 





1 


-1 



Not all anisotropic dynamics are compatible with a satisfactory Standard Cosmological 
Model but, as shown in the early Seventies, some can be represented, under suitable 
conditions, as a FRW model plus a gravitational waves packet ||341l f225l. 
The interest in the Mixmaster 113681 [3671 relies on having invariant geometry under the 
SO (3) group, shared with the closed FRW Universe. From a cosmological point of view, 
the relevance of this model arises also from the decomposition of the line element as 

ds 2 = ds 2 - S {a){b) G^ m dx i dx k (3.2.4) 

where ds^ denotes the line element of an isotropic Universe having positive constant cur- 
vature, Gfc'^' is a set of spatial tensors and S^^\(t) are amplitude functions, resulting 
small sufficiently far from the singularity. The tensors introduced in (|3.2.4|) satisfy the 
equations 

G («)(6);/ = _ (w 2_ 3)g WW G («)(2>)fc = 0/ G W(»)i = 0/ (3.2.5) 

in which the Laplacian is referred to the geometry of the sphere of unit radius. 

Let us choose a basis of dual vector fields co", preserved under isometries. Therefore, 
recalling (|3.1.19|) , the four-dimensional line element can be expressed as 

ds 2 = N 2 (t)dt 2 - rj ah (f ) co a ® co\ (3.2.6) 

parametrized by the proper time, where co a obey the Maurer-Cartan equations (|3.1.29[) . 
The explicit expression for the co a for Bianchi I, whose structure constants are Ci = 0, 

is 

co 1 = dx 1 , co 2 = dx 2 , to 3 = dx 3 , (3.2.7) 
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while for Bianchi IX, being C| c = e abc , will be specified in (|4.7.16b|) . 

In order to distinguish between expansion (change of volume) and anisotropy (change 
of shape) it is useful to parametrize the metric on the spatial slices as 

nab = R 2 e 2a ^) ab , (3.2.8) 

where Ro is the initial radius of the Universe and all the other parameters are functions 
of time only. The matrix /3 a j, satisfies the condition Tr/3 = 0, ensuring hypersurface three- 
volume dependance on the conformal factor a only as Vuniv ~ R(je 3a ||370fl . 

For a review of works devoted to the study of the dynamics of Bianchi models in differ- 
ent cosmological paradigms see the collection of articles If3^ 12001 144TI [1161 13481 14831 
ITMI4TI 15001 TO MM 13451 TO 
l3l3iraMl35^IT37lTOl2^l^ 



3.3 Field Equations 

As we derived, the Einstein equations for a homogeneous Universe can be written in the 
form of a system of ordinary differential equations which involve functions of time only, 
once introduced the tetradic basis built as in the previous Section IZ3l whose projections 
in empty space take the form 



R° 

K Q — 


1 1 

x j,a 1 b v a 
2~ K a ~^ K a K b 


(3.3.1a) 


R° = 




(3.3.1b) 


7?" — 
K b — 


-2^*<VW)-S. 


(3.3.1c) 


where the relations K ab = fj a }, and 


= rjact] cb hold, the dot denotes differentiation with 



respect to t, and the projection P fl /, = f]b c Pa °f the three-dimensional Ricci tensor becomes 

Pub — —-^{c^Ccda + C cd h Cdca ~ ^b^^acd + ^cd^ab + ^cd^ba^j • (3.3.2) 

The Einstein equations have reduced to a much simpler differential system, involving 
only ordinary derivatives with respect to the temporal variable t. 

In the following we will discuss the Kasner solution which will be generalized to the 
dynamics of Bianchi types VIII and IX. 

3.4 Kasner Solution 

The simplest and paradigmatic solution of the Einstein equations (|3.3.1a|) - (|3.3.1c|) in the 
framework of the Bianchi classification is the type I model, firstly obtained by Kasner[289] 
in 1921, which is appropriate to describe the gravitational field in empty space. 
The simultaneous vanishing of the three structure constants and of the parameter a im- 
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plies the vanishing of the three-dimensional Ricci tensor as well 

^P ab = 0. (3.4.1) 



e a = 5 a 

C c = 



Furthermore, since the three-dimensional metric tensor does not depend on space coordi- 
nates, also Ro a = 0, i.e. this model contains the standard Euclidean space as a particular 
case. Then the system (|3.3.1a|H|3.3.1c[) describes a uniform space and reduces to 



~k\A = (3.4.2a) 



1 

-L^(^) = 0. (3.4.2b) 

From (|3.4.2b|) we get the first integral 

y/ri4 =2A b a = const, (3.4.3) 
and contraction of indices a and b leads to 

K a a = l = ^-k a ai (3.4.4) 

and finally 

rj = Gt 2 G = const. . (3.4.5) 

Without loss of generality, a simple rescaling of the coordinates x a allows to put such a 
constant equal to unity, thus providing 

A» = l. (3.4.6) 

Substituting (|3.4.3[) in (|3.4.2a|) one obtains the relations among the constants \ b a 

\%k\ = \, (3.4.7) 

and lowering index b in (|3.4.3|l one gets a system of ordinary differential equations with 
respect to j ab 

2 

Vab = jKVcb ■ (3.4.8) 

The set of coefficients X c a can be considered as the matrix of a certain linear transformation, 
reducible to its diagonal form. In such a case, denoting its eigenvalues as p\, p2, P3 real 
and not equal to each other, and its eigenvectors as nW, n( 2 ), n' 3 ), the solution of (13.4.8I ) 
writes as 

nah = fnn^nW + t^n^n^ + ^ n f) n f) . (3.4.9) 

If we choose the frame of the eigenvectors (recall that e a a = S%) and denote them with 
x 1 , x 2 , x 3 , then the spatial line element reduces to 

dl 2 = ^{dx 1 ) 2 + t 2 ^(dx 2 ) 2 + t 2 K(dx 3 ) 2 . (3.4.10) 
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Here p\, p%, p$ are the so-called Kasner indices, satisfying the relations 



pi + vi + P3 = pi + pi + n = 1 ' 



(3.4.11) 



therefore only one of these numbers is independent. Except for the cases (0,0,1) and 
(— V3, 2 /3, 2 /3), such indices are never equal, one of them being negative and two positive; 
in the peculiar case pi = pi = 0, p3 = 1, the metric is reducible to a Galilean form by the 
transformation 

t sinh x 3 = £ , f cosh r 3 = t , 



i.e. with a fictitious singularity in a flat space-time. 
Once that Kasner indices have been ordered according to 

Pi < Pi < P3 , 
their corresponding variation ranges are 

"^<Pl<0, 0<p 2 <^, |<P3<1. 

In parametric form we have the representation 

—u , s 1 + u 



pi (") 



1 + u + u 2 r K ' 1 + u + u 2 
as the parameter u varies in the range (see Fig l3.1|) 

1 < u < +00. 



P3(") 



u(l + u) 
1 + u + u 2 



(3.4.12) 



(3.4.13) 



(3.4.14) 



(3.4.15) 



(3.4.16) 
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Figure 3.1: Evolution of Kasner indices in terms of the parameter l/u. The domain of u 
is [1 , oo); for lower values of u the inversion property (|3.4.17|) holds. 



The values u < 1 lead to the same range following the inversion property 

Pi Qj = Pl (") > P2 Q) = P3 (U) , p 3 = P2 (U) . (3.4.17) 
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The line element from (|3.4.9|) describes an anisotropic space where all the volumes lin- 
early grow with time, while linear distances grow along two directions and decrease 
along the third, differently from the Friedmann solution where all distances contract to- 
wards the singularity with the same behavior. This metric has only one non-eliminable 
singularity in t = with the only exception of the case p\ = p2 = 0, ps = 1 mentioned 
above, corresponding to the standard Euclidean space. 

For a discussion of the Bianchi I model in presence of several matter sources, see the 
following literature ||393|263l!5lll!!![l^ 
IZl^l3^l3Ml3^l4^ITM 
fl36ll453llTT9ll . 



3.5 The role of matter 

Here we discuss the time evolution of a uniform distribution of matter in the Bianchi type 
I space near the singularity; it will result that it behaves as a test fluid and thus it does not 
affect the properties of the solution. 

Let us take a uniform distribution of matter and assume that we can neglect its influence 
on the gravitational field. The hydrodynamics equations describe the evolution |319| as 

Here u l is the four-velocity and a is the entropy density; in the neighbourhood of the 
singularity it is necessary to use the ultra-relativistic equation of state p = e/3, and then 
we get c ~ e 3//4 . 

As soon as all the quantities are functions of time, we have 

| (abcu e 3/i ) = 0, 4e^ + ujl = . (3.5.2) 

From H3.5.2D , we obtain the two integrals of motion 

abcuoe 3/i = const. , w«e 1/4 = const. . (3.5.3) 

From (|3.5.3|l we see that all the covariant components u u are of the same order. Among 
the controvariant ones, the greatest as t — > is u 3 = u$. Retaining only the dominant 
contribution in the identity UjU 1 = 1, we have Mq ~ and, from (|3.5.3|) , 

e ~ -^tj, u a ~ Vab, (3.5.4) 
or, equivalently, for the Kasner solution 

e „ r 2( Pl+P2 ) = r 2(l- P3 ) f Ua „ t d-p 3 )/2 (3 _ 5 _ 5) 

As expected, e diverges as t — > for all the values of p^, except P3 = 1 (this is due to the 
non-physical character of the singularity in this case). 

The validity of the test fluid approximation is verified from a direct evaluation of the 
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components of the energy-momentum tensor Tf, whose dominant terms are 

T ° ~ eu 2 ~ H 1+ ^ , T\ ~ e ~ r 2 ^^ , (3.5.6a) 

7f ~ ew 2 w 2 ~ r^+^-P 8 ) , T 3 3 - eu 3 u 3 - r^ 1+ ^ . (3.5.6b) 

As t — > 0, all the components grow slower than i~ 2 , which is the behavior of the domi- 
nant terms in the Kasner analysis. Thus the fluid contribution can be disregarded in the 
Einstein equations. 

This test character of a perfect fluid on a Kasner background remains valid even in the 
following Mixmaster scenario, both in the homogeneous as well as in the inhomogeneous 
case [54, 57J. The reason for the validity of such extension relies on the piece-wise Kasner 
behavior of the oscillatory regime. For a discussion on the effects of ultra-relativistic 
matter, of scalar and electromagnetic field on the quasi-isotropic solution, see l|375t 13781 
13771 . 



3.6 The Dynamics of the Bianchi Models 

The Kasner solution properly approximates those cases when the Ricci tensor appearing 
in the Einstein equations P a p is of higher order in 1 / 1 with respect to all other terms in- 
volved. However, since one of the Kasner exponents is negative, terms of order higher 
than t~ 2 appear in the tensor P a p. In such a case the discussion of solutions has to be ex- 
tended to the general anisotropic case, in the search of a general behaviour of the Universe 
towards the initial singularity. In fact, the outlined Kasner regime 

relies on a restriction over the phase space of the solution (not discussed here in the de- 
tails, see |54], §3) which causes an instability with perturbations violating this condition. 
A general solution, by definition, 

• is completely stable, i.e. the effect of any perturbation is equivalent to a change of 
the initial conditions at some moment of time and 

• must satisfy arbitrary initial conditions, i.e. the perturbation cannot change the form 
of the solution. 

Nevertheless, the cited restriction over the Kasner solution makes it unstable with re- 
spect to perturbations destroying it: such perturbation promoting the transition to a new 
state cannot be considered small and lies outside the region of the infinitesimal ones. 

Let us introduce three spatial vectors e a = l(x 7 ) , m(x 7 ) , n(x 7 ) and take the matrix h x p 
diagonal in the form 

Kp = a 2 {t)Up + Jp-(t)m a Tnp + c 2 (i)n a «p ■ (3.6.1) 
Consequently, the Einstein equations in a synchronous reference system and for a generic 
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homogeneous cosmological model in empty space are given by the system 

(abc)' 1 



-R m 



-R n 



abc 
(abc)' 

abc 
(abc)' 



2a 2 b 2 c 2 
1 

2a 2 b 2 c 2 
1 



¥ 2 b A 



( ¥ b 2 
(Aa 2 



^ 2 
vc 2 ) 



vc 1 ) 



abc 2a 2 b 2 c 2 



v 2 c 4 - (Aa 2 -fib 2 ) : 









and 



(3.6.2a) 
(3.6.2b) 
(3.6.2c) 

(3.6.3) 



where the other off-diagonal components of the four-dimensional Ricci tensor identically 
vanish as a consequence of the choice of the diagonal form as in (|3.6.1|) . Eventually, the Oa 
components of the Einstein equations can be non-zero if some kind of matter is present, 
leading to an effect of rotation on the Kasner axes |54| . The constants A, ]i,v correspond 
to the structure constants C\\, C22, C33 respectively, introduced earlier in Section l3.1.3l In 
particular, we will study in details the cases of (A, }i, v) for Bianchi type II (1, 0, 0), VII (1, 
1, 0), VIII (1, 1, -1) and IX (1, 1, 1). 

All these equations are exact and contain functions of time only, without any restriction 
regarding the vicinity to the singular point t = 0. Through the notation 



« = In a , /3 = lnb, 7 
and the new temporal variable T defined by 

dt = abc dr , 

Eqs. (I3A21 and d3A3l simplify to 

2ol tt = 
2/3 TT = 

2 7tt = (Aa 2 -]ib 2 ) 



lnc 



( ¥ b 2 
(Aa 2 



vc 



vc 



,2\ 2 



aV 

ll 2 ^ 

v 2 c\ 



(3.6.4) 
(3.6.5) 

(3.6.6a) 
(3.6.6b) 
(3.6.6c) 



- (a + $ + 7 ) TT = a Tj 6 T + a T j T + fi T j T , (3.6.7) 

where subscript t denotes the derivative with respect to t. Manipulating the system 
l|3.6.6|) and using l|3.6.7|) , one obtains the first integral 

a T /6 T + oc T y T + j6 T 7 T = 

= - (A 2 a 4 + }i 2 b 4 + i/ 2 c 4 - 2A^a 2 b 2 - 2Ava 2 c 2 - 2]ivb 2 c 2 ^ (3.6.8) 

involving first derivatives only. The Kasner regime (I3.4.10I) discussed before is the so- 
lution of equations (|3.6.6|) corresponding to neglecting all terms on the right-hand side. 
However, such a situation cannot persist indefinitely as t — > since there are always some 
terms on the right-hand side of (|3.6.6|) which are increasing. 
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3.7 Application to Bianchi Types II and VII 

In this paragraph, the dynamics of the types II and VII is discussed in some detail: these 
spaces, in fact, present some features in common with the Mixmaster model, as it will be 
clear in the Sections I4.1H4.2I 



3.7.1 Bianchi type II 

Introducing the structure constants for the type II model, the system (|3.6.2|) reduces to 



(abc)- 
abc 


a 2 
2b 2 c 2 ' 


(3.7.1a) 


(abc)- 
abc 


a 2 
2b 2 c 2 ' 


(3.7.1b) 


(abc)- 
abc 


a 2 
2b 2 c 2 ' 


(3.7.1c) 


a b 
a b 


b - = o. 

c 


(3.7.2) 



In (|3.7.1|) the right-hand sides play the role of a perturbation to the Kasner regime; if at 
a certain instant of time f they could be neglected, then a Kasner dynamics would take 
place. This kind of evolution can be stable or not depending on the initial conditions; 
as shown earlier in Section 13.41 the Kasner dynamics has a time evolution which differs 
along the three directions, growing along two of them and decreasing along the other. For 
example, for a perturbation growing as a 4 ~ f 4 P", toward the singularity, two scenarios 
are possible: if the perturbation is associated with one of the two positive indices, it will 
continue decreasing till the singularity and the Kasner epoch is stable; on the other hand, 
if p a < 0, the perturbation grows and cannot be indefinitely neglected. In this case, 
the analysis of the full dynamical system is required, and this can be achieved with the 
logarithmic variables (|3.6.4|) - (|3.6.5|) and the system (|3.6.6|) becomes 

*tt = ~\e^ (3.7.3a) 

/3tt = 7tt = \e ia . (3.7.3b) 

The equation (|3.7.3a|) can be thought as the motion of a one-dimensional point-particle 
moving in an exponential potential: if the initial "velocity" da/ ' dx is equal to p„, then the 
potential will slow down, stop and accelerate again the point up to a new "velocity" — p a - 
From there on, the potential will remain negligible forever. 
Furthermore, the second set of equations (|3.7.3b|) implies that the conditions 

«TT + j6 T T = «TT + 7TT = (3.7.4) 
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hold. Considering the explicit solutions of (|3.7.4|) 

1 

oc(r) = - In (cisech (rc\ + c 2 )) (3.7.5a) 
1 

/3(t) = c 3 + tc 4 - - In (cisech (tc x + c 2 )) (3.7.5b) 
1 

7(t) = c 5 + tc 6 - - In (cisech (tci + c 2 )) , (3.7.5c) 

where ci ...,06 are integration constants and we see how this dynamical scheme de- 
scribes two connected Kasner epochs, where the perturbation has the role of changing 
the values of the Kasner indices. Let us assume that the Universe is initially described 
by a Kasner epoch for r — > oo, with indices orderd as p\ < p m < p„; the perturbation 
starts growing and the point bounces against the potential and a new Kasner epoch be- 
gins, where the old and the new indices (the primed ones) are related among them by the 
so-called BKL map (54) 

i _ U>\\ _/ .. 2 \Pl\ - Pm / Pn ~2\pi\ 



Pl -l-2\ Pl \' Vm ~ \-l\px\ ' Vn ~ \-2\pi\ • ( } 

In this new era, the negative power is no longer related to the Z-direction so that the previ- 
ously increasing perturbation is damped and eventually vanishes toward the singularity. 
We will see how l|3.7.6|) will be valid in general. 

For a detailed analysis on the main results concerning the Bianchi II model see the 
following literature 115^ l48l 13^ [3^ [TT^ 14^ [3^ 13^ I2MI [T38l . 



3.7.2 Bianchi type VII 

The analysis of Bianchi type VII can be performed analogously, leading to the Einstein 
equations 

(obey -a 4 + b 4 

abc 2a 2 b 2 c 2 
(abc)- a 4 — b 4 

abc 2a 2 b 2 c 2 ' 
(abc)- __ [a 1 - b 2 ) 2 

abc 2b 2 c 2 

and the constraint H3.7.2D holding unchanged. Comparison of H3.7.7aD with (13.7.11 ) allows 
a similar qualitative analysis: if the right-hand sides of (|3.7.7a|) are negligible at a certain 
instant of time, than the solution is Kasner-like and can be stable or unstable depending 
on initial conditions. If the negative index is associated with the n direction, than the per- 
turbative terms a 4 and b 4 , evolving as t 4 ? 1 and t 4 ? m , decrease up to the singularity and the 
Kasner solution turns out to be stable; in all other cases, one and only one of the pertur- 
bation terms starts growing, blasting the initial Kasner evolution and ending as before in 
a new Kasner epoch. 

The main difference between the types II and VII is that many other transitions can occur 
after the first one and this can happen, for example, if the new negative Kasner index is 
associated with the m direction, i.e. with b. In this case, the b 4 term would start growing 



(3.7.7a) 
(3.7.7b) 
(3.7.7c) 
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and a new transition would occur with the same law (|3.7.6[) . The problem of understand- 
ing if, when and how this mechanism can break up is unraveled considering the BKL 
map written in terms of the parameter u, i.e. 

Pi = Pi(") } [ Pi = P2(u-l) 

Vm = piiu) \ => I p' m = pi(u- 1) (3.7.8) 

Pn = P3(U) J { P'„ = P3(«-l) 

In this representation and from the properties (|3.4.15|) and (|3.4.17|) we see how the exact 
number of exchanges among the I- and m-directions equates the integer part K of the 
initial value u$ describing the dynamics. In fact, for the first Kq times, the negative index 
is exchanged among / and m, then it passes to the n direction, a new and final (toward 
the singularity) Kasner epoch begins, and no more oscillations take place. The collection 
of the total of K epochs is called a Kasner era: in this sense we can say that in the general 
case the type VII dynamics is composed by one era and a final epoch. 

For additional informations about the Bianchi VII model (as well as some interesting 
features regarding the Bianchi VI Universes) see the articles Il32^l420ll440l l68ll504l[38l[282l 
12831 13591 198HTTH. 
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4 Chaotic Dynamics of the Bianchi Types 
VIII and IX Models 



In this Section we provide the detailed construction of the oscillatory-like regime of the 
Mixmaster model while approaching the initial singularity both in the field equations 
formalism and in the ADM one. A relevant part is centered around the chaotic properties 
of its dynamics. The cosmological implementation is also discussed, and in the last three 
Sections the effects on chaos of matter fields and of the number of space dimensions is 
reviewed. 



4.1 Construction of the solution 

At this point we are going to address the solution of the system of equations (|3.6.2|) for the 
cases of Bianchi types VIII and IX cosmological models, following the standard approach 
of Belinsky, Khalatnikov and Lifshitz (BKL) [54 . 57^. Although the detailed discussion is 
devoted to the Bianchi IX model, it can be straightforwardly extended to the type VIII. 
Explicitely, the Einstein equations (|3.6.2|) reduce to 

a 4 (4.1.1a) 
b 4 (4.1.1b) 
c 4 , (4.1.1c) 

togheter with the constraint l|3.6.7|l unchanged, leading to the consequent constant of mo- 
tion 

cc T p T + a r7T + £ t7t = - (V + & 4 + c 4 - 2a 2 b 2 - 2a 2 c 2 - 2b 2 c 2 ) . (4.1.2) 

Let us therefore consider again the case in which, for instance, the negative power of the 
Pi's exponents corresponds to the function a{t) (that is to say p\ = p\)\ the perturbation 
of the Kasner regime results from the terms A 2 a 4 (remember that A = 1 for both models) 
while the other terms decrease with decreasing t, in fact 

{a ~ — | pi | In t 
1 / forf^O (4.1.3) 

a ~ -. r 
t\vi\ 



For a thorough discussion of works which concern the dynamics of such models see the literature 15171 

l447ll4iniM7llll^lT8aiT^I^ 

424 . 215 62 16 . 450 425, 76 77 371 . 14 . 189 . 415] 



2a TT 


= {b 2 


-c 2 ) 2 




= (a 2 


-c 2 ) 2 


27tt 


= (a 2 


-b 2 ) 2 
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and along the other directions 

P2>^P2 = \P2l ^|p 2 |ln f \ forf ^ Q 
p 3 > -> p 3 = |p 3 |, 7~|p 3 |ln£\ 

Preserving only the increasing terms on the right-hand side of equations (|3.6.6|) , we obtain 
a system identical to l|3.7.3a|) , whose solution describes the evolution of the metric from 
its initial state (|3.4.10[) . Let us fix, for instance, 

Pl = Pl, Pm = P2, Pn = P3, (4.1.5) 

so that 

a ~ , b~t? 2 , c ~ 0* , (4.1.6) 

and then 

abc = At 
1 

T = — In t + const. (4.1.7) 
A 

where A is a constant, so that the initial conditions for (|3.7.3aD can be formulated as 

a T = Api , |S T = Apr , j T = Ap 3 , (4.1.8) 
for t — > oo . 

The system (|3.7.3a|) with (|4.1.8|) is integrated to 



a 1 = — - , r (4.1.9a) 
cosh (2 | p x | At) 

b 2 = b 2 exp [2A (p 2 - \ p\ \)x\ cosh (2 | p x \ At) (4.1.9b) 

c 2 = c 2 exp [2A (p 3 - | p! |) t] cosh (2 | pi | At) (4.1.9c) 



where bo and Co are integration constants. 



4.2 The BKL oscillatory approach 

Let us consider the solutions (|4.1.9|) in the limit T — > oo: towards the singularity they 
simplify to 



that is to say, in terms of t, 



a ~ exp [— Apir] 
b ~ exp [A(p 2 + 2pi)T] 
c ~ exp [A(p 3 + 2pi) t] 
t ~ exp [A(l + 2pi) t] 



(4.2.1a) 
(4.2.1b) 
(4.2.1c) 
(4.2.1d) 



a ~ t p > , b ~ , c ~ f p » , flbc = A'f, (4.2.2) 
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where the primed exponents are related to the un-primed ones by 

J _ I Pi I v ' _ 2 I Pi I ~n f42 3a) 

P '-l-2| Pl |' Vm ~ \-2\vi\' ( } 

v'n = n ~l\l X ) , A' = (1-2 I Pl |) A, (4.2.3b) 



which, we note, is very similar to what found in the type II case with relations ( I3.7.6D plus 
the additional expression involving A. Summarizing these results, we see the effect of the 
perturbation over the Kasner regime: a Kasner epoch is replaced by another one so that 
the negative power of t is transferred from the 1 to the m direction, i.e. if in the original so- 
lution p\ is negative, in the new solution p' m < 0. The previously increasing perturbation 
A 2 fl 4 in (|3.6.2|) is damped and eventually vanishes. The other terms involving y? instead 
of A 2 will grow, therefore permitting the replacement of a Kasner epoch by another. Such 
rules of rotation in the perturbing property can be summarized with the rules (|3.7.8|) of 
the BKL map, with the greater of the two positive powers remaining positive|55J. 
The following interchanges are characterized by a sequence of bounces, with a change of 
the negative power between the directions 1 and m continuing as long as the integral part 
of the initial value of u is not exhausted, i.e. until u becomes less than one. In terms of the 
parameter u, the map (|3.7.6|) takes the form 

1 

u = u — 1 for u > 2 , u' = for u < 2 . (4.2.4) 

u — 1 _ v ' 

At that point, according to (|3.4.17|) , the value u < 1 is turned into u > 1, thus either 
the exponent pi or p m is negative and p n becomes the smaller one of the two positive 
numbers, say p n = pi- The next sequence of changes will bounce the negative power 
between the directions n and 1 or n and m. 

The phenomenon of increasing and decreasing of the various terms with transition 
from a Kasner era to another is repeated infinitely many times up to the singularity. Let 
us analyze the implications of the BKL map (|3.7.8[) and of the property l|3.4.17[) . 

If we write u = k + x as the initial value of the parameter u, with k and x being its 
integral and fractional part, respectively, the continuous exchange of shrinking and en- 
larging directions l|3.7.8|) proceeds until u < 1, i.e. it lasts for k° epochs, thus leading a 
Kasner era to an end. The new value of u is u' = 1 fx > 1 (|3.4.17|) and the subsequent set of 
exchanges will be l-n or m-n: for arbitrary initial values of u, the process will last forever 
and an infinite sequence of Kasner eras takes place. 

For an arbitrary, irrational initial value of u the changes l|3.7.8|) repeat indefinitely. In 
the case of an exact solution, the exponents p\, p m and p n loose their literal meaning, thus 
in general, it has no sense to consider any well defined, for example rational, value of u. 

The evolution of the model towards the singularity consists of successive periods, the 
eras, in which distances along two axes oscillate and along the third axis decrease mono- 
tonically while the volume decreases following a law approximately ~ t. The order in 
which the pairs of axes are interchanged and the order in which eras of different lengths 
(number of Kasner epochs contained in it) follow each other acquire a stochastic charac- 
ter. Successive eras 'condense' towards the singularity. Such general qualitative proper- 
ties are not changed in the case of space filled in with matter, however the meaning of the 
solution would change: the model so far discussed would be considered as the principal 
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terms of the limiting form of the metric as t — > 0. 



4.3 Stochastic properties and the Gaussian distribution 

a decreasing sequence of values of the parameter u corresponds to every s-th era there. 

This sequence, from the stc 
can introduce the notation 



This sequence, from the starting era has the form Umax, u^lx — 1/ ««« — 2, ... , . We 



u ( s ) = + x& (4.3.1) 

then 



"it = * (S) < 1 ' "i- = k(3) + ^ (S) ' ( 4 - 3 - 2 ) 



where u„ iax is the greatest value of u for an assigned era and £;( s ) = u^lx (square brackets 

denote the greatest integer less or equal to u„l x ). The number denotes the era length, 
i.e. the number of Kasner epochs contained in it. For the next era we obtain 



(8+1) _ J_ , (s + l) 

u max — c-i / 
X\ s > 



(4.3.3) 



•(») 

For large u, the Kasner exponents approach the values (0, 0, 1) with the limiting form 

11 1 

Vi ~ — / P2 ~ - , P3 ~ 1 - -j / ( 4 -3-4) 
u u u 

and the transition to the next era is governed by the fact that not all terms in the Einstein 
equations are negligible and some terms are comparable: in such a case, the transition is 
accompanied by a long regime of small oscillations [54 J lasting until the next era, whose 
details will be discussed in Section l4~4l after which a new series of Kasner epochs begins. 
The probability A of all possible values of x^ which lead to a dynamical evolution to- 
wards this specific case is strongly converging to a number A <C l[326j. If the initial value 
of is outside this special interval for A, the special case cannot occur; if x^ lies in 
this interval, a peculiar evolution in small oscillations take place, but after this period the 
model begins to regularly evolve with a new initial value x^°\ which can only acciden- 
tally fall in this peculiar interval (with probability A). The repetition of this situation can 
lead to these cases only with probabilities A, A 2 , . . ., which asymptotically approach zero. 

If the sequence begins with k^ + x^°\ the lengths k^ 2 \ . . . are the numbers appear- 
ing in the expansion for x^ in terms of the continuous fraction 

= , (4.3.5) 

)tW + ^ 

jfc(2) + ' 



m + ... 

which is finite if related to a rational number, but in general it is an infinite one [287J. 

For the infinite sequence of positive numbers u ordered as (|4.3.3|) and, admitting the 
expansion (|4.3.5|) , it is possible to note that 

i) a rational number would have a finite expansion; 
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ii) periodic expansion represents quadratic irrational numbers (i.e. numbers which are 
roots of quadratic equations with integral coefficients) 

iii) irrational numbers have infinite expansion. 

All terms k^',k^ 2 ' ,k^\ ... in the first two cases having the exceptional property to be 
bounded in magnitude are related to a set of numbers x^ < 1 of zero measure in the 
interval (0,1). 

An alternative to the numerical approach in terms of continuous fractions is the statistic 
distribution of the eras' sequence for the numbers x' ' over the interval (0, 1), governed 
by some probability law. For the series x( s > with increasing s these distributions tend to 
a stationary one zv(x), independent of s, in which the initial conditions are completely 
forgotten 

= (TTFjTn2 ' <4 ' 3 - 6) 

In fact, instead of a well defined initial value as in (14.3.11) with s = 0, let us consider a 
probability distribution for x^ over the interval (0, 1), Wo(x) for x^ = x. Then also the 
numbers x^ are distributed with some probability. Let w s (x)dx be the probability that 
the last term in the s-th series x^ = x lies in the interval dx. The last term of the previous 
series must lie in the interval between l/(k + 1) and 1/k, in order for the length of the 
s-th series to be k. 

The probability for the series to have length k is given by 



. i 

W,(k) = T w s - 1 (x)dx. (4.3.7) 



For each pair of subsequent series, we get the recurrent formula relating the distribution 

w s+ i(x) to w s (x) 



w s+1 (x)dx = f> 5 (j^Ti) 
or, simplifying the differential interval, 



1 



(k + x) 



(4.3.8) 



If for increasing n the zv s+n distribution (|4.3.9|) tends to a stationary one, independent of 
s,w(x) has to satisfy 

CO , / 1 \ 

w ( x ) = V rw . (4.3.10) 

A normalized solution to (I4.3.10D is clearly given by d4.3.6l) l54|; substituting it in (14.3.71 ) 
and evaluating the integral 

WW = /j^ W , I=s L ta iL±J, (4.3.U) 

we get the corresponding stationary distribution of the lengths of the series k. Finally, 
since k and x are not independent, they must admit a stationary joint probability distri- 
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bution 



w (k, x) = — — (4.3.12) 

y ' {k + x) (fc + x + l)ln2 v ; 



which, for u = k + x, rewrites as 

w{u)= . jpr. . , (4.3.13) 
u(u + l)ln2 

i.e. a stationary distribution for the parameter u. 

The existence of the Gauss map was firstly demonstrated in the work of Belinskii, Kha- 
latnikov and Lifshitz [54| , showing how a statistical approach [295 J describes the time 
evolution of the cosmological models near the singularity. These features opened the way 
to further investigations (see, for example |32|), in view also of the peculiar properties of 
the discrete map w leading to final form for a measure of the full degrees of freedom of the 
discrete Mixmaster dynamics expressed as a not separable function |3"T| , and such map 

• has positive metric- and topologic-entropy; 

• has the weak Bernoulli properties (i.e., the map cannot be finitely approximated); 

• is strongly mixing and ergodic [385 J. 

Thus, we see that the approach to the initial singularity is characterized by an infinite 
series of Kasner epochs described by the stochastic properties of the associated map. 
For a discussion of the chaotic behavior inherent the homogeneous early cosmologies, 

see llilyfll2l^lTMl4^ITWl l^ 
11201 147^14571 1! I35mil53iram 



4.4 Small Oscillations 

Let us investigate a particular case of the general (homogeneous) solution constructed 
above. We analyze an era during which two of the three functions a, b, c (for example 
a and b) oscillate so that their absolute values remain close to each other and the third 
function (in such case c) monotonically decreases, so that c can be neglected with respect 
to a and b. As before, we will discuss only the Bianchi IX model, since for the Bianchi VIII 
case the argumentations and results are qualitatively the same I54ll57l . Let us consider 
the equations we can obtain from (|3.6.6|) and (|3.6.8[) 

0, (4.4.1a) 
e^-e itt , (4.4.1b) 

- aT B T + l( e 2a -e 2 ?) 2 . (4.4.1c) 
The solution of (|4.4.1a|) is 

2a 2 

cc + B = -0( T -T )+21n(flo), (4.4.2) 
Co 



0L TT B TT = 

j T (tx T + B T ) = 
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where flo and £q are positive constants. For what follows we conveniently replace the time 
coordinate x with the new one £ defined as 

£ = Soexp(|t( T _ To )^ (4.43) 
in terms of which the equations (|4.4.1b|) and (|4.4.1c|) rewrite as 

X& +?Xi+2 sinh ( 2 ^) = ' < 4A4a) 

n = ~~k + i ( 2 ** + cosh(2;c) " x ) ' (4A4b) 

where we have introduced the notation x = * — jS and = d/L / Since T is defined 
in the interval [To, -co), from equation (|4.4.3|) we note that f G [£o/0). Since a general 
analytic solution of the system 04.4.41) is not available, we shall consider the two limiting 
cases £ ^> 1 and £ <C 1 only. 



Let us start with the £ ^> 1 region. In this approximation the solution of equation 
(|4.4.4a|) reads as 

2/4 

X = -=sin(£-£ ), (4.4.5) 

A being a constant and therefore obtaining 7 ~ A 2 (£ — £o)- As we can see, the name 
"small oscillations" arises from the behavior of the function X- The functions a and b, i.e. 
the expressions for the scale factors, are straightforwardly obtained as 

a,fc = a ^(l±^sin(£-£o)) , (4.4.6a) 

c = c exp[-A 2 (£ -£)] • (4.4.6b) 

The synchronous time coordinate t can be re-obtained back from the relation dt = abc dr 

as 

t = f exp[-A 2 (£ -£)] • (4.4.7) 
Of course, these solutions only apply when the condition Cq <C «o is satisfied. 



Let us discuss the region where £ <C 1. In such a limit, the function # reads as 

X = K\n£ + 6, 6 = const. (4.4.8) 

where K is a constant which, for consistency, is constrained in the interval K G (—1,1). 
We can therefore derive all other related quantities, and in particular 

fl ^(l+K)/ 2/ b „gl-K)/2 t c _ r (l-^)/ 4/ f „£(3+K*)/2_ (449) 

This is again a Kasner solution, with the negative power of t corresponding to c and the 
evolution is the same as the general one. Moreover, we can easily note how, for such a 
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Kasner epoch, the parameter u (introduced in (|3.4.15|l ) becomes 

for K > : u = , for K < : u = \— | . (4.4.10) 

1-K l+K y ' 

Summarizing, the system initially crosses a long time interval during which the functions 
a and b satisfy [a — V) I a < l/£ and perform small oscillations of constant period A£ = 
2n, while the function c decreases with t as c = cof/fo- When £ ~ 0(1)/ equations 
( |4.4.6aD and d4.4.6bl ) cease to be valid, thus after this period the c starts increasing and a 
and b decreasing. At the end, when condition c 2 / (ab) 2 ~ t~ 2 is realized, a new period of 
oscillations is reached (Kasner epochs) and the natural evolution of the system is restored. 

Let us stress how the matching of the constants of the above limit regions is possible. 
In particular we can express the constants (K, 8) in terms of (A, go)- Such procedure can 
be implemented also by replacing sinh(2^) with 2% < 1) in equation (|4.4.4a|) , and 
finally the solution is expressed in terms of Bessel functions which can be asymptotically 
expanded and compared with the explicit mentioned solutions JT5). 



4.5 Hamiltonian Formulation of the Dynamics 

So far we reviewed the formulation of General Relativity in terms of geometrodynamics, 
i.e. studying the equations governing the evolution of the metric tensor. In what follows 
we will consider a different interpretation of the terms involved, leading to a peculiar 
Hamiltonian formulation. This paragraph is devoted to specialize it to the Mixmaster 
dynamics, i.e. to the Bianchi VIII and IX models 155511557115701 . 

The space-time line element adapted to the 3+1 foliation can be written as in (|2.4.3|) I 370I 
where, for a spatially homogeneous model, the three-metric tensor h a a is parametrized by 
three functions q a = q a (t), in particular as 

h«p = J'SdKWfyx) , (4.5.1) 

although such symmetry is valid in the case of absence of matter only. In fact, the Einstein 
equations in vacuum Rq k = allow the choice of a diagonal metric h x p, thus leaving three 
degrees of freedom only. The presence of matter induces in general a non-diagonal metric 
(which will be analyzed later) while here we focus on the vacuum one. 

Given the parametrization (|4.5.1|) , considering that all the canonical variables, including 
the lapse function and the shift vector, are independent of the spatial coordinates and that 
this model identically satisfies the super-momentum constraint, we get the action as 

S = J dt (p a d t q a - NH) , (4.5.2) 

Va = Pa(t) being the momenta canonically conjugate to q a . The super-Hamiltonian H 
reads explicitly as 




(4.5.3) 
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where r\ = exp q a and the Ricci three-scalar ( 3 '.R plays the role of the potential in the 
dynamics. 

Let us introduce the "anisotropy parameters", which will be convenient in the following, 

as 

The functions defined as in (|4.5.4|) allow a clearer interpetation of the last term in the 
right-hand side of equation (|4.5.3|) as a potential form U for the dynamics. In fact, it can 
be re-written in the more expressive form 

U = tj®R = Yj X2 aV 1Q " ~ E hKn Qb+Qc > (4-5.5) 

a b^c 

where the constants X a specify the model under consideration, being X a = (1,1,-1) or 
A n = (1, 1, 1) for Bianchi VIII and IX, respectively. 

The main advantage of writing the potential as in (|4.5.5I) , arises when investigating its 
proprieties in the asymptotical behavior toward the cosmological singularity (rj — > 0). In 
fact, the second term in (|4.5.5[) becomes negligible, while the value of the first one results 
to be strongly sensitive to the sign of Q„s and the potential can be modeled by an infinite 
well as 

U = Y2®(Q«) (4-5.6) 



where 



»w= r (4.5.7) 

0, if x > . 



By (|4.5.6|) we see how the dynamics of the Universe resembles that of a particle moving in 
a dynamically-closed domain ITq, defined as the one where all the anisotropy parameters 
Q a are simultaneously positive. 



4.6 The ADM Reduction of the Dynamics 

The ADM reduction of the dynamics relies on the idea of identifying a temporal parame- 
ter as a functional of the geometric canonical variables, before applying any quantization 
procedure or, in other words, of solving the classical constraints before quantizing. For 
the moment, we will simply discuss the reduction of General Relativity to a pure canoni- 
cal form focusing on the elimination of non-dynamical variables. 

Let us count the degrees of freedom of the gravitational field. We have twenty phase- 
space variables, given in the 3+1 formalism by the set (N,n;N a ,n a ;/z a/3 ,n^) ( see Sec- 
tion 123]), subject to eight first-class constraints (n = 0,Yl a = 0, H = 0,H a = 0). Since 
each constraint eliminates two phase-space variables, we remain with four of them, cor- 
responding to two physical degrees of freedom of the gravitational field, i.e. to the two 
independent polarizations of a gravitational wave in the weak field limit. 

After eliminating N and N a , we have 12 x oo 3 variables (h a ^(x),Yl a -^(x)). Then, we 
can remove 4 x oo 3 variables thanks to the secondary constraints (|2.4.14|) . The remaining 
4 x oo 3 non-physical degrees of freedom result from the lapse function and the shift vector. 
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In analogy with the Yang-Mills theory, also in this case it is necessary to impose some sort 
of gauge, fixing the lapse function and the shift vector. 

In order to obtain a true canonical form for the canonical theory, we have to follow the 
key steps listed by Isham in [280 j. With respect to this formulation, we refer to a canonical 
description of the physical degrees of freedom only and with the following procedure we 
are able to discard the non-physical variables. 

i) Perform a canonical transformation 

(h^{x),Il«P{x)) -» (x A {x),P A {x);<p r {x),n r {x)) , (4.6.1) 
A = 1,2,3,4 r = l,2 

where X A { X ) define a particular choice of the space and time coordinates^, P A (x) 
are the corresponding canonically conjugate momenta and the four phase space 
variables {<p r (x) , n r {x)) represent the physical degrees of freedom of the system. 
We emphasize that these "physical" fields are not Dirac observables, in the sense 
defined in Section [2~4l 

Consequently, the symplectic structure is determined by 

{ X A (x),P B (x')} = 5 A 5\x-x'), (4.6.2a) 
{f(x), n s (x')} = S' s S 3 (x - x') , (4.6.2b) 

while all other Poisson brackets are vanishing. 

ii) Express the super-momentum and super-Hamiltonian in terms of the new fields 
and then write the Lagrangian density as 

L'(N, N\x A , Pa, f, 7t r ) = P A X A + ~ NH ' ~ N a H' a . (4.6.3) 

iii) Remove 4 x oo 3 variables arising from the constraints H = and H a = once 
solved with respect to P A the equation 

P A {x)+h A (x,x,<p,n) =0 (4.6.4) 

and re-inserting it in (|4.6.3[) . Removing the remaining 4 x oo 3 non-dynamical vari- 
ables we yield the so-called true Lagrangian density |280| 

&true = 7T r (p r ~ ^AX. A / ( 4 -6-5) 

where the lapse function and the shift vector do not play any role, but specifying the 
form of the functions x A ■ After solving the constraints, there is no longer informa- 
tion about the evolution of x A i n terms of the parametric time t. Thus, in equation 
(|4.6.5|) we have chosen the conditions x (x,t) = xt i x ) an d the true Hamiltonian is 

Htruc = J d 3 x X A {x)h A {x,xt,<p{t),n{t)) . (4.6.6) 

2 These fields can be interpreted as defining an embedding of Z in M via some parametric equations. 
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From (|4.6.6|) one can derive the equations of motion as 

dtp = {<p r ,Htrue}red (4.6.7a) 
d t n s = {7T Sr H true } red , (4.6.7b) 

where the notation {...} re d refers to the Poisson brackets evaluated in the reduced 
phase space with coordinates given by the physical modes p and n s . 

This is an operative method to classically solve the constraints, i.e. pulling out all the 
gauges, and obtaining a canonical description for the physical degrees of freedom only. 

We also note that such procedure violates the geometrical structure of General Relativ- 
ity, since it removes parts of the metric tensor. Nevertheless, it avoids any obstacle at a 
classical level, but there are several problems when implemented in the quantum frame- 
work B280U315[|26| . In particular, they are related to the one of the reduced-phase space 
quantization, for more details see |251|) . 



4.7 Misner variables and the Mixmaster model 



In order to implement the formalism to our purposes, we start introducing a new set of 
variables, as mentioned in expression (|3.2.8|) . In fact, we parametrize the tetradic projec- 
tion of the metric as 



\ / ab 



(4.7.1a) 



where /3 8 & is a three-dimensional matrix with null trace of the kind diag(/3n, ^22, (633) and 
the exponential matrix has to be intended as a power series of matrices, so that 



detO) =e 2 ^ = l 



(4.7.2) 



and 



rj = e 6x . (4.7.3) 

The BKL formalism used in Section l3T4l matches the present one (for the Bianchi I solution) 
once considering the expression 



Pab 



d (In//) 



ab 



din, 



From equations J4.7.1a| > and 04.7.41 ) the relation 



Pab = g 



&ab + 



da 



follows, so that the first Kasner condition (|3.4.11|) rewrites as 

l = EF-tr Pflb = l + Itr(|) 



(4.7.4) 



(4.7.5) 



(4.7.6) 



53 



4 Chaotic Dynamics of the Bianchi Types VIII and IX Models 



which is an identity since the trace tr/S fl b = 0. 
The second Kasner relation <|3.4.11|) rewrites as 



tr (p 2 ) = 1 (4.7.7) 



which, by virtue of (|4.7.5|) becomes 



and then 

f! 

\ da 

which is no longer an identity but a consequence of the Einstein equations in empty space. 



df5ah ^ = 6 (4.7.9) 



4.7.1 Misner Approach to the Mixmaster 

Once seen how to perform the convenient change of variables for the diagonal case, we 
step further to approach the Mixmaster using this formalism. 
The matrix fi^ has only two independent components defined as 

Pn = + (4.7.10a) 
j6 2 2 = (4.7.10b) 
/3 33 = -2/3+, (4.7.10c) 

and then the Kasner relation for Bianchi I reads as 

The variables f>± together with a are the Misner coordinates. The relation 04.7.11D in terms 
of the Kasner exponents now is 

d ? + = hi -3p 3 ) (4.7.12a) 



doc 2 



d^ 



da 2 
or equivalently with the u parameter 



IV3(PI-P2) (4.7.12b) 



da 2 " " 1 + u + u 2 

Such quantities represent the anisotropy velocity $ 



d 4±- = -l + l—h_ (4.7.13a) 

da 21 + u + u 2 v ' 



(4.7.13b) 
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which measures the variation of the anisotropy amount with respect to the expansion as 
parametrized by the cl parameter. The volume of the Universe behaves as e 3a , tends to 
zero towards the singularity and is directly related to the temporal parameter. 
Eventually the presence of matter as well the effects of the spatial curvature can lead the 
norm | \ | to a deviation from the Kasnerian unity. 

In order to develop a general metric for a homogeneous space-time we rewrite the line 
element in the general form 

ds 2 = N 2 (t)dt 2 - e 2u (e 2 f i ) ab co a ® co b . (4.7.15) 

The cosmological problem reduces to the equations involving the functions oc, N, B a i in 
terms of the independent temporal parameter t, independently of the spatial coordinates. 
Explicitly the dual 1-forms associated to the Bianchi types VIII and IX are, respectively 

co 1 = — sinh ip sinh 6 dtp + cosh tpdO 
(VIII) co 2 = -coshxpsinhedcp + sinhtpdO (4.7.16a) 
a; 3 = cosh Qdtp + dip 



co = sin ip sin 6d(p + cos xpd8 
(IX) co 2 = -cosxpsm9d(p + sintpdO (4.7.16b) 
a; 3 = cos 8dcp + dtp . 

For example, the — Einstein equation with N = 1 reads as 

3 (a 2 - /3V 2 - fSJ) + \ ( ^Rb) = *A , (4.7.17) 

where is the curvature scalar for the three-dimensional spatial surface correspond- 

ing to t = const, and the index B refers to the symmetry properties for the Bianchi cosmo- 
logical models. Such term contains the peculiar difference between the nine types of the 
Bianchi classification, to be evaluated through expressions (|3.3.2|) in terms of the structure 
constants. 

For the models referring to types VIII and IX such curvature scalar reads as 

^R VIU = V 2a (4^+-2tr e - 2 ^-tre 4 ^ (4.7.18a) 

^R IX = i e - 2a tr (le- 2 ? - e 4 ^ , (4.7.18b) 

respectively, and the trace operation has to be intended over the exponential of diagonal 
matrices, without ambiguity. 

Equation (|4.7.17|) with (|4.7.18b|) or (|4.7.18a[) can be interpreted as a contribution of the 
anisotropy energy, connected to the term T° , to the volume expansion a 2 , so that it ap- 
pears as a potential term together with the kinetic ones fi 2 ± . Close to the singularity, this 
term becomes negligible for small values of the anisotropy variables B±. 
Finally, equation (|4.7.17|) has to be regarded as a constraint over the field equations. 
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4.7.2 Lagrangian Formulation in the Misner Variables 

Despite the variables ft introduced in Section 14.51 are linked to the Misner ones by the 
linear transformations 



nevertheless we conveniently restate ab initio the Lagrangian approach. From the general 
expression of the action <|2.1.11|) in vacuum, we get the variational principle 

8 [ t2 Ldt = 0, (4.7.20) 

in which t\ and t^ {ti > t\) are two values of the temporal coordinate and integration is 
performed with respect to the spatial coordinates. 

In particular, integration for the Bianchi type VIII is considered over a spatial volume 
(47r) 2 , in order to have the same integration constant used for the type IX (and to keep a 
uniform formalism), having 



ft = 2a + 2/3+ + 2\/3/S 
ft = 201 + 2,8+ -2a/3^ 
ft = 2a - 4yS , 



(4.7.19a) 

(4.7.19b) 
(4.7.19c) 




(4.7.21) 



This way, the Lagrangian L is written a^l 




where ()' = — and is a function linear in having a potential role. The varia- 

tional principle rewrites explicitely as 




(4.7.23) 



in which H is given by 




(4.7.24) 



and the potential V as 




(4.7.25) 



where is specified for the two Bianchi models under study as 




) 



(4.7.26b) 



(4.7.26a) 



Here we choose the Einstein constant k — 16n, i.e. G = 1. 
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From Lagrangian (|4.7.22|) it is standard to derive the conjugate momenta 



4.7.3 Reduced ADM Hamiltonian 

In order to obtain Einstein equations, the variational principle requires SS to be null for 
arbitrary and independent variations of p±, p U/ f>±, ot, N. Variation with respect to N leads 
to the super-Hamiltonian constraint H = 0. As we have seen in Section l4~6l the procedure 
prescribes the choice of one of the field variables, or one of the momenta, as the temporal 
coordinate and subsequently one can solve the constraint (|4.7.24|) with respect to the cor- 
responding conjugate quantity. 

It is customary, as in this general approach, to set t = ot and solve H = obtaining 



"Kadm = ~p K = ^/pi + pi+V. (4.7.28) 

Within this equation a relation between the temporal gauge described by the function N 
and the dynamical quantity "Kadm is defined. 

Through (I4.7.28D we explicit p K in the action integral, so that the reduced variational 
principle in a canonical form reads as 

SS red = (4.7.29) 

f>red being written as 

Sred = J (p+dp+ + Ml 6 - " X ADM dK) (4.7.30) 
together with the equation defining the temporal gauge for k = 1 

12ne 3a 

N ADM = ^ ■ (4.7.31) 



4.7.4 Mixmaster Dynamics 

In the present Section we will write in general the approach to the Mixmaster dynamics 
and later on it will be applied to prove specific properties, such as chaoticity in a covariant 
approach, with respect to the temporal gauge, and subsequent statistical effects. 

The Hamiltonian introduced so far differs from the typical expression of classical me- 
chanics for the non positive definiteness of the kinetic term, i.e. the sign in front of p\, 
and for the peculiar form of the potential as a function of a (say time) and f>±, reduced to 
the study of a function of the kind V (f} + , ft- ) . 
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Figure 4.1: Equipotential lines of the Bianchi type VIII model in the f} + , jS_ plane. 



A Hamiltonian approach permits to derive the equations of motion as 



P± = 



dH 

dp* ' 
dH . 

dp± 



dH 

do. 
dH 



(4.7.32a) 
(4.7.32b) 



This set considered with the explicit form of the potential (see Figure |4~T1 and |4~2|) , can be 
interpreted as the motion of a "point-particle" in a potential. The term V is proportional 
to the curvature scalar and describes the anisotropy of the Universe, i.e. in the regions 
of the configuration space where it can be negligible the dynamics resembles the pure 
Kasner behaviour, corresponding to |/S'| = 1. 




Figure 4.2: Equipotential lines of the Bianchi type IX model in the j6_ plane. 
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In general, it is necessary a detailed study of the potential form which behaves as a 
potential well in the plane /5+, /3_ . 

Asymptotically close to the origin, i.e. /3± = 0, equipotential lines for the Bianchi type IX 
are circles 

U IX (fi + ,P-) ~ -3+ (/3 + 2 + J 6_ 2 ) +o( i 6 3 ), (4.7.33) 
while for Bianchi type VIII are ellipses 

U ym QS +/J 6_) ~ (40j3+ 2 + 24j6_ 2 ) - 8 (,6+ + 0_) + 5 + o (j6 3 ) . (4.7.34) 
The expressions for large values of /3 are the same for both types 

U (fi) ~ <T 8 P+ , j6+ — ► -oo (4.7.35) 



or 

U (/3) ~ 48£_ 2 e 4/3 + , j6+ — ► +oo (4.7.36) 

when 

|/3_ |< 1. (4.7.37) 

In the figures |4~T1 and l4~2l are represented some of the equipotential lines U (jS) = const., 
for which the potential value has an increment of a factor e 8 ~ 3 x 10 3 for AjS ~ 1. 
The Universe evolution is described as the motion of a point-like particle under the influ- 
ence of such potentials and it corresponds to bounces on the potential walls when evolv- 
ing towards the singularity. The behaviour of the anisotropy variables /3± in this regime 
consists of a Kasner epoch followed by a bounce and then a new epoch with different Kas- 
ner parameters, in correspondence with the description given in Section l4~2l according to 
the BKL approach. 



Let us describe in more detail the bounces performed by the billiard ball representing 
the Universe. From the asymptotic form l|4.7.35|) for the Bianchi IX potential, we get the 
equipotential line f> wa \\ cutting the region where the potential terms are significant. In 
particular we have 

^«// = |-^ln(3H 2 ), (4.7.38) 

where the super-Hamiltonian H is defined as in (I4.7.24D . Since we have H = const, inside 
the potential well, from (|4.7.38[) one gets \P' wa u\ =1/2, i.e. the /3-point moves twice as 
fast as the receding potential wall, i.e. the particle will collide against the wall and will be 
reflected from one straight-line motion (Bianchi I) to another one. 



A relation of reflection-type lays for the bounces |368|. In fact, considering the super- 
Hamiltonian (|4.7.24|l in the (|4.7.35[) -case, if 0, and Qt are the angles of incidence and of 
reflection of the particle off the potential wall, respectively, the relation 

1 

sin0/ - sin0,- = - sin(0; + 6 f ) , (4.7.39) 

holds and in terms of the parameter u introduced in section l3.7.21 it is nothing but Uf = 

Uj — 1. 
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4.8 Misner-Chitre-like variables 



A valuable framework of analysis of the Mixmaster evolution, able to join together the 
two points of view of the map approach and of the continuous dynamics evolution, relies 
on a Hamiltonian treatment of the equations in terms of Misner-Chitre variables |131|. 
This formulation allows to individualize the existence of an asymptotic (energy-like) con- 
stant of motion once performed an ADM reduction. By this result, the stochasticity of 
the Mixmaster can be treated either in terms of the statistical mechanics (by the micro- 
canonical ensemble) |276], either by its characterization as isomorphic to a billiard on a 
two-dimensional Lobachevsky space fT7\ and such scheme can be constructed indepen- 
dently of the choice of a time variable, simply providing very general Misner-Chitre-like 
(MCI) coordinates |273H272|. 



To this purpose, let us here re-define the anisotropy parameters Q a (following 112991 
13021 . (3Z1I3Q3)) as the functions 

i a _l . Fx a 

(4.8.1a) 
(4.8.1b) 
(4.8.1c) 



Qi = 


l 

3 


+ «+- 


h v 7 ^- 
3a 


Qi = 


l 

3 


+ "+- 


- V^jS- 
3a 


Q 3 = 


1 

3 


2/3+ 
3a 





excluding the pathological cases when two or three of them coincide. 
We then introduce the Misner-Chitre variables {t, £, 6} as 

a = -e T cosh£ (4.8.2a) 
,8+ = e T sinh£cos0 (4.8.2b) 
jS_ = e T sinh£sin0 (4.8.2c) 

where < £ < oo, < 9 < 2rc, and T plays the role of a "radial" coordinate coming 
out from the origin of the |6± space 11370 1 . In terms of l|4.8.2|l , it is possible to study the 
first interesting approximation of the potential (|4.7.26|) as independent of T towards the 
singularity, i.e. for a — > — oo. 

To discuss the contrasting results concerning chaoticity and dynamical properties which 
arose from numerics [2031 14431 l64l I257| , it is necessary to introduce a slight modification 
to the set (|4.8.2[) via the MCI coordinates {T(t), 6} through the transformations 

ol = -e Y ^l (4.8.3a) 
fi + = e T ^J^-\cose (4.8.3b) 



j6_ = e r ( T Y^ 2 -lsin0 (4.8.3c) 

where 1 < ^ < oo, and T(t) stands for a generic function of t: in the original work, Chitre 
took simply T(t) = T and set ^ = cosh^". 

This modified set of variables permits to write the anisotropy parameters (|4.8-H> Q a as 
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independent of the variable T as 

Qi = 
Qi = 
Q 3 = 



3 3g 
3 



i +2 



3? 
3£ 



COS0 + sinK 



cos ft -V3 sine* 



COS0 . 



(4.8.4a) 
(4.8.4b) 
(4.8.4c) 



A/Z dynamical quantities, if expressed in terms of (14.8.4D wz'ZZ be independent ofr too. 




10 



Figure 4.3: The limited region EIq^, 0) of the configuration space where the dynamics of 
the point Universe is restricted by mean of the curvature term which corresponds to an 
infinite potential well. In this region the conditions Q a > are fulfilled. 



4.8.1 The Hamilton Equations 

The main advantage relying on the reformulation of the dynamics as a chaotic scatter- 
ing process consists of replacing the discrete BKL map by a geodesic flow in a space of 
continuous variables H23I3E1-- EZE 
In terms of (|4.8.3[) , the Lagrangian (|4.7.22|) becomes 



N 



while D expresses as 



^(r 00,60), 



D = exp{-3^ r «} , 



(4.8.5) 



(4.8.6) 



and since it vanishes towards the singularity independently of its particular form, the 
only property required for T is to approach infinity in this limit. 
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The Lagrangian (|4.7.22|l leads to the conjugate momenta 



12D / T dT\ 2 



12D e 2Y , 

= -w^i^ ( 4 - 8 - 7b ) 

Pe = ^e 21 [e ~ 1) 9> (4.8.7c) 

which by a Legendre transformation make the variational principle assume the form [273J 

,-2T 



5 J Ipg + + p T r> - -^-Hj dt = 0, (4.8.8) 



where 



2 „ 2 



H = "t^tt + P? 2 (£ 2 " !) + #-r + 24ye2r - ( 4 - 8 - 9 ) 



At) 



4.8.2 Dynamics in the Reduced Phase Space 

Solving the super-Hamiltonian constraint we get the expression for "Kadm 



~Pr= j^-adm = ^Ve 2 + 2We^, (4.8.10) 
where 

^=(?-l)pf + J^. (4.8.11) 
In terms of this constraint, the principle (|4.8.8|) reduces to the simpler form 

5 J (pg' + p e 6' - I'-Kavm) dt = 0, (4.8.12) 



whose variation provides the Hamilton equations for £' and 9' 11273 

r" 



= 7T7 (r-l)PC (4.8.13a) 



0' = F/ ^ 

From (|4.8.7a|l we find the time-gauge relation 



(4.8.13b) 



N ADM (t) = ^-^r' , (4.8.14) 

thus our analysis remains fully independent of the choice of the time variable until the 
form of T and x' is fixed. 
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By choosing dT / dr = 1, the principle (|4.8.12[) reduces to the two-dimensional one [273] 

3 f [pg + p e e' - H ADM ) dt = 0, (4.8.15) 

where, remebering (|4.5.5|) 



<K ADM = Ve 2 + U, U = 2We 2T ; (4.8.16) 

moreover, the choice r' = 1 for the temporal gauge provides the lapse function as 

12D 

N A dm(t) = - e 2 \ (4.8.17) 

•H-ADM 

The reduced principle (|4.8.15|) finally gives the Hamilton equations 

tt 1 - 1) 



^ADM 

1 Ve 

^ADM (^ 2 — 1) 



2 Pe 2 



^ADM 
1 

Z'Hadm d8 



dU 



(4.8.18b) 
(4.8.18c) 



V'e = -^— (4.8.18d) 



where, because of the choice of the time gauge, ( )' = -j- . 

(A. L 



4.8.3 Billiard Induced from the Asymptotic Potential 

The Hamilton equations are equivalently viewed through the two time variables T and T 
and then, for this Section only, we choose the natural time gauge X 1 = 1, in order to write 
the variational principle (I4.8.12I) in terms of the time variable T as 



dO \ 

^-M ADM )dT = 0. (4.8.19) 



Nevertheless, for any choice of time variable T (i.e. r = t), there exists a corresponding 
function T (t) (i.e. a set of MCI variables leading to the scheme (|4.8.19l) ) defined by the 
(invertible) relation 

dr 12D v ; v ' 

The asymptotically vanishing of D is ensured by the Landau-Raichaudhury theorem 
near the initial singularity (which occurs by convention at T = 0, where T now denotes 
the synchronous time, i.e. dT = —N (t) dr), for T — > we have d\nD/dT > 0. In terms 
of the adopted variable t we have 

D^O => T(t)^oo, (4.8.21) 
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consequently, by (|4.8.6|) and (|4.8.20|) , also 

dmD _ d\nDdT _ rflnD 
dr dT dr dT 



N (t) (4.8.22) 



and therefore D monotonically vanishes even in the generic time gauge as soon as dT /dr > 
for increasing T, according to (|4.8.20|l . 

Approaching the initial singularity the limit D — > for the Mixmaster potential (|4.5.5[) 
implies an infinite potential well behavior, as discussed in Section 1431 

Furthermore, by (|4.8.10|) the important relation 



(4.8.23) 



dt ~ dt dT ~ dT 

holds, i.e. explicitly 



dY ZKadia V sr 

In this reduced Hamiltonian formulation, the functional T(t) simply plays the role of a 
parametric function of time and we recall how actually the anisotropy parameters Q a are 
functions of the variables £, 6 only (see (14.8.4I )). 

Therefore in the dinamically allowed domain TIq (see Fig |4.3l) the ADM Hamiltonian 
becomes (asymptotically) an integral of motion 

d^ADM _ n _ dE 

V{£,0}en Q { dT ~ dT (4.8.25) 
"Kadm = Ve 2 + 24 U ^e = E = const. . 



4.8.4 The Jacobi Metric and the Billiard Representation 

Since above we have shown that asymptotically to the singularity (r — > oo, i.e. a. — > — oo) 
d^ADM I dT = 0, i.e. 'Kadm = e = E = const., the variational principle l|4.8.19|) reduces to 

5 J {-p^ + p e d6 - EdT) = S J (p^dg + p d6) = 0, (4.8.26) 

where we dropped the third term on the left-hand side since it behaves as an exact differ- 
ential. 

By following the standard Jacobi procedure UTTI to reduce our variational principle to 
a geodesic one in terms of the configuration variables x a , we set x a> = g ab pb, and by the 
Hamilton equation (|4.8.13D we obtain the metric II273H276M 

8^ = j(e-l), g 66 = JgFZJ- (4-8-27) 
By (|4.8.27|l and by the fundamental constraint relation obtained rewriting (|4.8.11|) as 

(^ 2 -l)p ? 2 + ^ T = E 2 , (4.8.28) 



64 



4.8 Misner-Chitre-like variables 



we get 



al V 
gab* * 



r 

E 



Using the definition 

equation (|4.8.29|) is rewritten as 



dx" ds 
ds dt 



gabU a u h I 



which leads to the key relation 



dt 



(di. 



dt ' 
Y'E, 



Y'E. 



(4.8.29) 
(4.8.30) 

(4.8.31) 



Y'E 



Indeed the expression (I4.8.32D together with pgg' 
tional principle (|4.8.26|) in the geodesic form 



ds . (4.8.32) 
pe9' = ET' allows us to put the varia- 



■ Jl'Edt = 8 J ^Jg ab u a u h Y'E ds = 5 J \J G ab u"u b ds = 0, (4.8.33) 



where the metric G a \, 
therefore 



T'Egab satisfies the normalization condition G a i,u a u = 1 and 



ds_ = , ds_ 
dt dT 



(4.8.34) 



where we take the positive root since we require that the curvilinear coordinate s increases 
monotonically with increasing values of r, i.e. approaching the initial cosmological sin- 
gularity. 

Summarizing, the dynamical problem in the region YIq reduces to a geodesic flow on a 
two-dimensional Riemannian manifold described by the line element |9l l273| 



ds 2 



i)de 2 



(4.8.35) 



The above metric has negative curvature, since the associated curvature scalar is R = 
— 2 /e 2 ; therefore the point-Universe moves over a negatively curved bidimensional space 
on which the potential wall (|4.5.6|) cuts the region YIq, depicted in Figure 1431 
By a way completely independent of the temporal gauge we provided a satisfactory rep- 
resentation of the system as isomorphic to a billiard on a Lobachevsky plane Ifl7| . 

From a geometrical point of view, the domain defined by the potential walls is not 
strictly closed, since there are three directions corresponding to the three corners in the 
traditional Misner picture [368, 367] from which the point universe could in principle 
escape (see Fig[ 



However, as discussed in Section 14.21 for the Bianchi models, the only case in which 
an asymptotic solution of the field equations has this behaviour corresponds to having 
two scale factors equal to each other (i.e. 6 = 0); but, as shown by [50J, these cases are 
dynamically unstable and correspond to sets of zero measure in the space of the initial 
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conditions. Thus, it has no sense to speak of a probability to reach certain configurations 
and the domain is de facto dynamically closed. 

The bounces (billiard configuration) against the potential walls together with the geodesic 
flow instability on a closed domain of the Lobachevsky plane imply the point-Universe 
to have stochastic features, with a formalism true for any Bianchi type model. Indeed 
the types VIII and IX are the only Bianchi models having a compact configuration space, 
hence the claimed compactness of the domain bounded by the potential walls guaran- 
tees that the geodesic instability is upgraded to a real stochastic behaviour fl~8l 11501 [3841. 
On the other hand, the possibility to deal with a stochastic scattering is justified by the 
constant negative curvature of the Lobachevsky plane and therefore these two notions 
(compactness and curvature) are both necessary for our considerations |509, 60 J. 



4.9 The Invariant Liouville Measure 

Here we show how the derivation of an invariant measure for the Mixmaster model (per- 
formed by (1241 13031 1379 1 within the framework of the statistical mechanics) can be ex- 
tended to a generic time gauge [272J (more directly than in previous approaches relying 
on fractal methods as in I152l ll5l| ) provided the Misner-Chitre-like variables used so far. 
We have seen how the (ADM) reduction of the variational problem asymptotically close 
to the cosmological singularity permits to modelize the Mixmaster dynamics by a two- 
dimensional point-Universe randomizing in a closed domain with fixed "energy" (just 
the ADM kinetic energy) (|4.8.25|) ; the key point addressed here is that we consider an 
approximation dynamically induced by the asymptotic vanishing of the metric determi- 
nant. 

From the statistical mechanics point of view[278], such a stochastic motion within the 
closed domain TIq of the phase-space, induces a suitable ensemble representation which, 
in view of the existence of the "energy-like" constant of motion, has to have the nat- 
ural features of a microcanonical one. Therefore, the stochasticity of this system can be 
described in terms of the Liouville invariant measure 

dq = const x 8 (E — e) d^dOdp^dpo (4.9.1) 

characterizing the microcanonical ensemble, having denoted by 8 (x) the Dirac function. 
The particular value taken by the constant e (e = E) cannot influence the stochasticity 
property of the system and must be fixed by the initial conditions. This useless informa- 
tion from the statistical dynamics is removable by integrating over all admissible values 
of e. Introducing the natural variables (e, cp) in place of (p%, pg) by 



Pt = — f ==co S( p, p e = e^?-lsm<p, < <f> < In (4.9.2) 

the integration removes the Dirac function, leading to the uniform and normalized in- 
variant measure |272| 

1 

dji = d^d9d(p-^ . (4.9.3) 

The approximation on which our analysis is based (i.e. the potential wall model) is reli- 
able since it is dynamically induced, no matter what time variable T is adopted. Further- 
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more, such invariant measure turns out to be independent on the choice of the temporal 
gauge, as shown in II272| . 

The use of the invariant measure and of the Artins theorem |20| provides the complete 
equivalence between the BKL piece-wise description and the Misner-Chitre continuous 
oneim. 

According to the analysis presented in [379J, by virtue of (|4.8.13|) and (|4.8.34|) , the a- 
symptotic functions £ (T) , 6 (T) , <p (T) during the free geodesic motion are governed by 
the equations 



|U V e 2 -]cos(/) (4.9.4a) 



dd sin <p 



d(p £ sin cp 

whose parametric solution £ (T) has the form 



2 



(4.9.4b) 
(4.9.4c) 



sin (p 

T(cp) = larctanhfl p2 + fl2c ° s2 ^ +& (4.9.5b) 
yrj 2 \2 apcosf J K ' 



a 2 + sin <p a,b = const. G 1R . 

However, the global behaviour of £ along the whole geodesic flow is described by the 
invariant measure 04.9.31 ) and therefore the temporal behavior of T (t) acquires a stochas- 
tic character: if we assign one of the two functions T (t) or N (t) with an arbitrary analytic 
functional form, then the other one will exhibit a stochastic behaviour. Finally, by retain- 
ing the same dynamical scheme adopted in the construction of the invariant measure, we 
see how the one-to-one correspondence between any lapse function N (t) and the associ- 
ated set of MCI variables (I4.8.3D guarantees the covariance with respect to the time-gauge. 



4.10 Chaos covariance 

We have discussed the oscillatory regime whose properties characterize the behavior of 
the Bianchi types VIII and IX cosmological models in the BKL formalism |[54l 1571 13671 near 
a physical singularity, in which it is outlined the appearance of chaotic properties [402J : 
firstly, the dynamics evolution of the Kasner exponents characterized the sequence of the 
Kasner epochs, each one described by its own line element, with the epochs sequence 
nested in multiple eras. Secondly, the use of the parameter u and its relation to dynamical 
functions offered the statistical treatment connected to each Kasner era, finding an appro- 
priate expression for the distribution over the space of variation: the entire evolution has 
been decomposed in a discrete mapping in terms of the rational /irrational initial values 
attributed to u. 

The limits of this approach essentially reside in the non-continuous evolution toward the 
initial singularity and the lack of an assessment of chaoticity in accordance with the indi- 
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cators commonly used in the theory of dynamical systems, say in terms of the estimate of 
Lyapunov exponents. 

A wide literature faced over the years this subject in order to provide the best possible 
understanding of the resulting chaotic dynamics |49l[64| . 

The research activity developed overall in two different, but related, directions: 

(z) on one hand, the dynamical analysis was devoted to remove the limits of the BKL 
approach due to its discrete nature (by analytical treatments |33l 11241 11051 l82l 11521 
HMjESES) and by numerical simulations Il443ll64ll67[l74ll74ll. 

(if) on the other one, to get a better characterization of the Mixmaster chaos (especially 
in view of its properties of covariance |198H196ll480l l258l. 

The first line of investigation provided satisfactory representations of the Mixmaster dy- 
namics in terms of continuous variables [87J, mainly studying the properties of the BKL 
map and its reformulation as a Poincare one [32 j. 

In parallel to these studies, detailed numerical descriptions have been performed with 
the aim to test the precise validity of the analytical results 1443ll64ll75ll7llll56| . 

The efforts 1|82T , I1481[443H to develop a precise characterization of the chaoticity observed 
in the Mixmaster dynamics found non-trivial difficulties due to the impossibility, or in the 
best cases the ambiguity, to apply the standard chaos indicators to relativistic systems. 
However, the chaotic properties summarized so far were questioned when numerical 
evolution of the Mixmaster equations yielded zero Lyapunov exponents |104U258l 12031. 
Nevertheless, exponential divergence of initially nearby trajectories was found by other 
numerical studies yielding positive Lyapunov numbers. This issue was understood when 
in |64| and £198 j numerically and analytically was shown how such calculations depend 
on the choice of the time variable and parallely to the failure of the conservation of the 
Hamiltonian constraint in the numerical simulations by [514], although was still debated 
by 12581. 

In particular, the first clear distinction between the direct numerical study of the dy- 
namics and the map approximation, stating the appearance of chaos and its relation with 
the increase of entropy, has been introduced by [105J. The puzzle consisted of simulations 
providing even in the following years zero (see for example fl!96| ) Lyapunov numbers, 
claiming that the Mixmaster Universe is non-chaotic with respect to the intrinsic time (as- 
sociated with the function a introduced for the Hamiltonian formalism) but chaotic with 
respect to the synchronous time (i.e. the temporal parameter t). The non-zero claims [480 J 
about Lyapunov exponents, using different time variables, have been obtained reducing 
the Universe dynamics to a geodesic flow on a pseudo-Riemannian manifold: on average, 
local instability has been discussed for the BKL approximations. Moreover, a geometrized 
model of dynamics defining an average rate of separation of nearby trajectories in terms 
of a geodesic deviation equation in a Fermi basis has been interpreted for detection of 
chaotic behavior as a principal Lyapunov exponent. A non-definitive result was given: 
the principal Lyapunov exponents result always positive in the BKL approximations but, 
if the period of oscillations in the long phase (the evolution of long oscillations, i.e. when 
the particle enters the corners of the potential) is infinite, the principal Lyapunov expo- 
nent tends to zero. 

For example, the author of |64~| reports the dependence of the Lyapunov exponent on 
the choice of the time variable. Through numerical simulations, the Lyapunov exponents 
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were evaluated along some trajectories in the (/3 + , j8_) plane for different choices of the 
time variable, more precisely T (BKL), O (Misner) and A, the "mini-superspace" one, i.e. 
d\ = | — + p+ + p 2 - | 1 ^ 2 dr. It is shown that the same trajectory giving zero Lyapunov 
exponent for T or Q-time, fails for A. 

Such contrasting results can find a clear explanation realizing the non-covariant nature 
of these indicators and their inapplicability to hyperbolic manifolds. The existence of 
such difficulties prevented, up to now, to say a definitive word about the general picture 
concerning the covariance of the Mixmaster chaos, with particular reference to the possi- 
bility of removing the observed chaotic features by a suitable choice of the time variable, 
apart from the indications provided by fll52H15l|| . 

Interest in these covariance aspects has increased in recent years in view of the contra- 
dictory and often dubious results that emerged. 

The confusion which arises regarding the effect of a change of the time variable in this 
problem depends on some special properties of the Mixmaster model when represented 
as a dynamical system, in particular the vanishing of the Hamiltonian and its non-positive 
definite kinetic term (typical of a gravitational system). These peculiar features prevent 
the direct application of the most common criteria provided by the theory of dynamical 
systems for characterizing chaotic behavior (for a review, see [258]). 

Although a whole line of research opened up |418t |482J, the first widely accepted in- 
dications in favor of covariance were derived with a fractal formalism by |152l fl51l (see 
also |387J). Indeed, the requirement of a complete covariant description of the Mixmaster 
chaoticity when viewed in terms of continuous dynamical variables, due to the discrete 
nature of the fractal approach, leaves this subtle question open and prevents a general 
consensus in this sense from being reached. 



4.10.1 Invariant Lyapunov Exponent 

In order to characterize the dynamical instability of the billiard in terms of an invariant 
treatment (with respect to the choice of the coordinates £, 6), let us introduce the following 
(orthonormal) tetradic basis [273] 

1 sin</> 



„ _ (-I^ sin *, |-^£_) . ( 4,0,b, 

Indeed, the vector v l is nothing else than the geodesic field, i.e. it satisfies 

£_£ +lt *._0, (4,0.2, 
while the vector w* is parallely transported along the geodesic, according to the equation 

^ = ^ + r> V=0, (4.10.3) 
ds ds 
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where T l kl are the Christoffel symbols constructed by the reduced metric (|4.8.35|) . Project- 
ing the geodesic deviation equation along the vector w l (its component along the geodesic 
field v 1 does not provide any physical information about the system instability), the cor- 
responding connecting vector (tetradic) component Z satisfies the following equivalent 
equation 

w-i- (4 - m4) 

This expression, as a projection on the tetradic basis, is a scalar one and therefore com- 
pletely independent of the choice of the variables. Its general solution reads as 

Z (s) = dei + c 2 e~i , c 1/2 = const. , (4.10.5) 

and the corresponding invariant Lyapunov exponent [343J is defined as [409J 

-(z 2 +(f) 2 ) 
X v = sup lim — ^ '— , (4.10.6) 

s^oo 2S 

which, in terms of (|4.10.5|> , takes the value |273| 

A D = ^>0. (4.10.7) 

The limit (|4.10.6|) is well defined as soon as the curvilinear coordinate s approaches oo. 
In fact, from (|4.8.34|) we see that the singularity corresponds to the limit r — > oo, and this 
implies s — > oo. 

When the point-Universe bounces against the potential walls, it is reflected from a 
geodesic to another one, thus making each of them unstable. Though with the limit of 
our potential wall approximation, this result shows that, independently of the choice of 
the temporal gauge, the Mixmaster dynamics is isomorphic to a well described chaotic 
system. Equivalently, in terms of the BKL representation, the free geodesic motion cor- 
responds to the evolution during a Kasner epoch and the bounces against the potential 
walls to the transition between two of them. By itself, the positive Lyapunov number 
d4.10.7[) is not enough to ensure the system chaoticity, since its derivation remains valid 
for any Bianchi type model; the crucial point is that for the Mixmaster (types VIII and 
IX) the potential walls reduce the configuration space to a compact region (TLq), ensuring 
that the positivity of \ v implies a real chaotic behaviour (i.e. the geodesic motion fills the 
entire configuration space). 

Furthermore, it can be shown that the Mixmaster asymptotic dynamics and the struc- 
ture of the potential walls fulfill the hypotheses at the basis of the Wojtkowsky theorem 
|509tl60"l| ; this result ensures that the largest Lyapunov exponent has a positive sign almost 
everywhere^. 

Generalizing, for any choice of the time variable, we are able to give a stochastic repre- 
sentation of the Mixmaster model, provided the factorized coordinate transformation in 



4 In 1 386 1 it is shown that, given a dynamical system of the form dx/dt = F(x), the positiveness of the 
associated Lyapunov exponents are invariant under the following diffeomorphism: y = (p(x,t),dr = 
A(x, t)dt, as soon as several requirements hold, fulfilled bu the Mixmaster 1 60 1 
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the configuration space 



0+ 



OL 



-e T ^a(6,Z) 
e T{T) b + (0,£) 
e T ^b_ (0,£) , 



(4.10.8a) 
(4.10.8b) 
(4.10.8c) 



where T, a, b± denote generic functional forms of the variables x, 9, £. 

It is worth noting that this analysis relies on the use of a standard ADM reduction of 
the variational principle (which reduces the system by one degree of freedom) and overall 
because, adopting MCI variables, the asymptotic potential walls are fixed in time. This ef- 
fectively is the difference between the above approach and the one presented in [481 , 474J 
(see also for a critical analysis II102| ), though in those works the problem of the Mix- 
master chaos covariance is faced even with respect to the choice of generic configuration 
variables. 

4.10.2 On the Occurrence of Fractal Boundaries 

In order to give an invariant characterization of the dynamics chaoticity many methods 
along the years have been proposed, but not all approaches have reached an undoubtable 
consensus. A very interesting one, relying on techniques considering fractality of the 
basin of initial conditions evolution has been proposed in 1997 in [152J and opened a 
whole line of debate. The conflict among different approaches has been tackled by using 
an observer-independent fractal method, though leaving some questions open about the 
conjectures lying at the basis of it. 

The asymptotic behaviour towards the initial singularity of a Bianchi type IX trajectory 
depends on whether or not we have a rational or irrational initial condition for the param- 
eter u in the BKL map. In such a scheme, the effect of the Gauss map has been considered 
together with the evolution of the equations of motion, in order to "uncover " dynamical 
properties about the possible outcoming configurations with the varying of the corre- 
sponding initial conditions. 

Nevertheless, such approach led to some doubts regarding the reliability of the method 
itself. In fact, let us observe that rational numbers initial conditions are dense and yet 
constitute a set of zero measure and moreover correspond to fictitious singularities [54, 
367. 368 J. The nature of this initial set needs to be compared with the one regarding the 
complete set of initial conditions, with finite measure over a finite interval: the conclusions 
obtained after the dynamical evolution are not necessarily complementary between the 
two initial assumptions. 

In ||151U152T) , Cornish and Levin used a coordinate-independent fractal method to show 
that the Mixmaster Universe is indeed chaotic. By exploiting techniques originally devel- 
oped to study chaotic scattering, they gained a new perspective on the evolution of the 
Mixmaster cosmology, finding a fractal structure, namely the strange repellor (see Fig l4.4|) , 
that well describes chaos. A strange repellor is the collection of all Universes periodic in 
(u, v), and an aperiodic one will tipically experience a transient age of chaos if it brushes 
against the repellor. The fractal pattern was exposed in both the exact Einstein equations 
and in the discrete map used to approximate the solution. The most important feature 
gained is a fractal approach independent of the adopted time coordinate and the chaos 
reflected in the fractal weave of Mixmaster Universes is unambiguous. 
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Figure 4.4: The numerically generated basin boundaries in the [u, v) plane are built of 
Universes which ride the repellor for many orbits before being thrown off. Similar fractal 
basins can be found by viewing alternative slices through the phase space, such as the 
(fl, p) plane. The overall morphology of the basins is altered little by demanding more 
strongly anisotropic outcomes. (From \152\ ) 



The approach used in ||152t I151| is based on the method firstly stated in [84J where it 
is shown how fractal boundaries can occur for some solutions involving chaotic systems. 
The space of initial conditions is spanned giving rise to different exit behaviours whose 
borders have fractal properties: this constitutes a conjecture as a typical property of chaotic 
Hamiltonian dynamics with multiple exit modes. 

For the case of the Bianchi IX model potential (see Figure I4.2[) the openings are obtained 
widening the three corners, on the basis that the point representing the evolution spends 
much of the time there nearby. 
This method has three essential fallacies: 

i) the case-points chosen as representatives within this framework are the ones whose 
dynamics proceeds never reaching the singularity; 

ii) the "most frequent" dynamical evolution is the one in which the point enters the 
corner with the velocity not pamllelly oriented towards the corner's bisecting line 
and, after some oscillations, it is sent back in the middle of the potential. This effect 
is altered when opening the potential corners; 

iii) the artificial opening up of the potential corners adopted in the basin boundary 
approach could be creating the fractal nature of it. 

In particular, the third observation is supported by the existence of strange attractors 
that are not chaotic, as counter-exampled by [224 J and discussed by [248 J. The choice of 
the method adopted to characterize the property of chaos or its absence is very relevant, 
especially when based on the presence of fractal boundaries in the dynamics underlying 
Bianchi IX models. This is important to be checked, first of all, because the result of 
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||152l I151]| relies on the conjecture as in [84J that opening gates in a chaotic Hamiltonian 
system can result in the presence of fractal basin boundaries (which needs in principle to 
be checked in the case of Bianchi IX), not satisfying the necessity of a general statement 
concerning chaos: even the opening of the corners does not solve the question about what 
happens when taking the limit of closing them and if there is an universal behaviour (for 
general systems). 

Secondly it is needed to integrate the Bianchi IX flow and this operation is not necessarily 
commuting with the statement regarding the remaining (and equally relevant) part of the 
set of initial conditions constituted by the irrational numbers, which needs to be checked. 

Motter and Letelier, in )|387| in the criticism to the paper of |I1521 1151| , claim the same 
results with more accurate comprehension of the global chaotic transient and afford cal- 
culations involving a more stable constraint check and a higher order integrator. Again 
the same criteria used by ||152l lT51 1 is followed to get the same results. The informations 
obtained following the Farey map approach are not relevant (only rational values of u 
are led to the three peculiar outcomes) because the corresponding invariant set contains 
almost every point of phase space. But they claim it is possible to get strict indications of 
chaos with the Hamiltonian exit method >MlH6Tl[T68lH67ll85| . 

Firstly one has to fix the width (< — > u ex j t ) of the open corners, then let the system 
evolve. The future invariant set leads to a box-counting dimension Do (estimated from 
the uncertainty exponent method B402fl ) coherent with previous results, which is, by con- 
struction, a function of the width itself. The value of Do found, equal numerically to 1.87, 
is dependent on the change done to the original potential, and converges to the value of 
2, which is an indicator of non-chaoticity ||169|l . Any of such fundamental properties, if 
outlined in a specific case, must be jointed through a limit procedure to the general one. 

Pianigiani and Yorke in [410J study the evolution of a ball on a billiard table with 
smooth obstacles so that all trajectories are unstable with respect to initial data. This 
is a system energy conserving and then they open a small hole on such table in order to 
allow the ball to go through. Such two differences have not been taken account of in the 
Bianchi IX analysis. 

In the work by |461|1 Schneider et al. it is supposed to show the existence of a chaotic 
saddle, whose signature is the chaotic basin. 

The paper by 1169] declares the absence of such points, in a model with A = 0, hence 
we infer the inapplicability of that method to discover a supposed unknown feature of a 
dynamical system. 

They stress too that the limit (not unnatural) for A — > doesn't matches: it doesn't permit 
to characterize the chaos in mixmaster vacuum model: a continuous change in a param- 
eter of the theory heavily affects the method's applicability, mainly while the study of 
Bianchi IX dynamics is of interest towards the initial singularity, where the BKL approach 
applies: in such approximation, the domain walls close to a circle. 
Hence there are objections which are subject for interesting further investigation 

(z ) has the opening of a polygonal domain the same effect as the opening a circle (which 
has curvature)? 

(ii ) Is the system truly independent from temporal reparameterization? 
(Hi ) Is the opening independent of temporal (either spatial) reparameterization? 
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(iv ) Could exist a temporal reparameterization whose effect is to close the artificial open- 
ings? 

(v ) How to interpret this eventuality? 

Even if this is not relevant for dynamical system in classical mechanics, in General Rel- 
ativity it would; for example, Cornish and Levin claim that they open a non compact 
domain, differently from Misner. 

For all the criticism here outlined we consider an analytical approach crucial to distin- 
guish among chaos indicators relying on numerical properties not well-manageable via 
numerical simulations. 



4.11 Isotropization Mechanisms 

The isotropic FRW model can accurately describe the evolution of the Universe until the 
decoupling time, i.e. 10~ 3 — 10~ 2 seconds after the Big-Bang |310| . On the other hand, 
the description of its very early stages requires more general models, like at least the 
homogeneous ones. Therefore we are interested to investigate some mechanisms al- 
lowing a transition between these two cosmological epochs. When the anisotropy of 
the Universe is sufficiently suppressed, we can speak of a quasi-isotropization of the 
model |472[ 11391 11401 1146 1 (for a detailed discussion of the isotropization mechanism see 
IlIM 11821 14351 M 11911 1^ 1471 127^ 

and such mechanism can be regarded as a "bridge" between the two stages. 

In this paragraph we will discuss the origin of a background spaced when a real self- 
interacting scalar field <p is taken into account, following the work [305 j. Let us extend the 
Misner-like variables ||302| (see (|4.5.1|) ) 

q u = A ap +e[f ^3 = $ (4A\.\) 

where ; = 1, 2, a = 1, 2, 3 and AJ satisfy 

£>» = 0, £A<A« k = 65 jk . (4.11.2) 

a a 

As usual, a parametrizes the isotropic change of the metric with the singularity appearing 
as a — > -co and the ft the anisotropics of the model. When these variables approach some 
constants as a limit we can speak of a quasi-isotropization of the model (for a discussion 
on the quasi-isotropic solution see | |325H293[ 12921 12771). 
The action expressed in terms of them writes as 



it 



Prdtf? + P a d t a - I £P r 2 + 6U 



(4.11.3) 



5 The chaotic nature of the evolution toward the singularity implies that the geometry, and therefore all the 
geometric quantities, should be described in an average sense only. With this respect, the Universe does 
not possess a stable background near the singularity 1 299 304] . 
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where r = 1,2,3, h = exp(3a) and the potential term <|4.8.16)> is now defined as U = 
h(W - (3) R) with W((p) = y2[h'idi(pdj(p + V(<p)]. From the action d4.ll.3ll an inflationary 
solution comes out imposing the constraint 



-U ~ V(d>) ~ const. > (3) R 
h 



(4.11.4) 



which can be realized by an appropriate process of spontaneous symmetry breaking, ex- 
haustively studied in fl305l 137711 . Let us consider the situation where U = hA, where 
A = const. 

The Hamilton-Jacobi equation is 



?(f) 2 -(!) 2+6e>< p< 3 » )A = - 

whose solution can be expressed as 



(4.11.5) 



S(p,oc)~K r p r +^K a + ^]n 



K a — K 



K a + K 



where K a (K r ,oc) = ±y£r K} + 6Aexp(3a), with some generic constants K 
and K r . The equation of motion for a is readily obtained as from (|4.11.3[) 



3a: 
dt 



NP n , 



3exp(3a/2) 



(4.11.6) 



(4.11.7) 



Choosing a as the time coordinate, i.e. d t a = 1, the time gauge condition becomes N = 
—3 exp(3a/2) / P u . Since the lapse function is positive defined we must also have P a < 0. 

Accordingly to the Hamilton-Jacobi method, one has firstly to differentiate with respect 
to K r and then, equating the results to arbitrary constants, one finds the solutions describ- 
ing the trajectories of the system as 



SS 

W r 



ft 



>)=& + 



Kr 

3\K\ 



In 



K„ — K 



(4.11.8) 



where fi r are new arbitrary constants. Let us investigate the two limits of interest. First 
of all, let us note that for h — > oo (i.e. a. — > oo, _K a — > oo) the solution (|4.11.8|) transforms 
into the inflationary one obtained in M472L corresponding to the quasi-isotropization of 
the model as the functions $ r approach the constants f$ r Q . On the opposite limit, i.e. for 
h — > (a — > — oo, K a — > K), the solution d4.ll.8l > provides the generalized Kasner one as 
expected, simply modified by the presence of the scalar field 



P r («) = Po 



K Tf 



(4.11.9) 



txo being the remaining constants. 
The existence of the solution (|4.11.8)> shows how the inflationary scenariol 230U2311 13271 



462J can provide the necessary dynamical "bridge" between the fully anisotropic and the 
quasi-isotropic epochs of the Universe evolution. In fact, during that time the anisotropies 
/3± are dumped and the only effective dynamical variable is a, i.e. the one related to the 
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isotropic volume of the Universe. This shows how the dominant term during the inflation 
is Ae 3a and any term in the spatial curvature becomes more and more negligible although 
increasing like (at most) e 2a . 

For a sample of works dealing with the Bianchi models dynamics involving a cosmo- 
logical term, see II4M 15071 ITM[Tl5ll4T5ll4T4ll45^ 



4.12 Cosmological Implementation of the Bianchi Models 

In this paragraph we shall examine the cosmological issues of the Bianchi models. In 
particular, we will focus attention on the question about the corresponding theoretical 
predictions regarding the relic Cosmological Background Radiation (CMB) anisotropy 
Ifl66ll323l . Confirmed observations 112391 14171 1469ll342ll392ll show that the large-scale relic 
radiation anisotropy AT/T is about 

AT/T <2-l(T 5 , (4.12.1) 

in terms of the equivalent black body radiation temperature T. Therefore the Universe 
becomes transparent to the relic radiation at an epoch when the expansion anisotropy is 
"small". A comparison between the theory and the experimental observations requires 
the analysis of a quasi-isotropic stage of the Bianchi models. 

The discussion which follows is mainly based on the textbook fl516| and references 
therein (see also |47ll). Initially we approach the Bianchi I model and then we focus 
attention on the more general Bianchi IX model. 



4.12.1 Expected Anisotropy of CMB for Bianchi I 

Let us analyze the Bianchi I model, characterized by a flat, comoving three-dimensional 
space. 

We will include in the dynamics an ordered magnetic field having oriented fluxes of 
relativistic particles since they inevitably develop as a consequence of the processes dur- 
ing the early stages of anisotropic expansion. Let the magnetic field be oriented along 
the z-axis, W be its energy density and p x , p y , p z components of the pressure of the free 
particles along the corresponding axes. Let us also assume e> W and e 3> p x , Py, p z and 
moreover a~ 1 da/dt ~ b~ 1 db/dt ~ c~ 1 dc/ dt, e being the energy density of ordinary mat- 
ter. Such relations ensure that, to first approximation, the expansion isotropically takes 
place and thus the equations of motion read as 



J_d_ 

abc dt 
1 d 
abc dt 



da 
li 
da 
~di 



b — abc 
dt J 

■ c~ x %- | abc 
dt J 



K (Vx ~ Py) 
K[{p x -p z )+2W] 



(4.12.2a) 
(4.12.2b) 



Suppose that P is the pressure associated to e and consider the particular case of ultra- 



relativistic matter P = e/3 (corresponding to a ~ b ~ c ~ t ' ). In this case, we can 
derive the behavior of the deformation anisotropies (defined as the left-hand sides of the 
following equations (|4.12.3|l ) when the expansion is almost periodic and (|4.12.2|) rewrite 
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as 



a 1 da/dt — b 1 db/dt 
a~ 1 da/dt 

ar x dal dt — c~ 1 dc/dt 
a~ 1 da/dt 



e \ t 



HPx-Pz) { 6W | ft* 
e e V t 



1/2 



(4.12.3a) 
(4.12.3b) 



t* being a constant. In absence of matter the anisotropy decays as £~ 1//2 . On the other 
hand, when the expansion is almost periodic, the ratios (p x — Py)/e as well as W/e re- 
main constant. Therefore, in presence of a magnetic field or of a flux of relativistic parti- 
cles, the deformation anisotropy is conserved during the era when P = e/3. In conclu- 
sion, the anisotropy of the stress-energy tensor slows down the isotropization of solution. 
This statement can be generalized for all (reasonable) kinds of matter and in particular 
for the case P = 0. 

Let us discuss the equation for the anisotropy of the relic radiation temperature and 
how it is related to the expansion anisotropy, and consider an observer receiving the radi- 
ation from the direction of the two axes with scale factors a and b. The observed difference 
of temperature in those two directions is 



AT 

T" 



f— »0O 



(4.12.4) 



21 



where Z\ is the redshift associated to the time when the Universe is transparent to the 
radiation and the radiation field is isotropic^. 

If the sharply anisotropic stage ends before the time when the radiation density p r 
equates the baryon one p m , then the deformation anisotropy (|4.12.3a[) is conserved. After 
this stage, in the P = epoch the anisotropy decreases according to 



6p 



ams .-2/3 



a p m 

Panis being the density of the anisotropic neutrino flux (p an is/Pr 
transparency redshift as Z\ = 8, one obtains 



(4.12.5) 

0.1 — 1). Assuming the 
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(4.12.6) 



Hence, expression (|4.12.6|) gives the expected value for the relic radiation anisotropy al- 
though it is valid only in the case p m = p cr n (flat space). 



4.12.2 Expected Anisotropy of CMB for Bianchi IX 

In the homogeneous cosmological models with curved comoving space, as Bianchi IX 
is, the deformation anisotropy behaves in a way similar to the one described above. In 
fact, it is constant during the radiation dominated era up to the time when p r = p m . 
Nevertheless, the difference relies in the role played by the non-interacting particles in 
Bianchi I, which actually is taken by the spatial curvature. Then, after this epoch, the 

6 After this time the radiation freely propagates. 
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anisotropy decreases as £ /3 , while during the time when p r = p m the matter becomes 
transparent to the background radiation and the photons freely propagate. 

Therefore, let us underline two main aspects of the relic radiation anisotropy in the 
anisotropic models. The first feature is related to the conservation of the deformation 
anisotropy during the radiation dominated era, and it is responsible for the amplitude of 
AT/T. The second characteristic is connected to the curvature of the comoving space and 
the consequent motion of the photons which gives rise to the different angular distribu- 
tions on the sky of the background anisotropy. 

The anisotropy of the relic radiation temperature can now be calculated with the same 
arguments given above for Bianchi I and it depends on three characteristic temporal steps: 

i) t<p, the time of the beginning of the isotropic stage, when a ~ b; 

ii) t c , when the equation of the state changes, i.e. when it effectively becomes P = 
and after this time all the deviations from exact isotropy decay as a power law; 

iii) t e , when matter becomes transparent. 

Explicitly, the anisotropy of the temperature background reads as B516|) 



It is worth noting that the above formula (|4.12.7|) is alarmingly unstable to small numerical 
changes in underlying parameters. This feature is widely discussed in [36| where an 
analysis of the cosmological evolution of matter sources with small anisotropic pressures 
is performed (see also |37| for applications to the magnetic case and [35] in which a more 
general analysis is shown). Furthermore the dependence of (J4.12.7|) on the isotropization 
time t,p is extremely weak. In fact, taking t c ~ t e and assuming AT/T ~ O(10~ 3 ), one 
has ttp ~ 0(fp), where tp ~ O(10 _44 s) is the Planck time. Therefore isotropization would 
take place in the region when the applicability of classical General Relativity is expected 
to fail. On the other hand, AT/T is of order of the values experimentally observed only 
if the factor (t c / t e ) 2 ^ 3 is small, i.e. if the matter becomes transparent much later than the 
end of the radiation dominated epoch (t e 3> t c ) and this is possible only if the amount of 
ionized intergalactic gas is so relevant that the critical parameter is Q ~ 0(1)- 

With respect to the angular distribution, we stress that the light propagates along the 
principal direction of the deformation tensor and therefore the anisotropy has a quadrupole 
character. 

As last point, we have to stress that in above analysis no mention to inflation is made. 
Indeed it would a major damping effect on any primordial anisotropy, massively reducing 
it below (14.12.71) though. About this point we refer to I36ll37ll35|. 



In this Section we face the influence of a scalar field when approaching the cosmological 
singularity. As shown in the initial studies p2ll53| . such field can suppress the Mixmaster 
oscillations during the evolution toward the singularity. 




(4.12.7) 



4.13 The Role of a Scalar Field 
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Following the approach given in [69], let us consider the Mixmaster Universe in the pres- 
ence of a self-interacting scalar field (p. As we have seen previously, the Einstein equations 
are obtained from the variation of the Hamiltonian constraint H = 0, where (N oc e 3a ) 

H = H K + H V , (4.13.1) 

Hk and Hy being the kinematic and potential part of the Hamiltonian, respectively. In 
particular we have 

Hk = -pI + P 2 + + P 2 - + Pl' (4- 13 - 2 ) 

and 



H v = e 4a 



- 2 (>+ + e^- 2 ^- + ^++2^-)] + (4^3.3) 

where with V((p) we denote a generic potential of the scalar field. Working with such 
Misner variables the cosmological singularity appears as a — > — oo. Therefore, unless 
V(<p) contains terms exponentially growing with a, the very last term in (|4.13.3|) can be 
neglected at early times, i.e. e 6a V((p) —> as a —> — oo. 

Let us consider the kinematic part Hk- Its variation yields equations whose solution 
reads as 

J>± = fi° ± + v ± \a\ 



b° + v 4) \dl\, 



(4.13.4) 



where v± = p±/\p a \ and = p<p/\p K \- Therefore the constraint Hk = becomes 

v\ + v 1 _+v\ = \. (4.13.5) 

Let us firstly analyze the case without the scalar field, i.e. (p = 0. Introducing polar 
coordinates in the anisotropy planes as v + = cos 6 and y_ = sin 9, through equations 
H4.13.4fr the potential (|4.13.3fr rewritee as 

H V ~ e - 4 l a l( 1 + 2 cos0) + e -4|«|(l-cose-V3slne) + g -4|*|(l-cose+%/3sine)^ (4.13.6) 

where we maintained the dominant terms only, i.e. the first three terms in (|4.13.3fr . Except 
for the set of zero measure of values 9 = (0,27r/3, 47r/3), any generic value of 9 will cause 
the growth of one of the terms on the r.h.s. of (|4.13.6fr . 

Let us consider the case cp / and hence v 2 , > 0. Equation (|4.13.5fr is replaced by 

v\ + v 2 _ = l-vl<l, (4.13.7) 

thus none of the terms in fl4.13.6fr will grow if the following conditions are satisfied 

'l + 2y+ > 

l-v+- V3v_ > (4.13.8) 
1 -v+ + y/3v- > 0, 
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situation which is realized if v\ < 1/2 and v 2 _ < 1/12, which occur if 2/3 < vi < 1. In 
[74J it is described how p a decreases at each bounce, and therefore for any initial value of 
Pip it is always vi > 2/3. 

As we have seen, the approach to the singularity of the Bianchi IX model is described 
by a particle moving in a potential with exponentially closed walls bounding a triangular 
domain. During the evolution, the particle bounces against the walls providing an infinite 
number of oscillations toward the singularity. The scalar field influences such dynamics 
so that for values of v± satisfying (|4.13.8|) , there are not further bounces and the solution 
approaches (|4.13.4|) . In other words there will be an instant of time after which the point- 
Universe does not ever reach the potential walls and no more oscillations appears. In this 
sense the scalar field can suppress the chaotic Mixmaster dynamics toward the classical 
cosmological singularity. 



4.14 Multidimensional Homogeneous Universes 



When the number of spatial dimensions is greater than three, homogeneous models can 
loose their chaotic dynamics, as already arises in the four-dimensional case. 
The question of chaos in higher dimensional cosmologies has been widely investigated, 
and many authors f42l 12061 12381 [237 1 showed that none of higher-dimensional extensions 
of the Bianchi IX model possesses proper chaotic features and the crucial difference is 
given by the finite number of oscillations characterizing the dynamics near the singularity. 
Without loss of generality we will follow the analysis proposed by Halpern[238, 237J and 
limit our discussion to the case of a homogeneous model with four spatial dimensions. 

The work of Fee[195] classifies the four-dimensional homogeneous spaces in 15 types, 
named GO — G14, and it is based on the analysis of the corresponding Lie groups. The 
line element can be written using the Cartan basis of left-invariant forms and explicitly 
reads as (N = 1) 

ds 2 = dt 2 - \ s {t)co r ® a; 5 . (4.14.1) 
1 

The 1-forms of obey the relation dco' = -C pq co p A co q , where the C r pq are the four- 
dimensional structure constants. The discussion remains quite general even limiting it 
to the case of a diagonal matrix 4 ^ rs 

% s = dia S {a 2 ,b 2 ,c 2 ,d 2 ) . (4.14.2) 
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The Einstein equations are obtained with the standard procedure as 



R° 



Rl 



4 



R 



ii b 
a b 


c d 
h~ + -j = 0, 
c d 


(4.14.3a) 


{abed) 


+ s\ = o, 


(4.14.3b) 


{abed) 
abed 


+ S 2 2 = , 


(4.14.3c) 


{abed) 
abed 


+ S 3 3 = 0, 


(4.14.3d) 


{abed) 
abed 


+ s 4 4 = , 


(4.14.3e) 


\x n 


c m = o • 


(4.14.3f) 



where x n (n = 1,2,3,4) denote the scale factors a,b,c,d, respectively, and the S n n are 
functions of them and of the structure constants. 

The analysis developed for the standard Bianchi type IX can be straightforwardly gener- 
alized to the five-dimensional case, obtaining the system 



a TT = -A 2 S 1 l , f> TT = -A 2 S 2 2 , (4.14.4a) 
A — -A 2 ^ 4 



7tt = -A 2 S 3 3 , S TT = -A 2 S 4 



= 2a Tj 6 T + 2a T 7 T + 2oc T 3 T + 2,6 T 7 T + 2,S T £ T + 2y r S r , (4.14.4b) 

and 

R° n = . (4.14.4c) 

The dynamical scheme (|4.14.4|) is valid for any of the 15 models using the corresponding 

en 



Among the five-dimensional homogeneous space-times, G13 is the analogous of the 
Bianchi type IX, having the same set of structure constants. 
The Einstein equations can be written as 



2«TT 




{b 2 - c 2 ) 2 


-a 4 


d 2 , 


(4.14.5a) 


2/3 TT 




'{a 2 - c 2 ) 2 


-fo 4 


d 2 , 


(4.14.5b) 


27tt 




\b 2 - a 2 ) 2 


-c 4 


d 2 , 


(4.14.5c) 


5 TT 


= 0, 






(4.14.5d) 



togheter with (|4.14.4b|) . If we assume that the BKL approximation is valid, i.e. that 
the right-hand sides of equations d4.14.5l) are negligible, then the asymptotic solution for 
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-oo is the five-dimensional and Kasner-like line element 

4 



ds 2 = dt 2 ~Y J t 2pr {dx r ) 2 , (4.14.6) 



r=l 

with the Kasner exponents p r satisfying the generalized Kasner relations 

4 4 



r=l r=\ 

This regime can only hold until the BKL approximation works; however, as T approaches 
the singularity, one or more of the terms may increase. Let us assume p\ as the smallest 
index; then a = exp(a) is the largest contribution and we can neglect all other terms, 
obtaining 

1 

a TT = --exp(4a + 2£) , 

/3tt = 7tt = ^exp(4 a + 25), (4.14.8) 
5 TT = 0. 

As soon as the asymptotic limits for x — > ±oo are considered, from the solution to (|4.14.8|) 
we obtain the map 

, = Pi +P4 / = P2 + 2gi + p 4 

Pl 1 + 2pi + p 4 ' Vl 1 + 2pi + p 4 ' 

P 3+2 P1 +P4 , = VA 

l+2 Pl + p 4 F4 l + 2 Pl + p 4 ^ ^ 

abed = A't , A' = (1 + 2pi + p 4 )A . (4.14.10) 

The difference of this dynamical scheme with the four-dimensional case relies in the con- 
ditions needed to undergo a transition: analysing the behavior of the potential terms in 
(|4.14.5|) we see how two of the four parameters must satisfy the inequality 

1 - 3pl - 3p\ - 2pi p 2 + 2pi + 2p 2 > , (4.14.11) 

and one of the following ones 

3 Pl + Pl + Pl - P2 - PlP2 < , 

3p| + p? + P2"Pi -P1P2 < 0, (4.14.12) 
3 Pi + P2 - 5pi - 5p 2 + 5pip 2 + 2 < . 

Figure 14.51 shows the existence of a region where condition (|4.14.11|) is satisfied but none 
of d4.14.12D is, thus the Universe undergoes a certain number of transitions and Kasner 
epochs and eras; as soon as the Kasner indices p\, p 2 assume values in the shaded region, 
then no more transitions can take place and the evolution remains Kasner-like until the 
singular point is reached. 
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Figure 4.5: The shaded region corresponds to all of the couples (pi, p-j) not satisfying 
(14.14.12)) : as soon as a sey (pi, p-i) takes values in that portion, the BKL mechanism breaks 
down and the Universe experiences the last Kasner epoch till the singular point. 



Type G14 case is quite similar to G13: for this model, the structure constants are the 
same as Bianchi type VIII and, under the same hypotheses, only a finite sequence of 
epochs occurs. 

The results of this analysis for the five-dimensional homogeneous space-times can be 
extended to higher dimensions and reveal how chaos is a dimensional phenomenon, for the 
homogeneous case, limited to the four-dimensional space-time. 



4.15 The Role of a Vector Field 



In this section we investigate the effects of an Abelian vector field on the dynamics of 
a generic (n + 1) -dimensional homogeneous model in the BKL scheme; the chaos is re- 
stored for any number of dimensions, and a BKL-like map, exhibiting a peculiar depen- 
dence on the dimension number, is worked out |59| . These results have also been inserted 
in more general treatment by Damour and Hennaux lil64|| . 

A generic (n + 1) -dimensional space-time coupled to an Abelian vector field Au = 
(<p, Am), with a = (1, 2, . . . , n) in the ADM framework is described by the action 

S = J d n xdt (n a ^/z aj 8 + n a |;A a + (pDJT -NH- N a H^ , (4.15.1) 



where 
H 



1 



T (nS) 



i 



(4.15.2a) 
(4.15.2b) 



-V^ + IT^j, 

denote the super-Hamiltonian and the super-momentum respectively, while F a p is the 
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spatial electromagnetic tensor, and the relation D a = d a + A a holds. Moreover, TP and 
n a P are con j U g a te momenta to the electromagnetic field and to the n-metric, respec- 
tively, which result to be a vector and a tensorial density of weight 1/2, since their explicit 
expressions contain the square root of the spatial metric determinant. The variation with 
respect to the lapse function N yields the super-Hamiltonian constraint H = 0, while with 
respect to q> it provides the constraint 9 a n a = 0. 

We will deal with a source-less Abelian vector field and in this case one can consider the 
transverse (or Lorentz) components for A K and Il a only. Therefore, we choose the gauge 
conditions cp = and D a Yl a = 0, enough to prevent the longitudinal parts of the vector 
field from taking part to the action. 

It is worth noting how, in the general case, i.e. either in presence of the sources, or in 
the case of non- Abelian vector fields, this simplification can no longer take place in such 
explicit form and the terms cp(d a + A a )Il a must be considered in the action principle. 

A BKL-like analysis can be developed|59] as well as done previously, following some 
steps: after introducing a set of Kasner vectors l a and the Kasner-like expanding factors 
exp(^ fl ), the dynamics is dominated by a potential of the form X^e^Aj, where A fl are the 
projection of the momenta of the Abelian field along the Kasner vectors. With the same 
spirit of the Mixmaster analysis, an unstable n-dimensional Kasner-like evolution arises, 
nevertheless the potential term inhibits the solution to last up to the singularity and, as 
usual, induces the BKL-like transition to another epoch. Given the relation exp(q") = f", 
the map that links two consecutive epochs is 

P'^—^ V' a = V i±4^, (4.15.3a) 

X; = Ai , X; = A fl ( 1 - 2 ( n ~ 1 )Pi ^ (4.15.3b) 

1 V (n-2)p a + np 1 J 

An interesting new feature, resembling that of the inhomogeneous Mixmaster (as we will 
discuss later), is the rotation of the Kasner vectors, 

£' a =£ a + cr a £ x , (4.15.4a) 

a a = = -2 ( (n " 1)Pl h . (4.15.4b) 

Ai (« - 2) y a + npi Ai 

which completes our dynamical scheme. 

The homogeneous Universe in this case approaches the initial singularity described by 
a metric tensor with oscillating scale factors and rotating Kasner vectors. Passing from 
one Kasner epoch to another, the negative Kasner index pi is exchanged between differ- 
ent directions (for istance t\ and £2) and, at the same time, these directions rotate in the 
space according to the rule ((4.15.4b)) . The presence of a vector field is crucial because, 
independently of the considered model, it induces a dynamically closed domain on the 
configuration space|60|. 

In correspondence to these oscillations of the scale factors, the Kasner vectors J a rotate 
and the quantities a a remain constant during a Kasner epoch to lowest order in q a ; thus, 
the vanishing of the determinant h approaching the singularity does not significantly af- 
fect the rotation law ((4.15.4b)) . 
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There are two most interesting features of the resulting dynamics: the map exhibits a 
dimensional-dependence, and it reduces to the standard BKL one for the four-dimensional 
case. 
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This Section faces the question of the quantum evolution of the Mixmaster model. In par- 
ticular, after a discussion on the Wheeler-DeWitt approach and the so-called problem of 
time, the first quantization of the model is achieved in both Misner variables and MCI 
ones. Furthermore, the model is also analyzed in the Loop Quantum Gravity framework 
(with particular attention to the disappearance of chaos) and in the Generalized Uncer- 
tainty Principle one. A brief discussion of the so-called Quantum Chaos closes the Section. 
In the rest of the paper, we set h = 1. 



5.1 The Wheeler-DeWitt Equation 

In this Section we will briefly review the Wheeler-DeWitt approach to the quantum grav- 
ity formulation, in the metric formalism. For the detailed literature about this topic we 
refer to lMl3TBl l26l ll90ll244ll . 

This scheme relies on the Dirac approach to a first-class constrained system 12511 1181| , 
i.e. the quantum theory is constructed without solving the constraints. Of course this 
method carries some un-physical information that will be removed imposing some con- 
ditions to select the physical states. In particular, if G a is a first class constraint, a phys- 
ical state must remain unchanged when performing a transformation generated by G a . 
Thus, the physical states are the ones annihilated by the quantum operator constraints, 
i.e. G a |Y)=0. 

The quantization of General Relativity in a canonical formalism prescribes to imple- 
ment the Poisson algebra (|2.4.11aH2.4.11c|) in the form of the canonical commutation rela- 
tions 



haf}(x,t),h 7 s(x' r t) 
f\^{x,t),fV\x',t) 



o 





iK5l/ 5) 5\x-x'). 



(5.1.1a) 
(5.1.1b) 
(5.1.1c) 



We note that equation (|5.1.1a|) is a kind of microcausality condition for the three-metric 
field, though the functional form of the constraint is independent of any foliation of space- 
time: thus this confirms that the points of the three-manifold E are space-like separated. 

In the next step, we impose the constraint equations (|2.4.14[) as operators to select the 
physically allowed states 



H [x;h,Il) T = 
H a (x;h,n) T = 0. 



(5.1.2a) 
(5.1.2b) 
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Since the dynamics of General Relativity is fully contained in the classical constraints, we 
do not have to analyze the dynamical equations, whose meaning is readily obtained. As 
we have seen, the Hamiltonian for GR (|2.4.12|) reads as 

% = J^d 3 x(NH + N x H x ), (5.1.3) 

and therefore, considering (|5.1.2a|) and (|5.1.2b|) , in a putative Schrodinger-like equation 

i Y t = $C¥ t = (5.1.4) 
at 

the wave functional Y (also known as "the wave function of the Universe" M245I1 ) turns 
out to be independent of "time". This is the so-called "frozen formalism", because it ap- 
parently implies that nothing evolves in a quantum theory of gravity. By other words, 
we have an identification of the quantum Hamiltonian constraint as the zero-energy 
Schrodinger equation "KY = 0. This is known as the problem of time and has deep 
consequences on the interpretation of the wave function of the Universe and deserves to 
be treated in some details in Section EH 

When writing the expression l|5.1.3|) , we assumed the primary constraints 

c(x,t) = n(x,t) = o, c a (x,t) = n x (x,t) = o, (5.1.5) 

implemented at a quantum level and therefore the wave functional Y = Y(h K p,N,N a ) 
becomes function of the three-metric only, i.e. Y = Y(/z a ^). 

Let us now explicitly discuss the meaning of constraints (|5.1.2a|) and (|5.1.2b|) . First of 
all, a representation of the canonical algebra can be chosen to be 

-~- 5Y 

h afi Y = h a *Y, LT^Y = -ix— . (5.1.6) 

This is the widely used representation of the canonical approach to quantum gravity in 
the metric formalism. However, the above equations do not define proper self-adjoint 
operators because of the absence of any Lebesgue measure on E [280 j. 

Let us firstly address the constraint (|5.1.2b|) , which is the so-called diffeomorphisms 
one, because the wave functional Y(/z a ^) depends on a whole class of three-geometries 
{h a p} (invariant under three-diffeomorphisms) and not only on the three-metric, i.e. Y = 
Y({/i a ^}). Therefore the configuration space for quantum gravity will be the Wheeler 
superspace|370J. In literature this is referred to as the "kinematical constraint". The 
dynamics is generated via the scalar constraint (|5.1.2a[) , providing the famous Wheeler- 
DeWitt equation fT77l [1781 IT79J , which explicitly reads as 

H(x)Y = -% p 5 {x) /'J - WRY = 0, (5.1.7) 

where 9a^s is the supermetric (|2.4.10c)> . This equation is at the heart of the Dirac con- 
straint quantization approach and the key aspects of the canonical quantum gravity are 
all connected to it. 

There are several problems in the WDW approach to quantum gravity, both mathemat- 
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ical and conceptual ones, for which we refer to a wide literature [315 , 280 . 485 . 486J on the 
topic. We will discuss the two most relevant among them: 

i) since equation fl5.1.7| ) contains products of functional differential operators evalu- 
ated at the same spatial point, it is therefore hopelessly divergent; moreover the 
distributions in the denominator are not clearly defined. 

ii) Understanding the physical meaning of the WDW equation requires a clarification 
of the notion of time at a quantum level, and eventually this is related to the classical 
(3 + l)-splitting performed before the quantization procedure 13801 1362ll44| . 

5.2 The Problem of Time 

A major conceptual problem in quantum gravity is the issue of what time is and how it 
has to be treated once a formalism is adopted (for a detailed discussion see [280], while 
for the role of conformal three-geometries we refer to [513]). This task is deeply connected 
with the special role assigned to temporal concepts in all theories of physics different from 
GR. For example, in Newtonian physics, as well as in non relativistic quantum mechanics, 
time is an external parameter to the system itself and is treated as a background degree of 
freedom. In ordinary quantum field theory the situation is similar since the Minkowski 
background is fixed and the Newtonian time is replaced by the time measured in a set 
of relativistic inertial frames. Such notion of "time" plays a crucial role in the conceptual 
foundations of the quantum theory. In fact, in the conventional Copenhagen interpreta- 
tion of quantum mechanics, an observable is a quantity whose value can be measured 
at fixed time. Moreover, the scalar product is conserved under the time evolution and 
the quantum fields have to satisfy the microcausality conditions. Finally, it can be shown 
[495] how a perfect clock, in the sense of a quantum observable T whose values mono- 
tonically grow with abstract time t, is not compatible with the physical requirement of a 
positive spectrum of the energy, and this is a peculiar feature of the quantum theory. 

The problem of time arises also in the canonical formulation of the quantum theory 
of gravity, as happens in any diffeomorphism-invariant quantum field theory, and the 
Schrodinger equation is replaced by a Wheeler-DeWitt one, where the time coordinate is 
not present in the formalism. 

The proposals to address this fundamental problem are often related with the intro- 
duction of a reference system. In fact, this can be achieved in two different ways: the 
first one consists in adding a dynamical fluid or fields to the vacuum gravitational pic- 
ture (see Il99l 14331 14321 [44]), while the other in fixing the frame in a geometrical way 
(1 3161 l99l 13801 1362| ). However, both of them lead to an equivalent evolutionary quantum 
dynamics |433t I432T l99l 1361 1 . Below we will discuss the Brown and Kuchar mechanism, 
the so-called evolutionary quantum gravity and, finally, the multi-time approach. 

5.2.1 The Brown and Kuchar mechanism 

This approach |99| is devoted to find a medium leading to a Schrodinger equation when 
applying the Dirac quantization to a constrained system. In particular, an incoherent 
dust, i.e. one with the gravitational interaction only, is included in the dynamics. This 
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procedure leads to the new constraints H-j-(X) and H| a (X) in which the dust plays the 
role of time and the true Hamiltonian does not depend on the dust variables. 

Let us introduce the variables T, Z a and the corresponding conjugate momenta M, W a , 
so that the values of Z K be the comoving coordinates of the dust particles and T be the 
proper time along their worldlines. In this scheme the new constraints read as 

H T = P(X) + h(X, Kp, LT^) = (5.2.1a) 

H u = P«(X) + /7 a (X,T,z a A /3 ,n^) = (5.2.1b) 

where 

h = -tJg(X) , G(X) = H\ grav ) - h a ^H^ grav yH( gmv ^ , (5.2.2a) 

K = Zilip + ^G(X)3 /s Tzf , (5.2.2b) 

H(gmv) an d H(grav)a being the usual scalar and momentum constraints, respectively, P the 
projection of the rest mass current of the dust into the four velocity of the observers, and 
P a = —PW a . This way the Hamiltonian h does not depend on the dust. 

The quantization of this model is performed in the canonical way and yields, from 
equation (|5.2.1a|) , the Hamiltonian constraint operator 

H T = P(X) + h{X, Kp, tl'P) = . (5.2.3) 
This mechanism leads to a Schrodinger equation for the wave functional Y = Y(T,h) 

i— = hY. (5.2.4) 

The central point of this procedure is the independence of the effective Hamiltonian h(X) 
on the dust; this allows a well posed spectral analysis formulation because h commutes 
with itself. Furthermore, the Schrodinger equation can be split into a dust- (time-) depen- 
dent part and a truly gravitational one. 



5.2.2 Evolutionary Quantum Gravity 

In this section we remark some of the fundamental aspects of the evolutionary quan- 
tum gravity as presented in 113801 13621 ). As a first step we analyze the implication of a 
Schrodinger formulation of the quantum dynamics for the gravitational field|235], and 
then we establish a dualism between time evolution and matter fields. 

Let us assume that the quantum evolution of the gravitational field is governed by the 
smeared Schrodinger equation 

id t Y = 0iY = J d 3 x(NH)Y, (5.2.5) 

where the wave functional Y is defined on the superspace, i.e. it is annihilated by the 
super-momentum operator H a . Let us take the following expansion for the wave func- 
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tional 

rt 



Y = J De X {e,{h^})ex V -ijf dt' J^d 3 x{Ne) 



(5.2.6) 



where De is the Lebesgue measure on the space of the functions e(x). This expansion 
reduces the Schrodinger dynamics to an eigenvalue problem of the form 

H X = ex, H*X = 0, (5.2.7) 

which outlines the appearance of a non-zero super-Hamiltonian eigenvalue. 

In order to reconstruct the classical limit of the dynamical constraints \5.2.7\ . we replace 
the wave functional X by its corresponding zero-order WKB approximation x ~ e ' S • m 
this case the eigenvalue problem {5.2.7) reduces to its classical counterpart 

HJS = e = -iVhToo, HJ a S = (5.2.8) 

where H J and H f a denote operators which, when applied to the phase S, reproduce the 
super-Hamiltonian and super-momentum Hamilton-Jacobi terms, respectively. We see 
that the classical limit of the Schrodinger quantum dynamics is characterized by the ap- 
pearance of a new matter contribution (associated with the non-zero eigenvalue e) whose 
energy density reads as 

^ Too = -0r (5 - 2 - 9) 

where by Tu we refer to the new matter energy-momentum tensor. 

Since the spectrum of the super-Hamiltonian has, in general, a negative component, 
then we can infer that, when the gravitational field is in its ground state, this matter 
comes out in the classical limit with a positive energy density. The explicit form of (|5.2.9j) 
is that of a dust fluid co-moving with the slicing of the three-hypersurfaces, i.e. the normal 
field n' becomes the four-velocity of the appearing fluid (in other words, we deal with an 
energy-momentum tensor Tu = ptiiUj). 

We stress how the space of the solutions can be turned into the Hilbert one and therefore 
a notion of probability density naturally arises, from the squared modulus of the wave 
functional. 

Let us now consider the opposite sector, i.e. a gravitational system in the presence of a 
macroscopic matter source. In particular, we chose a perfect fluid with a generic equation 
of state p = (7 — l)p (p being the pressure and 7 the polytropic index). The energy- 
momentum tensor, associated to this system reads as 

Tij = ypUiUj - (7 - l)pgij ■ (5.2.10) 

To fix the constraints when matter is included in the dynamics, let us make use of the 
relations 



H 
'27h 



Gijrtni = -K :r -^= (5.2.11a) 



Gijti^yj = k^=, (5.2.11b) 
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where d K y' are the tangent vectors to the three-hypersurfaces, i.e. nid^y' = 0. Equations 
(|5.2.11a|) and (|5.2.11b|) , by (|5.2.10[) and identifying u\ with (i.e. the physical space is 
filled by the fluid), rewrite as 

H Hi = 0, (5.2.12) 



and furthermore, we get the equations 

Gijd u y l d0i = G a/3 = x(y - l)ph u p . (5.2.13) 

The conservation law V/lj = implies the additional two conditions 

7 Vi (pu ! ) = (7 - l)u%p (5.2.14a) 
u^'V/M; = ( 1 - ~\ (di lnp - m,m ; 3 ; - Lrijo) . (5.2.14b) 

With the space-time slicing, looking at the dynamics into the fluid frame (i.e. n' = S' ), 
by the relation n 1 = (1/N, —N a /N), the co-moving constraint implies the synchronous 
nature of the reference frame. Since a synchronous reference is also a geodesic one, the 
right-hand side of equation (|5.2.14b|) must identically vanish and, for a generic inhomo- 
geneous case, this implies 7 = 1. Hence, equations (|5.2.14a|) yields p = — e(x)/2yh and 
substituted into (|5.2.12|) , we get the same Hamiltonian constraints associated to the Evo- 
lutionary Quantum Gravity given above in eq. (|5.2.7| ), as soon as the function e is turned 
into the eigenvalue e. In this respect, while e is positive by definition, the corresponding 
eigenvalue can also take negative values because of the structure of H. 

Thus, we conclude that a dust fluid is a good choice to realize a clock in Quantum 
Gravity, because it induces a non-zero super-Hamiltonian eigenvalue into the dynamics; 
furthermore, for vanishing pressure (7 = 1), the equation (|5.2.13|) reduces to the proper 
vacuum evolution equation for h a a, thus outlining a real dualism between time evolution 
and the presence of a dust fluid. 

For the cosmological applications of the above approach in the isotropic sector see 
|154H381 | and in a generic cosmological model see [382 J and [44 J where it is shown how, 
from a phenomenological point of view, an evolutionary quantum cosmology overlaps 
the Wheeler-DeWitt approach. 

5.2.3 The Multi-Time Approach 

The multi-time approach 113141 [280] represents an alternative interesting way to get a 
Schrodinger quantum dynamics, even if deeply different from the one described above, 
since it is based on the ADM reduction of the dynamics. In fact, starting from equation 
J4A41 . 

P A {x)+h A {x, X ,(p,n) =0, (5.2.15) 

and performing a canonical quantization of the model, we obtain a Schrodinger-like equa- 
tion 

i T ^ = h A Y / (5.2.16) 
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where %a is the operatorial version of the classical Hamiltonian density. Even if this ap- 
proach and the evolutionary one (discussed in the previous subsection) seem to overlap 
each other when solving the problem of the frozen formalism, this is not the case. As mat- 
ter of fact, the evolutionary quantum dynamics approach is based on a full quantization 
of the system. By the other hand, the multi-time approach relies on a quantization of only 
some degrees of freedom. In fact, the constraints are classically solved before implement- 
ing the quantization procedure, thus violating the geometrical nature of the gravitational 
field in view of real physical degrees of freedom. This fundamental difference between 
the two approaches is evident, for example, in a cosmological context. In fact, when we 
quantize a minisuperspace model in the ADM formalism, the scale factor of the Universe 
is usually chosen as a "time" coordinate, and therefore the dynamics is consequently ex- 
pressed. By the other hand, in an evolutionary approach the scale factor is treated on the 
same footing of other variables, i.e. the anisotropics, and the evolution of the system is 
considered with respect to a privileged reference frame. 

5.3 The Minisuperspace Representation 

Although the full theory ("Quantum General Relativity") is far from being reached, many 
approaches are properly treated in the context of the so-called minisuperspace, as for ex- 
ample quantum cosmology. In fact, only a finite number of gravitational degrees of free- 
dom are invoked in the quantum theory, and the remaining ones are frozen out imposing 
some symmetries on the spatial metric. These space-times are, for instance, the homoge- 
neous cosmological models. 

In this sense quantum cosmology is a toy model for quantum gravity (with finite de- 
grees of freedom) which is a simple arena to test ideas and constructions introduced in 
the full theory (a genuine quantum field theory). In particular, since on a classical level 
the Universe dynamics is described by such symmetric models, their quantization is re- 
quired to answer the fundamental questions like the fate of the classical singularity, the 
inflationary expansion and the chaotic behavior of the Universe toward the singularity. 

Moreover, as we will see in Section 16.11 in the general context of inhomogeneous cos- 
mology, the spatial derivatives in the Ricci scalar are negligible with respect to the tem- 
poral ones, toward the singularity. This is the well-known BKL scenario |54l |57| where, 
as the singularity is approached (on a classical level), the spatial geometry can be viewed 
as a collection of small independent patches, in general Bianchi IX like models. There- 
fore a minisuperspace reduction of the dynamics is important also in the description of 
a generic Universe toward the classical singularity when restricted to each cosmological 
horizon. 

5.4 On the Scalar Field as a Relational Time 

Let us now discuss in some details the role of a matter field, in particular that of a 
scalar field <p, used as a definition of time for the quantum dynamics of the gravitational 
field. As we have seen, when the canonical quantization procedure is applied, the usual 
Schrodinger equation is replaced by a Wheeler-DeWitt one in which the time coordinate 
is dropped from the formalism. One possible solution to this fundamental problem is 
adding a matter field in the dynamics and then let evolve the physical degrees of freedom 
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of the gravitational system with respect to it. This way the dynamics is described on a 
relational point of view, i.e. the matter field behaves as a relational time. Such an idea is 
essentially based on the absence of time at a fundamental level and therefore a field can 



evolve with respect to another one only (for more details on this approach see [435]). 

In particular, in quantum cosmology, such choice appears as the most natural one. In 
fact, near the classical singularity, a monotonic behavior of a massless (p always appears 
as a function of the scale factoiQ (more precisely the variable which describes the isotropic 
expansion of the Universe). 

Let us consider the case of the Bianchi IX model in presence of a massless scalar field, 
whose Hamiltonian constraint in the Misner variables (a = e a , j5±) has the form 



2 -2 



4k kh ' a 



V(j8±) + ^«0, (5.4.1) 



V(B±) being the relative potential term. 

When this system is canonically quantized, the associated Wheeler-DeWitt equation de- 
scribes how the wave function Y = Y(a, B±,<p) evolves with (p. More precisely, from 
d5Ail we obtairl 



dl Y = SY, E = k 



« 4 



(5.4.2) 



behaving as a Klein-Gordon equation with (p playing the role of (relational) time and S 
of the spatial Laplacian. An explicit Hilbert space arises after performing the natural 
decomposition of the solution into positive and negative frequencies parts. In particular, 
the positive frequency sector offers the Schodinger-like equation — z3<pY = a/EY. 

We can analyze in which sense the scalar field can be regarded as a good time parameter 
for the dynamics. First of all, near the cosmological singularity (a — > 0) the potential 
term V(/3±) can be neglectecH. In such approximation the (classical) equations of motion 
obtained from (|5.4.1|) read as 

a = , ^ , Pa = ~ , ^ - P + = P-=0, (5.4.3) 
where (...) = d(...)/d<p and pi = p\ + The solutions of system (|5.4.3j) have the form 
.(♦) ^Bexp | -*Jj , Mt) .^--^-]. (5.4.4) 

A and B being integration constants and = const. 

1 In fact, such behavior well approximates the one of an inflaton field when its potential is negligible at 
enough high temperature] 140 [. For a detailed discussion of the inflationary scenario within the frame- 
work of homogeneous cosmologies, see 1431H5151l226ll266l IT32ll252ll449l |204l l62l l288H117ll353l . 

2 The choice of the normal ordering is not important for the following discussion, thus we adopt the simplest 
one. 

3 For a discussion on the consistency condition ensuring that the quasi-classical limit of the Universe dy- 
namics is reached before the potential term becomes important, see l44l 
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We have recovered a monotonic dependence of the scalar field (p with respect to the 
isotropic variable of the Universe a and therefore the massless (p shows to be a good 
(relational) time for the gravitational dynamics. 

As we will see in Section 16.31 the dynamics toward the cosmological singularity of a 
generic inhomogeneous Universe is described, point by point, by the one of a Bianchi IX 
model. More precisely, the spatial points dynamically decouple toward the singularity 
and the spatial geometry can be viewed as a collection of small patches, each one inde- 
pendently evolving as a Bianchi IX model. From this point of view, the monotonic relation 
(|5.4.4|) is a proper general feature of the gravitational field, better clarifying the choice of 
the scalar field as a relational time. 

5.5 Interpretation of the Universe Wave Function 

In this Section we face the problem of the probabilistic interpretation of the Universe wave 
function, in agreement with the analysis developed in |502l 1501 II (for a different point of 
view see \29S\ 1274, 223J). In fact, in quantum cosmology, the Universe is described by 
a single wave function Y providing puzzling interpretations when analyzing the differ- 
ences between ordinary quantum mechanics and quantum cosmology. 

In quantum mechanics, given a wave function Y^;, t ) describing a system, the proba- 
bility to find the system in a configuration-space element dClq at time t is given by 

dP = \W{qi,t)\ 2 dn q , (5.5.1) 

and it is positive semi-definite, i.e. dP > 0. On the other hand, the Universe wave function 
(which is a solution of the Wheeler-DeWitt equation) depends on the three-geometries, on 
the possible matter fields and no dependence on time appears. Therefore, the associated 
probability (here we denote whole the set of the superspace variables by h) 

dP=\Y{h)\ 2 VSd n h, 9 = det(S a/57<5 ) , (5.5.2) 

(S being the supermetric) is not normalizable, because its integral over the whole super- 
space is diverging. Such behavior can be considered as the analogue of the quantum 
mechanical feature 

J \Y(q if t)\dCi q dt = oo. (5.5.3) 

In fact, in quantum cosmology, the "time" is included among the set of variables h. 

An alternative definition of the Universe probability can be given in terms of con- 
serve^ current \V77\ (for notation see Vilenkin ||502|| ) 

Vf/ = 0, f = (Y*V ; Y - YV/Y*) . (5.5.4) 

This approach is limited as well since the corresponding probability to find the Universe 
in a surface element dS a is 

dP = fdSL u (5-5.5) 

4 This idea relies on the consideration that the WDW equation is nothing but a Klein-Gordon equation with 
variable mass. 
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and it can be negative, similarly to the problem of negative probabilities in the Klein- 
Gordon equation. 

In order to correctly define a probability, i.e. a positive semi-definite one, we can pro- 
ceed in two different ways: first we can consider a Universe, in which all the variables 
are treated semiclassically, then we can analyze a Universe where a quantum subsystem 
is taken into account, i.e. when some variables are pure quantum ones. 

In the case of semiclassical variables, the wave function Y is given by 

Y = A(h)e isw , (5.5.6) 

which admits a WKB expansion and leads to a conserved current f = | A| 2 V'S. The clas- 
sical action S(h) describes a congruence of classical trajectories and we shall define prob- 
ability distribution on the (ft — 1) -dimensional equal-time surfaces. Requiring only single 
crossing between the trajectories and these equal-time surfaces, the probability (|5.5.5[) re- 
sults to be positive semi-definite. 

Let us now consider the case in which not all the superspace variables are semiclassi- 
cal. We shall assume that the quantum variables are labeled by q and their effect on the 
dynamics of the semiclassical variables h can be neglected. Therefore, the WDW equa- 
tion can be decomposed in a semiclassical part Ho, which corresponds to the one of the 
previous case, and a "quantum" one H q . The fact that the quantum subsystem is small is 
shown by the existence of a parameter A (proportional to h) such that H^Y/ HoY = 0(A). 
Thus, since Ho = 0(A~ 2 ), therefore H q = 0(A _1 ). The superspace metric can then be 
expanded in terms of A as 

S//(fc/?) = S|(fc) + 0(A) / (5.5.7) 
and the Universe wave function can be written as 

Y = A(h)e iS ^x(Kq), (5.5.8) 

where the function % has to satisfy 

[V§ + 2[V (ln A)] V + 2i(V S) V -H q ] X = 0, (5.5.9) 

where the operator Vo is built using the metric Q°(h). The first two terms are of higher 
order in A with respect the third and can be neglected, resulting into 

2i(V S)VoA = H qX , (5.5.10) 

which can be written as a Schrodinger equation as idtX = NH q x (we recall that VoS 
coincides with the classical conjugate momentum to the semiclassical variable). Thus, we 
can get two different currents, one for the components in the classical subspace and one 
for those in the quantum one with the corresponding probability distribution written as 

p(h,q,t) = p (h,t) \x{qMt),t)\ 2 > (5-5.11) 

where po(h, t) and \x\ 2 are the probability distribution for the classical and the quantum 
variables, respectively. Considering the surface element on equal-time surfaces = 
cULodClq, for dZo defined from the metric 9°(h), the probability distribution (|5.5.11|) results 
to be normalizable. 
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By other words, we have recovered the standard interpretation of the wave function 
for a small subsystem of the Universe (only). This result agrees with the intrinsic ap- 
proximate interpretation of the Universe wave function. In fact, in the interpretation of 
quantum mechanics, all realistic measuring devices have some quantum uncertainty. In 
particular, the bigger is the apparatus the smaller are the quantum fluctuations. In this 
sense, we are able to give a physical interpretation of the wave function of the Universe 
only in a domain in which some variables are semiclassical. 

Finally, we recall two assumptions underlying this model: 

(z) the analysis has been developed in the minisuperspace homogeneous models only; 

(ii) the fundamental requirement of existence of a family of equal-time surfaces is taken 
as a general feature. 



5.6 Quantization in the Misner Picture 

In this Section we provide a quantum representation of the dynamics, relying on the adi- 
abatic approximation ensured by the potential term, reduced to an infinite well and, ac- 
cording to C. W. Misner |368j|, we model the potential as an infinite square box with the 
same measure as in the original triangular picture. The volume-dependence of the wave 
function acquires increasing amplitude and frequency of oscillations as the Big Bang is 
approached and the occupation number grows, respectively. 

For a review of the canonical quantization of Bianchi cosmologies in the WDW frame- 
work see ll^mi4l^lTT2l l7^ TOl40^ 

14561 12531 12021 12011 |27l EEU [1351 HE [353 EH- For a discussion on the quantization in a 
supersymmetric form of the Bianchi IX model, see fl218H22ifl . 

In a canonical framework, by replacing the canonical variables with the correspond- 
ing operators and implementing the Hamiltonian constraints we get the state function 
describing the system xp = ip{a, i°> + , 

Adopting the canonical representation in the configuration space we address the WDW 
equation 

5 2 5 2 S 2 ' 
+ 



Hip 



5oc 2 5fS% Sfi 2 



xp-e R Vxp = 0. 



(5.6.1) 



We can find a solution in the form 



tp = Y^T n (a.)(p n (<x,, p +r ft-) , 



(5.6.2) 



where the coefficients T n are a -dependent amplitudes. 

Thus, equation 05.6.11 ) is reduced to the ADM eigenvalue problem 



a 2 9 2 



e ia V 



cp„ = E 2 {ot)cp n . 



(5.6.3) 



According to M368L we approximate the triangular infinite walls of the potential by a box 
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having the same measure to find the eigenvalues E n 

( 4 \ 1/2 |n| 
£„(*) = tt^J LA, (5.6.4) 

where n 2 = n 2 + + and n + , n_ G N are the two independent quantum numbers corre- 
sponding to the variables /3+, /3_, respectively. 

Substituting the expression for t/> in equation (I5.6.1I I we get the differential equation for 

T„ 

Y2{dlT n )cp„ + ^n{dl(p n ) + 2j^{d ct T n ){d a( p n ) + J2E 2 n T n cp„ = , (5.6.5) 

n n n n 

which, in the limit of the Misner adiabatic approximation of neglecting d a (p n (i.e. (p ~ 
<p(P+, simplifies to 

-j^ + V* = (5.6.6) 

where 



(5.6.7) 



The above equation is solved by r„(«) in the form 

r n (a) = Cii/asin ^y^lna^ +C2v / *cos ^y^lna^ , (5.6.8) 



where i/p^ = -y/A;^ — 1. From (|5.6.8[) the self-consistence of the adiabatic approximation 
is ensured. Figure IBTTI shows the behavior of T n (ot) for various values of the parameter k n . 
Such wave function behaves like an oscillating profile whose frequency increases with 
occupation number n and approaching the Big Bang, while the amplitude depends on 
the a variable only. 

By this treatment, one finds [368 J the interesting result that n on average is constant 
toward the singularity and then if the initial state of evolution of the Universe is classical, 
extrapolating backwards it maintains a semiclassical character. 



5.7 The quantum Universe in the Poincare half-plane 

The Misner representation provides a good insight in some qualitative aspects of the Mix- 
master model quantum dynamics, and allows some physical considerations on the evo- 
lution toward the singularity. In this picture the potential walls move with time, and this 
is an obstacle toward a full implementation of a Schrodinger like quantization scheme. 
These difficulties can be by-passed as soon as MCI variables are adopted, characterized 
by static potential walls; in particular, we will choose the so-called Poincare variables 
(m, v), defined as 

1 + U + U 2 + V 2 

£ = 7= , (5.7.1a) 

V3v 

v / 3(l + 2u) 

9 = - arctan - \ v — - = ' , . (5.7.1b) 
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In the vicinity of the initial singularity, we have seen that the potential term behaves as 
a potential well and as soon as we restric the dynamics to Uq, "Kadm = e and we can 
rewrite (|4.8.15|) and <(4.8.16|> as 

5Sn Q = 5 J dT( n t + p e 6 - H ADM ) = (5.7.2a) 



H ADM = Vy/pl + pi. (5.7.2b) 

The asymptotic dynamics is defined in a portion YIq of the Lobatchevsky plane, delimited 
by inequalities 

Qi(u,v) = —u/d > (5.7.3a) 
Q 2 (u,v) = (l + u)/d > (5.7.3b) 
Q 3 (u,v) = (u(u + l) +v 2 )/d > (5.7.3c) 
d = l + u + u 2 + v 2 , (5.7.3d) 

whose boundaries are composed by geodesies of the plane, i.e. two vertical lines and one 
semicircle centred on the absolute v = 0. 

The billiard has a finite measure, and its open region at infinity together with the two 
points on the absolute (0,0) and ( — 1,0) correspond to the three cuspids of the potential 
in Fig. [EH 

It is easy to show that, in the u, v plane , l493l becomes 

1 dudv dtp 

du = n-TT-. (5.7 A) 
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Figure 5.2: Tq(u,v) is the available portion of the configuration space in the Poincare 
upper half-plane. It is bounded by three geodesic u = 0,u = — 1, and (u + 1/2) 2 + v 1 = 
1/4, it has a finite measure p. = n (it is "dynamically" bounded). 



5.8 Continuity Equation and the Liouville theorem 

Since the Mixmaster provides an energy-like constant of motion toward the singularity 
the-point Universe randomizes within a closed domain and we can characterize the dy- 
namics as a microcanonical ensemble, as discussed in Section l4~8l 

The physical properties of a stationary ensemble are described by a distribution func- 
tion p = p(u, v, p u , p-o) |268], representing the probability of finding the system within an 
infinitesimal interval of the phase-space (u, v, p u , p v ), and it obeys the continuity equation 

d{up) d(vp) d(p u p) d{p v p) 

1 5 1 5 1 5 = U, (5.8.1) 

were the dot denotes the time derivative and the Hamilton equations associated to d5.7.2bl) 
read as 

v 1 

u = — p u , pu = , (5.8.2a) 

2 

v e 
v = —p v , p v = -- . (5.8.2b) 

e v 

From (I5.8.1I) and (15.8.21) we obtain 

v 2 p u dp v 2 p v dp e dp 



du e dv v dp v 



. (5.8.3) 



The continuity equation provides an appropriate representation when we are sufficiently 
close to the initial singularity only, and the infinite-potential-wall approximation works. 
Such model corresponds to deal with the energy-like constant of motion, and fixes the 
microcanonical nature of the ensemble. From a dynamical point of view, this picture 
naturally arises because the Universe volume element monotonically vanishes (for non 
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stationary corrections to this scheme in the MCI variables, see [379J). 
Since we are interested to the distribution function in the (u,v) space, we will reduce 
the dependence on the momenta by integrating p(u,v,p u ,p v ) in the momentum space. 
Assuming p to be a regular, vanishing at infinity in the phase-space, limited function, we 
can integrate over (|5.8.3|) getting the equation for w = w{u, v; k) 

2C 2 v 2 w 

0, (5.8.4) 




V£ 2 " {Cvf 

where the constant C appears, due to the analytic expression of the HJ solution, fixed 
by the initial conditions. However, we deal with a distribution function that cannot con- 
sider these initial conditions, and must be ruled out from the final result. We obtain the 
following solution in terms of the generic function g 



w(u, v; C) = — ^ — 7 - . (5.8.5) 

dVE — C z v z 

The distribution function cannot contain the constant C, and the final result is obtained 
after the integration over it. We define the reduced distribution w{u, v) as 

w(u,v) = / w(u,v;k)dk , (5.8.6) 

J A 

where the integration is taken over the classical available domain for p u = C 

A = 



E E 

v' v 



(5.8.7) 



In (|5.7.4|) we proved demonstrated that the measure associated to it is the Liouville one; 
the measure w mc (after integration over the admissible values of cf>) corresponds to the 
case g = const. 

/§ 1 ji 

\ dC = -j . (5.8.8) 



Summarizing, we have derived the generic expression of the distribution function fixing 
its form for the microcanonical ensemble. This choice, in view of the energy-like constant 
of motion 'Kadm, is appropriate to describe the Mixmaster system restricted to the con- 
figuration space. This analysis reproduces in the Poincare half-plane the same result as 
the stationary invariant measure described in Section l4!9l 



For completeness we report the explicit solution to (|2.4.2H) fo this model in the restricted 
domain TIq as the Hamilton-Jacobi function for the point-Universe 



% {u,v) =Cu + y / e 2 -C 2 v 2 -eln 2 V ^ +D, (5.8.9) 

where C is the separation constant, and D is an integration one. 
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5.9 Schrodinger dynamics 

The Schrodinger quantum picture is obtained in the standard way, i.e. by promoting 
the classical variables to operators and imposing some boundary condition to the wave 
function. 

The first step is achieved as 

0|) =v |) , u |) =u |) , 

fc-*-'^' (5 - 9 - 1} 
while for the second one we will require the Dirichlet boundary conditions 

0(an Q ) = . (5.9.2) 

The quantum dynamics for the state function <I> = <E>(u, v, t) obeys the Schrodinger equa- 
tion 



Here we have adopted a generic operator-ordering for the position and momentum parametrized 
by the constant a|310], as soon as we have no indication on it providing a first problem. 
The other one is linked to the multi-time approach, i.e. to the non-locality of the Hamilto- 
nian operator: when solving the Hamilton constraint with respect to one of the momenta, 
the ADM Hamiltonian contains a square root and consequently it might define a non-local 
dynamics. 

The question of the correct operator-ordering is addressed the next Section comparing 
the classic evolution versus the WKB limit of the quantum-dynamics and requiring the 
overlapping of the two. On the other hand, we will assume the operators 'Kadm an d 
•^"ADM having the same set of eigenfunctions with eigenvalues E and E 2 , respectively]! 

Under these assumptions, we will solve the eigenvalue problem for the squared ADM 
Hamiltonian given by 



W 2 Y 

where Y = Y(u,v,E). 



du 2 dv \ dv 



Y = E 2 Y , (5.9.4) 



5.10 Semiclassical WKB limit 

In order to study the WKB limit of equation (|5.9.4|) , we separate the wave function into its 
phase and modulus^] 

Y(u, v, E) = ^r(u r v,E)e'^ u ' v ' E)/h . (5.10.1) 

5 The problems discussed by |296| do not arise here because in the domain YIq the ADM Hamiltonian has a 

positive sign (the potential vanishes asymptotically). 
6 In this Section we restore h in the notation because we deal with the semiclassical limit. 
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In Eq. (|5.10.1|) the function r(u,v) represents the probability density, and the quasi- 
classical regime appears in the limit h — *■ 0; substituting (|5.10. 1|> in (|5.9.4|) and retaining 
only the lowest order in h, we obtain the system 



v 



da 



du 

dr da dr da 
du du dv dv 



E 2 , 

a da d 2 a d 2 a 



' v dv ^ dv 1 



du 1 



. 



(5.10.2a) 
(5.10.2b) 



In view of the HJ equation and of Hamiltonian (|5.7.2b|) , we can identify the phase a as the 
functional So- 

Taking (|5.8.9[) into account, Eq. (|5.10.2b|) reduces to 




. 



(5.10.3) 



Comparing (|5.10.3|) with (|5.8.4|) , we see that they coincide for a = 2 only [61 J. This corre- 
spondence is expectable for a suitable choice of the configurational variables; however, it 
is remarkable that it arises for the chosen operator-ordering only. Here arises the impor- 
tance of the correspondence, wich fixes a particular quantum dynamics for the system. 
Summarizing, we have demonstrated from our study that it is possible to get a WKB cor- 
respondence between the quasi-classical regime and the ensemble dynamics in the con- 
figuration space, and we provided the operator-ordering when quantizing the Mixmaster 
model 

' X (5.10.4) 



~2 ~2 



dv 



dv 



5.11 The Spectrum of the Mixmaster 



5.11.1 Eigenfunctions and the vacuum state 



Once fixed the operator ordering a = 2, the eigenvalue equation (I5.9.4I) rewrites as 

}2 



, 3 2 2 d 2 „ 9 



Y(u,v,E) = 



(5.11.1) 



By redefining Y(m, v, E) = ip(u,v,E) / v, we can reduce (|5.11.1|) to the eigenvalue problem 
for the Laplace-Beltrami operator in the Poincare plane [484 j 

/ d 2 d 2 \ 

\7 LB xp(u, v ,E) = v 2 l^ 1 + ^\ tp(u, v, E) = E s xp(u, v, E) , (5.11.2) 



which is central in the harmonic analysis on symmetric spaces and has been widely in- 
vestigated in terms of its invariance under SL(2,C) and its eigenstates and eigenvalues 
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are 

xp s (u,v) = av s + bv l - s + y/v a n K s _ 1/1 {2n\n\v)e 2ninu , a,b,a n £ C , (5.11.3) 

\7 LB ip s (u,v) = s(s - l)f s (u,v) , (5.11.4) 

where K s _i/ 2 (2nnv) are the modified Bessel functions of the third kind|l] and s denotes 
the index of the eigenfunction. This is a continuous spectrum and the sum runs over 
every real value of n. 

The eigenfunctions for our model, then, read as 

YCu,v,E) = av s - 1 + bv- s + Y an K s-V^\ n \ v \ ^inu ^ (5 n 5) 

with eigenvalues 

£ 2 = s(l-s). (5.11.6) 

To impose Dirichlet boundary conditions for the wave functions we will require a van- 
ishing behaviour on the edges of the geodesic triangle of Fig. 15.21 Let us approximate|61] 
the domain with the one in Fig. I5.3t the value of the horizontal line v = 1 / n provides the 
same measure for the exact as well as for the approximate domain 

dudv f dudv 

-=3T= / — = n - (5 - 1L7) 

llQ ^ J Approx domain u 



2 5 - 



J - 



> 1.5 - 




Figure 5.3: The approximate domain where we impose the boundary conditions. The 
choice v = I /re for the straight line preserves the measure }i = n (from 161 J). 



The difficulty to deal with the exact boundary conditions relies in the sophisticated 
number theory is linked to these functions while the circle, that bounds from below the 
domain, furthermore mixes solutions with different indices s. 

The Laplace-Beltrami operator and the exact boundary conditions are invariant under 
parity transformation u — > — u — 1; however, the full symmetry group is Cy , and this 
can be seen in the disk representation of the Lobachewsky plane. C$ v has two one- 
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dimensional irreducible representations and one two-dimensional representation[158, 220J. 

The eigenstates transforming according to one of the two-dimensional representations 
are twofold degenerate, while the others are non-degenerate. The latters can be divided in 
two classes, either satisfying Neumann boundary conditions, or Dirichlet ones. We focus 
our attention to second case. The choice of the line v = 1/rc approximates symmetry lines 
of the original billiard and corresponds to one-dimensional irreducible representations|419 , 
Il58ll220ll . 

The conditions on the vertical lines u = 0, u = — 1 require to disregard the first two terms 
in (|5.11.5|) ; furthermore, we get the condition on the last term 

oo 

£ gtomu ^ £ sin(TrnM) , (5.11.8) 

n^O n=\ 

for integer n. As soon as we restrict to only one of the two one-dimensional representa- 
tions, we get 

oo 

£ e 2mnu _^ £ S in(27rnw) , (5.11.9) 

n=l 

while the condition on the horizontal line implies 

E a n K s - 1/2 (2n) sm(2nnu) = , Vw G [-1, 0] , (5.11.10) 

)!>0 

which in general is satisfied by requiring K s _i/ 2 (2ft) = only, for every n. This last 
condition, together with the form of the spectrum (|5.11.6[) , ensure the discreteness of the 
energy levels, thanks to discreteness of the zeros of the Bessel functions. 
The functions K v (x) are real and positive for real argument and real index, therefore the 
index must be imaginary, i.e. s = + it. In this case, these functions have (only) real 
zeros, and the corresponding eigenvalues turn out to be real and positive. 




Figure 5.4: The intersections between the straight lines and the curves represents the roots 
of the equation K# (n) = 0, where K is the modified Bessel function (from [61 1).. 
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E 2 = t 2 + -. (5.11.11) 

The eigenf unctions (|5.11.5|) exponentially vanish as infinite values of V are approached. 
The conditions (|5.11.10|) cannot be analytically solved for all the values of n and t, and 
the roots must be numerically worked out for each n. There are several results on their 
distribution that allow one to find at least the first levels: a theorem [404| on the zeros 
of these functions states that K lv (vx) = 0o0<x<l; furthermore, the energy levels 
(|5.11.11|) monotonically depend on the values of the zeros. Thus, one can search the lowest 
levels by solving Eq. (|5.11.10|) for the firsts n; in the next Section we will discuss some 
properties of the spectrum, while now we will discuss the ground state only. 

A minimum energy exists, as follows from the quadratic structure of the spectrum and 
from the properties of the Bessel zeros, and its value is Eq = 19.8317z 2 , and correspond- 
ingly the eigenfunction is plotted in Fig. 15.51 together with the probability distribution in 
Fig. I5.6t the eigenstate is normalized through the normalization constant N = 739.466. 
The existence of such a ground state has been numerically derived, but it can be inferred 
on the basis of general considerations about the Hamiltonian structure; the Hamiltonian, 
indeed, contains a term v 2 pl, which has positive definite spectrum and does not admit 
vanishing eigenvalues. 




Figure 5.5: The ground state wave func- Figure 5.6: The probability distribution 
tion of the Mixmaster model is sketched in correspondence to the ground state 
(from HSU). °^ * ne theory, according to the boundary 

conditions (from |6TH ). 



5.11.2 Distribution of the Energy Levels 

In Table IBTTI we report the first ten "energy" levels evaluated solving (|5.11.10|) F1 

To study the distribution of the highest energy levels, we need to consider the asymp- 
totic behavior of the zeros for the modified Bessel functions of the third kind. We will 
discuss the asymptotic regions of the (t,n) plane in the two cases t~^> n and t ~ n ^> 1. 



7 For a detailed numerical investigation of the energy spectrum of the standard Laplace-Beltrami operator, 
especially with respect to the high-energy levels, see|158. 220 1, where the effects on the level spacing of 
deforming the circular boundary condition towards the straight line are numerically analyzed. 
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Table 5.1: The first ten energy eigenvalues, ordered from the lowest one. 

EW + i 

19.831 

40.357 

49.474 

63.405 

87.729 

89.250 
116.329 
128.234 
138.739 
146.080 



(z) For t 3> n, the Bessel functions admit the representation 



K it (n) 



2ne- tn/2 

(f2_„2)l/4 



ic=0 : 



COS 



fc=0 



(-1)' 



U2k+1 



(5.11.12) 



where a = 7r /4 — Vt 2 — n 2 + facosh(f /n), p = n/t and are the polynomials 
w (f) = 1 , 

1 1 (3 11 13) 

»*+i(0 = ^(l - ^KO + s i, (1 " 5j2 > k ^ dt ■ 

Retaining in the expression above only terms of order o( n /t), the zeros are fixed by 
the relation 



sin 



71 



f + f(log(2)-log(p)) 



cos 



lit L4 



7T 



In the limit n/t <C 1, Eq. (|5.11.14[) can be recast as 
flog(f/n) = ln ^ t 



f + f (ln(2) -ln(p)) 



Ztt 



productlog ( 



= . 
(5.11.14) 

(5.11.15) 



where productlog(z) is a generalized function giving the solution of the equation 
z = we w and, for real and positive domain, is a monotonic function. In d5.il. 151) I is 
an integer number much greater than 1 in order to verify n/f < 1. 



(ii) In case the difference between In and t is o(n 1/3 ) for f, n ^> 1, we can evaluate the 
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first zeros k S/V by the relations [29 



~ /1A -(2r-l)/3 

fcs,v~v + J](-l) r s r (fl s ) (-) , (5.11.16) 



r=0 



where a s is the s-th zero of Ai((2/z) 1/3 ), Ai(x) is the Airy function lrTBl and s, are 
some polynomials. From this expansion it results that, to lowest order 

t = In + 0.030n 1/3 . (5.11.17) 

Eq. (|5.11.17|) provides the lowest zero (and therefore the energy) for a fixed value of 
n and also the relation for the eigenvalues for high occupation numbers 

E 2 ~ An 2 + 0.12n 4/3 . (5.11.18) 

Let us to discuss i.e. the completeness of the spectrum and the definition of a scalar 
product. 

The problem of completeness can be faced by studying firstly the sine functions and then 
the Bessel ones. On the interval [—1,0], the set sm(27inu) is not a complete basis, but as 
soon as we request the wave function to satisfy the symmetry of the problem, it becomes 
complete. 

Let us take a value n > 0, thus functions (|5.11.5[) have the form <E>(m, v) = sm(2nnu)g(v), 
which substituted in d5.ll.lD provides v 2 (d 2 + (27rn) 2 )) g(v) = s(l — s)g(v), whose solu- 
tions are exactly the Bessel functions. |419J 

This property togheter with the condition on the line v = 1 / n form a Sturm-Liouville 
problem and we deal with a complete of eigenf unctions. 

Therefore, such eigenfunctions define a space of functions where we can introduce a 
scalar product, naturally induced by the metric of the Poincare plane|419J 



= f ip( X ,y)(p*( X ,y)^L , (5.11.19) 



where * denotes complex conjugation. 

Now we briefly discuss if the presence of a non-local function, like the square-root 
of a differential operator, can give rise to non-local phenomena. Following the work of 
Puzio|419], a wavepacket which is non-zero in a finite region of the domain (v < M) and 
far from infinity fails to run to infinity in a finite time, i.e., the probability P(v > M) to 
find the packet far away exponentially vanishes. 



/0 poo / 
/ \vJdl+d 2 v W(u,v 
-1JM v 

4M 2 y|(supY) 2 |_° i 



.r,dvdu 



dv du 

< 4MM- (sup T P / / e 



M V 2 



= 4M 2 y|(supY) 2 (lJ^ + Ei(-2M)) < 

< W ^ (sup Y) 2 Me- 2M (5.11.20) 
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where sup(Y) is the maximum value of the wavepacket in the domain v < M and Ei(z) = 
— f e~ I tdt is the exponential integral function. We can thus conclude that nevertheless 
the square root is a non-local function, non-local phenomena don't appear (like the case 
of a wavepacket starting from a localized zone and falling out to infinity). 



5.12 Basic Elements of Loop Quantum Gravity 

In this Section we will review some of the basic aspects of Loop Quantum Gravity (LQG) 
in a pedagogical manner for a non-expert reader. For a detailed and deeper approach, we 
refer to the textbooks I4~35ll485l and to some dedicated reviews 14341 l4§6l l22l 14671 . while 
for a critical point of view we recommend [393 J. 

LQG is an attempt to rigorously quantize General Relativity in a background indepen- 
dent manner, trying to define a quantum field theory just on a differential manifold M 
and not on a background space-time (M,go), i.e. independently of the choice of a fixed 
background metric go. However, a way to give a formulation of the theory independently 
of the possible topological changes of the underling manifold has not been achieved so 
far in LQG. 

In order to see how fundamental fixed background metric is, let us consider ordi- 
nary quantum field theory, whose whole framework, namely the Wightman axioms[234J, 
breaks down as soon as the metric is no longer considered fixed but with a dynamical 
structure according to General Relativity. In fact we could only construct a rigorous quan- 
tum theory in the Minkowski background^, which implies a preferred notion of causality 
(locality) and a symmetry group, i.e. the Poincare one. 

LQG works in a Hamiltonian approach^ and it is able to overcome some of the problems 
of the older geometro-dynamical approach (described in the previous Sections) using new 
conceptual and technical ingredients. 

Let us rewrite GR in the form of a SU(2) Yang Mills theory, so that the phase space is 
endowed with coordinates q by a SU(2) connection and by an electric field EJ. These 
are related to the ADM variables (h, K) by 

4 = oi + 7*«/A Ef = Vhef, (5.12.1) 

where to{ = l /2e 0l p 1 co 1 ^ is the spin connection, K K p the extrinsic curvature and 7 > the 
so-called Immirzi parameter which does not affect the classical dynamics but brings to 
inequivalent quantum predictions [436 J. The symplectic geometry is determined by the 
only non-trivial Poisson brackets 

{Ai(x),E?(x')} = KjSjS&ix - x') . (5.12.2) 

Let us rewrite the constraints (|2.4.14|) in terms of these new variables, which in the phase 



8 The only quantum fields in four dimensions fully understood to-date are the free, or perturbatively inter- 
acting, fields. 

9 For the covariant approach to LQG, i.e. the spin foams models, see |400 408 1. 

"Here and in the following paragraphs on the LQG/C we will adopt i,j,k, ... = 1,2,3 as SU(2) indices. 
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phase satisfy the usual spatial diffeomorphism constraint 

H a = F^Ef « 0, (5.12.3a) 
where F = dA + A A A is the curvature of A, and the scalar one 



H 



\/\fet(E)\ 



F l ccp-{l + r)eimnKty »0, (5.12.3b) 



where 7K4 = A„ — o;^. In the connection formalism, with respect to the metric approach, 
we have the additional Gauss constraint 

9i = VaE? = d K Ef + e ijk A{El « (5.12.4) 

which gets rid of the SU(2) degrees of freedom. Since an observable is a gauge invariant 
function, it has to commute with all constraints. 



As a second step we have to smear |435H485ll the fields A\ and E°- to overcome some of 
the problems arising in the direct approach. In fact, we switch from the connections A l a 
to the holonomies as the basic variables. Given a curve on E, i.e. an edge i, a holonomy is 
defined as 

h £ [A] = Texp (J^ A^j = IPexp Q AfadxA , (5.12.5) 

where T denotes the path order and the t, = crjli form a basis of SU(2) and (T{ are the 
Pauli matrices. The holonomies hi are elements of SU(2) and define the parallel trans- 
port of the connection A\ along the edge £. They are gauge invariant and have a one- 
dimensional support (rather than all E). The variable conjugate to the hi\A\ is the flux 
vector 

P l s [E] = J s *E ! ' = J e ah E itt dxP A dx^ (5.12.6) 

through any two-dimensional surface S C E, whose support is a two-dimensional sub- 
manifold of E. In order to compute the Poisson brackets between these variables, let us 
consider an edge £ that intersects a surface S in one point, thus obtaining 

{h,[A\,P' s [E]} = K -la{e,S)T%[A], (5.12.7) 

where a = in the case of the edge not intersecting the surface and a = ±1 when the 
orientation of the edge and surface are the same or the opposite, respectively. We note 
that the commutation relation d5.12.7l ) is non-canonical for the presence of hi [A] . 

The quantum kinematics can be constructed promoting such variables to quantum op- 
erators obeying appropriate commutation relations. The essential feature of LQG is pro- 
moting the holonomies hi [A] to operators rather than the connections A' a themselves. 

Let us investigate the kinematical Hilbert space "Kun of LQG, i.e. that of the spin net- 
works. These are defined as a graph T consisting of a finite number of edges and vertices 
with a given collection of spin quantum numbers ji = 1, 3 /2, . . ., one for each edge, and 
of other quantum numbers I, the intertwiners, one for each vertex. The wave function of 
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the spin network, the so-called cylindric function, can thus be written as 

Y T/lj) [A]=xp(h £l [A],...,h {n [A}). (5.12.8) 

If the wave functions ip are SlT(2)-gauge invariant, they satisfy the Gauss constraint and 
viceversa. These functions are called cylindrical because they have a one-dimensional sup- 
port, i.e. they probe the gauge connection on one-dimensional networks only. When pro- 
moted to quantum operators, the fundamental variables of the theory (/^[/L], P S [E]) act 
on the wave functions (|5.12.8|) as 

ht[A]Y Tif (A) = h e [A]Y Tif (A) (5.12.9a) 
Ph[E]V Ti ip(A) =i{Pm,V Ti ip(A)} . (5.12.9b) 

A key point in the LQG approach to the quantum gravity problem is the kinematical 
scalar product between two cylindric functions, since the main results of discretization of 
areas and volumes are based on it, which is defined as 

(Y r |Yp) = {° *<l*r' {5 - 12A0) 

where the integrals J kg are performed with the SU(2) Haar measure|470]. The definition 
is based on a strong uniqueness theorem [ 324 1 . The inner product vanishes if the graphs T 
and r' do not coincide and it is invariant under spatial diffeomorphisms, even if the states 
Yi and Y2 themselves are not, because the information that the two graphs coincide is 
diffeomorphism invariant. The information on the position of the graphs, carried by the 
wave function, disappears in the scalar product (|5.12.10|) . 
Let us stress three relevant aspects: 

(z) in contrast with the lattice gauge theory, where all quantities depend on the scale 
parameter, the "discretuum" of LQG is built via the construction of the scalar prod- 
uct (|5.12.10|) . This feature is responsible for the failure of the Stone-von Neumann 
theorem [234J: LQG and WDW lead to different results, since the two quantizations 
procedures are not equivalent. 

(z'z) The obtained Hilbert space is not-separable, as it does not admit a countable basis, 
because the set of all spin networks is not numerable and two non-coincident spin 
networks are orthogonal with respect to (|5.12.10|) . 

(z'z'z) States with negative norm are absent without imposing the constraints (|5.12.3a|) , 
(|5.12.3b|) , (|5.12.4[) , in contrast with the usual gauge theories, where the negative 
norm states can be eliminated only after imposing the constraints. However, this 
kinematical Hilbert space is not relevant to solve the quantum constraints of LQG. 
Anyway, the area and volume operators (which we do not discuss here) are com- 
puted at this level and the discretization of the corresponding spectrum is related to 
this Hilbert space properties^ 



11 The spectrum of the area operator Ag is computed applying it to a wave function Y defined in a given 
graph. Such a graph is refined in a way that the elementary surface S/ is pierced only by the edge of the 
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The last step is to impose the constraint at a quantum level to obtain the physical states. 
These constraints are implemented firstly to express them in terms of holonomies and 
fluxes, and secondly to investigate their properties. 

The Gauss constraint (|5.12.4|) imposes several restrictions on the quantum numbers I of 
the cylindric functions. The diffeomorphism constraint d5.12.3al ) is more difficult to deal 
with and will not be treated in terms of an operatorial one. In fact, a diffeomorphism 
generator does not exist as an operator and diffeomorphism invariant states (with the 
exception of the empty spin network state Y = 1) do not exist in H-un- This constraint 
is imposed implementing a "group average method" in a way properly adapted to the 
scalar product (|5.12.10[) I2141. A key step arises from solving the constraints on a larger 
space, i.e. the dual Cyl* of the cylindric functions Cyl, such that 

Cyl c "Kun C Cyl* , (5.12.12) 

which is the so-called Gel'fand triple [234J . An important feature of the new Hilbert space 
'Kdiff is of being separable, differently from t K\~ in . The final challenge will be finding 
a space annihilated by all the constraints and defining a physical inner product which 
yields the final physical Hilbert space. 

As mentioned in the WDW formalism, the main problem of all quantum gravity theo- 
ries is to impose at a quantum level the scalar constraint. This difficulty appears also in 
the LQG approach and will be reflected also in the minisuperspace theory, i.e. the Loop 
Quantum Cosmology (LQC). 

Let us rewrite the constraint J5.12.3b|) in terms of variables corresponding to a well 
defined quantum operator. This can be done using classical identities from which we can 
express the triads, the extrinsic curvature and the field strength in terms of holonomies 
and well-defined operators as the volume one. Then, similarly to quantum field theory, 
the Hamiltonian has to be regularized with a parameter e shrunk to zero at the end of 
the computation. This is a highly non trivial step and poses serious challenges. In fact, 
the regularisation e enters in the spin networks picture via a plaquette P(e), attached to a 
vertex and thus it modifies the underlying graph, but two wave functions supported on 
different networks are orthogonal by the scalar product (|5.12.10|l . In other words, for any 
cylindric function Y, the limit e — > of H e Y does not exist on Cyl. Usually one transfers 
the action of the scalar constraint to the dual space, adopting a weaker notion of limit. 
More specifically, the limit e — > is defined by 

{H*xW = lim(x|H e Y) (5.12.13) 

for all Y G Cyl and ^ G ^* C Cyl*. The space V* must be selected taking into account 
physical motivations. A natural choice is V* = l Kdiff> but this is not unique and further- 
more poses several problems |393j. 

The necessary use of the weaker limit (|5.12.13|) is reflected on the quantum constraints 

network, therefore obtaining 

A s Y = K 7 £jj p (j p + l)Y. (5.12.11) 

v 

The minimal accessible length appears to be of order of the Planck one and the spectrum discreteness is 
related to the structure of the spin networks. In fact, a further refinement of the area operator does not 
contribute to the result. 
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algebra. In fact, in LQG a weaker notion of algebra closure (with respect to the so-called 
off-shell closure) is formulated, requiring that equation 

[H*,H'*]x = (5.12.14) 

holds. The closure of the algebra is required only after the imposition of the constraints, 
and in this sense it is weaker. This is relevant when addressing the quantum space-time 
covariance, because only the off-shell closure exhibits this quantum gravity property and, 
relaxing the algebra closure, some crucial information may get lost. 

Let us conclude noting that all the main result of LQG and LQC are related to the 
discretization of the spectrum of area and volume operators. The nature of this cut-off 
could be ambiguous|5, 455J, since it can be dynamically generated or imposed by hand 
when dealing with the compact SU(2) rather than the true gauge group of GR, i.e. the 
Lorentz one. In fact, in a covariant formalism, without imposing a time-gauge, the area 
spectrum is continuous [6], but for a discussion on this point of view we demand to 
|328] and references therein. For the recent interesting approaches to solve the scalar 
constraints problems mentioned above, see |21 1 H2121 12131 14871 1240fl . 



5.13 Isotropic Loop Quantum Cosmology 

The Loop Quantum Cosmology is a minisuperspace model, i.e. a truncation of the phase 
space of classical GR, which is quantized according to the methods of LQG. Therefore 
LQC is not the cosmological sector of LQG, i.e. the inhomogeneous fluctuations are 
switched off by hand rather than being suppressed quantum mechanically. However 
the minisuperspace is an important arena to test a theory and in LQC the most spectacu- 
lar results appear: the absence of the classical singularity |88| , replaced by a Big-Bounce 
|23l [24] in the isotropic settings, a geometrical inflation [89] and the suppression of the 
Mixmaster chaotic behavior toward the singularity [93]. 

In this Section we will analyze some basic aspects of the so-called isotropic Loop Quan- 
tum Cosmology and the fate of the classical singularity. For exhaustive discussions on 
this argument see f9"0ll2"T| . 

In this case, the phase space of GR is two-dimensional, since the scalar factor a = a(f) 
is the only degree of freedom of an isotropic model. Therefore, the variables (|5.12.1|) are 

reduced to 

A = c "coWt, E = pV°h °ei t 1 (5.13.1) 

where a fiducial metric in E is fixed by h a y and thus by the triad °e, and co-triad °cu'. Then, 
the phase space has coordinates (c, p), which are conjugate variables satisfying {c, p} = 
This connection formalism is related to the metric one via the relations 

\p\=a 2 , c = ^{k + ja), (5.13.2) 

where k = 0, ±1 is the usual curvature parameter. Since the Gauss and the diffeomor- 
phism constraints are already satisfied (with (|5.13.1[) they identically vanish), the scalar 
constraint 

U^V) = -\\J\p\ (Jp(c - T) 2 + T 2 ^ 0, (5.13.3) 
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where r ot k, is the only remaining. This is nothing but the Friedmann equation and re- 
duces to a simple form in the flat case: Hj[^q = — -^ I c 2 \/\p\. We stress that p G (—00, +00) 
and the classical singularity appears for p = 0. Furthermore the changes in the sign of 
p correspond to the changes in the orientation of the physical triad ef, related to °ef via 
sign(p)/7p. 

The quantization of this model follows the lines of LQG and therefore we have to con- 
struct SU(2) holonomies and fluxes, in order to promote these variables to quantum op- 
erators, which become 

h e (c) =cos(^) +2(e a x)^sin(^£) , P s (p) = A s p , (5.13.4) 

where A$ is a factor determined by the background metric and y. G (— oo, +oo) is a real 
continuous parameter, along which the holonomies are computed. The elements of the 
holonomies can be recovered from the almost periodic functions N^(c) = e' ? ' C//2 and then 
the cylindirc functions of this reduced model are given by 

g(c)=££;^ c/2 eCyl s , (5.13.5) 
;' 

Cyl$ being the space of the symmetric cylindric functions. The holonomy flux algebra is 
generated by e 1 ^ ^ 1 and by p, i.e. we are in a hybrid representation between the Heisen- 
berg (p) and Weyl (e MX ) ones. 

The Hilbert kinematical space is now obtained requiring that the N^(c) form an or- 
thoromal basis, i.e. (N U \N U >) = Sum'/ in analog}@ with the scalar product (|5.12.10|) of 
the full theory. From general theoretical considerations, the Hilbert space is necessarily 
'Ks = ^ 2 {^Bohr,d}i), where RBohr is a compact Abelian group (the Bohr compactification 
of the real line). This is the space of the almost periodic function. 

Let us stress that the construction of the Hilbert space is the key point of all the theory 
and of its main results. In fact, using RBohr one can introduce a new representation of 
the Weyl algebra. This is inequivalent to the standard Schrodinger representation and 
therefore leads to different results from the WDW theory. 

In this Hilbert space the operators N« and p act by multiplication and derivation, re- 
spectively. As usual, let us introduce the bra-ket notation N fl (c) = (c\ji). The volume 
operator (V = |p| 3/2 ) has a continuous spectrum V\y) = V u \fi) <*■ |p| 3/2 Zp|p), in contrast 
to Lod^l. In the reduced theory the spectrum is discrete in a weaker sense: all the eigen- 
vectors are normalizable. Hence the Hilbert space can be expanded as a direct sum, rather 
than as a direct integral, of the one-dimensional eigenspaces of p. 

In view of the analysis of the singularity, we have to address the inverse scale factor 
operator which is a fundamental one because, at a classical level, the inverse scale factor 
sign(p)/ yjp] diverges toward the singularity. Let us express it in terms of holonomies 
and fluxes and then proceed to the quantization. We note that the following classical 



Note that this is a Kronecker-delta rather than the usual Dirac one. 
13 This feature can be attributed to the high degree of symmetry. In fact, in LQG the spin networks are 
characterized by a pair (£,j) consisting of a continuous edge (.. and a discrete spin j. Due to the symmetry, 
such pair now collapses to a single continuous label \i. 
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identity 

^ = ^Tr(£A-{^,^}), (5.13.6) 

holds, where the holonomy hj is evaluated along any given edge. In fact, since we have 
hihj 1 , the choice of the edge is not important, i.e. we do not introduce a regulator and the 
expression (|5.13.6)> is exact. With the scalar constraint the situation will be different. 

We can proceed with the quantization of operator (|5.13.6)> in a canonical way. The eigen- 
values are given by 

(^ Wot _L(^/3_^3) W , (5.13.7) 

where V„ is the volume operator eigenvalue defined above, whose fundamental proper- 
ties are to be bounded from above and to coincide with the operator 1/ \f\p~\ for ^> 1. 
The upper bound is obtained for the value ]i = 1 and therefore^ \ p\ m li 2 ~ lp 1 - 

The physical situation emerging is intriguing. Although the volume operator admits a 
continuous spectrum and a zero volume eigenstate (the \ji = 0) state), the inverse scalar 
factor is non diverging at the classical singularity, but is bounded from above. This can 
indicate that, at a kinematical level, the classical singularity is avoided in a quantum ap- 
proach. The semiclassical picture, i.e. the WDW behavior of the inverse scalar factor, 
appears for 3> 1 and therefore far from the fully quantum regime. Such a behavior 
contrasts with the WDW formalism where the inverse scale factor is unbounded from 
above and the differences reside in the non-standard Hilbert space adopted. In fact, dif- 
ferently from the WDW theory, all eigenvectors of p are normalizable in LQC, including 
the one with zero eigenvalue. Therefore, in order to define an inverse scale factor opera- 
tor, one has to formulate the alternative procedure discussed above. 

Let us stress that a key ingredient is the absence of the inhomogeneous fluctations 
which conversely are present in the full theory. In fact 1 100 J computed in full LQG that 
the analogue of the inverse scale factor is unbounded from above on zero volume eigen- 
states. The boundedness of the inverse scale factor is neither necessary nor sufficient for 
the curvature singularity avoidance. 

As a last point we have to impose the scalar constraint at a quantum level to discuss 
the fate of the singularity from a dynamical point of view. 

In order to follow the lines of the the full theory, the starting point will be the constraint 
(|5.12.3b|) and not the reduced one (|5.13.3|l . In the expression (|5.13.3|) , the connection c 
itself is present rather than the holonomies and the operator H( C 'P) is not well defined at 
a quantum level. Let us investigate, for simplicity, the flat case. The two parts of (|5.13.3|) 
are then proportional to each other. Mimicking the procedure followed in the full theory 
|485] we rewrite the constraint in terms of holonomies and fluxes as 

< P) = FT E^ Tr (WiQf)' 1 ' V}) + O(cVo) , (5.13.8) 

Kr y v-o ijk K ' 

where h f t '° denotes the holonomy computed around a square ol^ with each side having 



We recall that in this model the classical Ricci scalar curvature is given by R ~ I /a 2 and therefore, from 
equation ( 15.13.71 . at the classical singularity it doe not diverge, assuming the value R ~ 1/Zp. 
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length jio, i.e. hij = h\hjh i l h- 1 . Differently from expression (|5.13.6|) , the dependence on 
jio does not drop out and now plays the role of a regulator. However, at a classical level, 
we can take the limit jio — > and verify that the resulting expression coincides with the 
classical Hamiltonian (|5.13.3)> . As in the full theory, the problems arise at a quantum level. 

Let us investigate the action of the quantum operator ftjfc on the eigenstates which is 

A i P) \^ = 133 (Vkc " V^oXI^ + ^o) - 2\ H ) + \n-4jio)) • (5-13.9) 
K i ro 

The limit }io —> fails to exist and the classical regulated version of equation (|5.13.3|) be- 
haves as H^'^ — > H( c, p\ Since it contains c, a limit at a quantum level does not exist 
because the operator c itself does not exist in the Hilbert space J{$ ■ As a matter of fact, 
in this reduced theory, there is no way to remove the regulator. This is resolved in com- 
parison with the full theory. If we assume that the predictions of LQG are true, then the 
regulator jiQ can be shrunk until the minimal admissible length given by the area operator 
spectrum. In this sense, the }Iq — > limit is physically meaningless. However, how this 
reduced theory (LQC) can see a minimal length coming out from LQG is not fully under- 
stood, since it is not the cosmological sector, but the usual cosmological minisuperspace 
phase space quantized through the LQG methods. 

Let us investigate the physical states. As in the full theory, the physical states are those 
annihilated by all the constraints and live in some bigger space Cyl* s and they do not need 
to be normalizable. A generic stat^l (Y| S Cy/g can be expanded as 

(Y| = X>0*/*)04 (5.13.10) 

f 

where (p represents a generic matter field and satisfies the constraint equation 

(Y| (h^' p) +hJ ) + = (5.13.11) 
and therefore the equation for \p{ji,<p) 

(V^o - V^Mv + 4^ ,tf>) - 2( V w " (5.13.12) 

has to hold. This is nothing but a recurrence relation for the coefficients tp(fi,(p) which 
ensures that (Y| is a physical state. Even though ]i is a continuous variable, it is given 
by (|5.13.12|) as an algebraic rather then by a differential equation as, for example, in the 
WDW theory. 

Let us now briefly discuss what happens to the classical singularity at a dynamical 
level. It corresponds to the state \\i = 0) and as we can see from equation (|5.13.12|l , 
starting at ]i = — 4Nf/o we can compute the coefficients i/?(4^o (n — N), (p) for n > 1, all but 
ip(fi = 0, </>), because the generic coefficient vanishes if and only if n = N. Although the 
quantum evolution seems to break down right at the classical singularity, this is not the 
case, since the coefficient xp(}i = 0, (p) is decoupled from the others thanks to Vp — VL^ = 

15 The notation (Y| is adopted to the eventuality of non-renormalizable states. 
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and to H^ ip(ji = 0,(p) =0 realizes. Therefore the coefficients in (|5.13.12|) are such that 
one can unambiguosly evolve the states through the singularity even though \p(ji = 0,(p) 
is not determined and the classical singularity is solved in the LQC framework. 

We conclude reporting that the previous results have been significantly extended |24l 
with a rigorous formulation of the physical Hilbert space, of the Dirac observables and of 
the semi-classical states. Furthermore, taking a massless scalar field as a relational clock 
for the system, the classical Big-Bang is replaced by a quantum Big-Bounce displaying in 
detail an intuitive picture of the Universe evolution in the Planckian era is displayed. 



5.14 Mixmaster Universe in LQC 



We review the basic aspects of the Mixmaster Universe in the LQC framework in this 
Section, while for the detailed analysis we demand to the literature |93ll95ll94| . 

As we have seen, the Bianchi IX evolution toward the singularity sees infinite sequences 
of Kasner epochs characterized by a series of permutations as well as by possible rotations 
of the expanding and contracting spatial directions. However, this infinite number of 
bounces within the potential, at the basis of chaos, is a consequence of an unbounded 
growth of the spatial curvature. When the theory offers a cut-off length and the curvature 
is bounded, the Bianchi IX model naturally shows a finite number of oscillations and in 
LQC a quantum suppression of the chaotic behavior appears close to the singularity. 

Let us formulate in the connection formalism the vacuum Bianchi IX model. The spatial 
metric dl 2 = a 2 (co 1 ) 2 can be taken diagonal, leaving three degrees of freedom only. The 
basic variables for a homogeneous model are 

A\ = c (l) A\w{, Ef = pWAjXf, (5.14.1) 

where CO are the left-invariant 1-forms, X; are the left-invariant fields dual to to (a/(Xj) = 
Sj), the SO(3)-matrix A contains the pure gauge degrees of freedom and I,J,K,..., run 
as 1, 2, 3. The physical information are expressed in terms of the gauge invariants Cm and 

which satisfies the Poisson brackets {cj, pJ} = jkSI. These variables are related to 
the scale factors flj, the spin connections Tj and the extrinsic curvature Kj = — x lia\ by the 
relations 

p l = sign(a 7 ) |fl 7 fl K |sign(fli), cj = Tj - 7X7, (5.14.2) 

where 

r i(fi + ££_A) i(E; + 4_£«. (5 . 14 .3) 

2 \a K aj aja K J 2\pJ p K (p 1 ) 2 J 
The classical dynamics is governed by the scalar constraint (|5.12.3b|) which reads as 



K 



(TjY K - Tj) flj - -aiaja K + cyclic 



. (5.14.4) 



The particular case of the isotropic model, the open (k = 1) FRW, can be recovered setting 
d\ = «2 = fl 3 = a - On the other hand, the Bianchi I model can be obtained for T/ = 0. The 
potential term from (|5.14.4[) is given by 



W(p) = 2 ( p l p ] (TiTj - T K ) + cyclic ) (5.14.5) 
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which has infinite walls at small p due to the divergence of the spin connection compo- 
nents. 

The loop quantization is performed straightforwardly as in the isotropic case. In fact, 
an orthonormal basis is given by the p 7 -eigenstates 1^1,^2,^3) = <8> ^2) <8> ^3) and 
the Hilbert space is taken as a direct product of the isotropic ones, it is separable and is a 
subspace of the kinematical non-separable Hilbert space. As usual, the cylindric functions 
are given by a superposition of the functions (c\y) ~ exp( ! > c /2) and the basic quantum 
operators are the gauge invariant triad operators p 1 (fluxes) and the holonomies 

hi{c) = cos (I) + 2A< t ; - sin (|) . (5.14.6) 

The p 1 and hj act as derivative and multiplication operators, respectively. In particular, 
the volume operator defined from p 1 as V = y/\p l p 2 p 3 \, acts on the eigenstates | ^3) 
as 

v\m, n, n) = {2 K if /2 \J\m^m\ \n> y*j n) > (5.14.7) 

and it has a continuous spectrum. 

From these basic operators we can obtain, similarly to the isotropic case, the inverse 
triad operator. Since it is diverging toward the classical singularity, we are interested 
to its behavior at a quantum level. Conceptually its construction in the homogeneous 
cosmological sector is the same as in the isotropic one, the only difference residing in the 
computations and in the appearance of some quantization ambiguities, as in particular 
the half-integer j and a continuous parameter / G (0, 1). Nevertheless, all the results are 
independent on them. This behavior mimics the quantization ambiguities present either 
in LQG or in the isotropic sector of LQC. While in the full theory we find an ambiguous 
choice of the spin number j associated to a given edge of the spin network, in the isotropic 
LQC it is reflected on the choice of the Hamiltonian regulator Nevertheless, this is but 
a parameter that we can neither shrunk to zero, neither fix in some way in the context 
of LQC theory itself but it outcomes from keeping some prediction of another theory 
(usually LQG). In the inverse scale factor of isotropic LQC none quantization ambiguity 
appears. 

Since we can express | p 1 | _1 in terms of holonomies and volumes via a classical identity, 
it can be canonically quantized as follows 

\v l V x \v\,m,n) = Ajjfjjim)]^,^,^) , (5.14.8) 

where 




(5.14.9) 



and Ai\ is a function of the quantization ambiguities, i.e. j and I. The values fj f i(}i) 

decrease for \i < 2j and we have /y,/(0) = and therefore = 0) = 0. Thus, 

the inverse triad operator annihilates the state corresponding to the classical singularity 

\m = o). 

Fundamental properties of the eigenvalues of the operator (|5.14.8|) can be extracted 
from the asymptotic expansions of Fn(^j) = Aj r ifii(}ii). For large Fu(}ii) = F/(vj), 
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with Vi = jii/2], and no dependence on; appears, in particular 

1/vj ifvj>l (/*!>/') 



vi 
l + l 



Vi-' (5.14.10) 



j if vj < 1 (f*i < ;') . 



The classical behavior for the inverse triad components is obtained for Vi > 1 and the 
loop quantum modifications arise for Rewriting the spin connection in the triad 

representation, the potential (|5.14.5|) is given by 

W h , (v) = 2( 7 k/) 2 (vjv; (TiTj - T K ) + cyclic) , (5.14.11) 

where 

t i( v i) = 2 [ v Ksign(vi)F z (vj) +v;sign(v K )F;(i/ K ) - vjv K F? (vj)] . (5.14.12) 

The loop quantization of the scalar constraint is different in the homogeneous case. 
In fact, unlike the full theory, the spin connections are tensors due to the homogeneity 
symmetry and cannot vanish as in LQG, thus the holonomies will depend also on them 
and the quantum scalar operator leads to a partial difference equation, like in isotropic 
LQC, for which we refer to I95fl. 

Let us discuss the Mixmaster Universe in the LQC framework, analyzing its behavior at 
a semiclassical level, considering the modifications induced in the dynamics by the loop 
quantization. Since an effective Hamiltonian must arise from the underlying quantum 
evolution, we will proceed in two steps. 

z) For a slow varying solution, the differences equation is specialized, in the contin- 
uum regime \i\ ~ p 1 / jk 3> 1,0 thus obtaining the Wheeler-DeWitt equation 

^ p y ^VIPW| S(£) + +W jJ (p)y/\pYf\S(p) = 
where (p denotes a generic matter field. 



ii) We take the WKB limit of the wave function T(p) = y / \p 1 p 2 p 3 \S(p), i.e. T ~ e tA/h , 
leading to the Hamilton -Jacobi equation for the phase A to zero-th h order. 

We have obtained the classical dynamics plus the quantum loop corrections. The key 
point for the classical analysis of the effective Hamiltonian is that the classical region 
fij > 1 can be separated in two subregions. In fact, remembering the p-dependence in 
the WDW equation (|5.14.13|) given by ifr/2; ~ p 1 / jyK, we obtain the condition jij » 1 
for p 1 < jyK (j > Hi > 1) and p 1 > jjK > /' > 1). 

The second subregion 3> j 1 3> 1) is the purely classical one, i.e. where the Misner 
picture is still valid. From expansion (|5.14.10[) the eigenvalues of the inverse triad opera- 
tor correspond to the classical values. On the other hand, the region where / > /Jj > 1 is 



16 The eigenvalues of the triad operator p 1 are given by p I \}ii,}i2't i 3) — P 1 \V-ltV-2>V-3)i where p 1 ~ JKfij. 
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characterized by loop quantum modifications. In fact, the inverse triad operator eigenval- 
ues have a power law dependence. The quantum modifications of the classical dynamics 
are controlled by the parameter j: if it is large enough, one can move the quantum effects 
within the effective potential into the semi-classical domain. On the other hand, the WKB 
limit (h — > 0) is strictly valid in the Misner region (jij 3> ; 3> 1), because in the first region 
a dependence on (77c) appears in the potential term. We assume the validity of such 
approximation in the (j 3> ]i\ ^> l)-region, because the inverse triads vanish as p 1 — > 0. 

A qualitative study of the modified classical dynamics, arises from analyzing the po- 
tential term, whose explicit expression is more complex then the original one, making the 
analysis of the dynamics tricky. The volume variable is regarded as a time variable and 
in general the dependence on it does not factorize. Nonetheless, in the second region, the 
Misner potential is restored. 

Considering the particular case c\- = (the Taub Universal ), we can qualitatively 
study the effective (loop) potential. In this case we get 1/2 = V3 = v and v\ = p = 
V 2 / ((2j) 3 v) and therefore T2 = T3. The wall is seen for v ^> 1 so that F/(t/) ~ v^ 1 , but p 
not negligible, i.e. the relation F; (p) ~ p 2 holds and the potential wall (|5.14.11|) becomes 



where V denotes the volume. For p » 1, Flip) — P 1 and the classical wall e 4a_8 / 3 + 
is restored. The key difference with the Misner case relies in the wall finite height. As 
the volume decreases, the wall moves inwards and its height decreases as well. In the 
following evolution the wall completely disappears. In particular it gets its maximum for 
j8+ = —a and vanishes as e 12u <x V 4 toward the classical singularity (oc — ► —00). 

This peculiar behavior shows that the Mixmaster evolution stops at a given time and 
therefore chaos disappears. In fact, the point-Universe, when the volume is so small that 
the quantum modifications arise, will never bounce against the potential wall and the 
Kasner epochs will continue without any replacement. This behavior predicted by the 
LQC framework produces (qualitatively) the same results of the influence of a scalar filed 
on the Universe dynamics: in fact, also in that case, at a given time the particle performs 
the last bounce and then it freely moves (see Section l4.13j) . 

For a discussion of the Bianchi cosmologies in this new canonical first order formalism 

see iiwira i^ i^i2Mi^ro iM ir^iTTfl 

5.15 On the GUP and the Minisuperspace Dynamics 

This paragraph is devoted to explain some results obtained in a recent approach |45j |46) 
(see also 14961 14011 1497| ) to quantum cosmology, in which the notion of a minimal length 
naturally appears. In particular, the purpose is to quantize a cosmological model by using 
a modified Heisenberg algebra, which reproduces a Generalized Uncertainty Principle 



V 4 



Ffp^-lpF^p)) , 



(5.14.14) 




(GUP) 




(5.15.1) 



The Taub model is nothing but a particular case of the Bianchi IX one in which /5_ = 1 446 [. 
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where /3 is a "deformation" parameter. The uncertainty principle (|5.15.1|) can be obtained 
by considering an algebra generated by q and p obeying the commutation relation 

[q,p] =/(l + ft> 2 ) . (5.15.2) 

Such deformed Heisenberg uncertainty principle appeared in studies on string theory 
|227U508l l7| and leads to a fundamental minimal scale. From this point of view, a minimal 
observable length is a consequence of the limiting scale given by the string one. However, 
we stress that the minimal scale predicted by the GUP is different from the one predicted 
by other approaches. In fact, equation (|5.15.1|) implies a finite minimal uncertainty in 
position Aq m j n = y / j6[291J. This way, we will introduce a minimal scale in the quantum 
dynamics of a cosmological model. 

The appearance of a non-zero uncertainty in position poses some difficulty in the con- 
struction of a Hilbert space, as position eigenstates cannot be physically constructed. An 
eigenstate of an observable necessarily has to have a vanishing uncertainty. Although 
it is possible to construct position eigenvectors, they are formal but not physical states. 
In order to recover information on position, we consider the so-called quasiposition wave 
functions 

m - N C (i4W'> exr ('7r rcian ( v/ft ')) * (p) - <5 - 153) 

where £ is the quasiposition defined by the main value of the position q on appropriate 
functions, i.e. (q) = 1 and N is a normalization factor. The quasiposition wave function 
(|5.15.3|l represents the probability amplitude to find a particle maximally localized around 
the position £, i.e. with standard deviation Aq m j n . For more details on quantum mechanics 
in this framework we recommend ||291H290l| . 

The GUP approach relies on a modification of the canonical prescription for quantiza- 
tion and can be applied to any dynamical system. Moreover, such formalism allows us to 
analyze some peculiar features of string theory in the minisuperspace dynamics. 

Let us investigate the consequences of a Heisenberg deformed algebra (|5.15.2|) of the 
quantum dynamics of the flat (k = 0) FRW model in presence of a massless scalar field (p. 
We are interested on the fate of the classical singularity. In the following we summarize 
the discussion and results reported in [45. 46J. 
The Hamiltonian constraint for this model has the form 



2 

H g rav + H <p = -9k P 2 x x + y ~ ' x = a 3 , (5.15.4) 



where a is the scale factor. In the classical theory, the phase space is four-dimensional, 
with coordinates (x, p x ; <p, p§), and the physical volume of the Universe vanishes at x = 
and the singularity appears. Moreover each classical trajectory can be specified in the 
(x,<p) -plane, i.e. (p can be considered as a relational time. The dynamical trajectories read 
as 

x 



1 

<b = ± — = ln 



+ <po, (5.15.5) 



where Xq and (po are integration constants. In equation d5.15.5l) , the plus sign corresponds 
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to a Universe expanding from the Big-Bang, while the minus sign to a contracting one into 
the Big-Crunch. The classical cosmological singularity is for <p = ±00 and every classical 
solution, in this model, reaches it. 

The canonical approach (the WDW theory) does not solve the singularity problem. In 
fact, we can construct a state localized at some initial time. Then, in the backward evolu- 
tion, the peak of the wave packet will move along the classical trajectory (|5.15.5|) falling 
in the classical singularity! 86 1, so that the latter is not tamed by quantum effects. 

This picture is radically changed in the GUP framework and the modifications can be 
realized in two steps. Firstly, we can show how the probability density |Y(£, t)\ 2 to find 
the Universe around £ ~ (around the Planckian region) can be expanded as 

\Y(C,t)\ 2 ^\A(t)\ 2 + ?\B(t)\ 2 , (5.15.6) 

where f is a dimensionless time t = 3\/x(p and the wave packets 

poo 

Y(£f) = / de g(e)V £ (Oe t£t (5.15.7) 
Jo 

are such that the state is initially peaked at late time, i.e. the weight function g(e) is a 
Gaussian distribution centered at some e* < 1 (at energy much less then the Planck one 
1/lp). Of course, Y e (£) rapresents the quasiposition eigenfunctions (|5.15.3|) of this problem. 

Near the Planckian region, the probability density to find the Universe is | A(t) \ 2 , which 
is very well approximated by a Lorentzian function peaked in t = 0, which corresponds 
to the classical time when x(t) = Xq. For Xq ~ 0(Zp), the probability density to find 
the Universe in a Planck volume has a maximum around the corresponding classical val- 
ues and vanishes for t — ► —00, where the classical singularity appears. In this sense the 
classical cosmological singularity is solved. 

The more interesting differences between the WDW and the GUP approaches relies on 
the wave packets dynamics. If we consider a wave packet initially peaked at late time, 
once numerically evolved backward, the integration results in a probability density, at 
different fixed values of £, still approximated by a Lorentzian function. The width of 
this function remains the same as the states evolve from large £ (~ 10 3 ) to £ = 0, while 
its peaks move along the classically expanding trajectory d5.15.5P for values of £ larger 
than ~ 4. Near the Planck region, i.e. when £ G [0,4], we observe a modification of the 
trajectory of the peaks since they follow a power-law up to £ = 0, reached in a finite time 
interval and they escape from the classical trajectory toward the singularity (see Figure 
I5.7|) . The peaks of the Lorentzian at fixed time f slowly evolve remaining close to the 
Planck region. Such behavior outlines how the Universe has a stationary approach to the 
cutoff volume, accordingly to Figure \57\ and is different from other approaches. 

In fact, it was shown how the classical Big-Bang is replaced by a Big-Bounce in the 
framework of LQC |24| . Qualitatively, one expects that the bounce (and the consequently 
repulsive features of the gravitational field in the Planck regime) are consequence of a 
Planckian cut-off length, but this is not the case. We can observe from Figure 15.71 that 
there is not a bounce for the quantum Universe. The main difference between the two 
approaches resides in the quantum modification of the classical trajectory. In fact, in LQC 
we observe a "quantum bridge" between the expanding and contracting Universes, while 
in our approach the probability density of finding the Universe reaches the Planck region 
in a stationary way. 
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Figure 5.7: The peaks of the probability density |Y(£,£)| 2 are plotted as functions of t 
and ln(£). The points (resulting from numerical computation) are fitted by a logarithm 
0.050 ln(£) + 0.225 for £ > 4 and by a power law 0.067£ 1060 for £ e [0, 4] (from HBI ). 



Finally, we emphasize two points: i) the above result on the singularity-free behavior of 
the GUP (FRW) Universe is confirmed by the more general Taub model (for a discussion 
about the Taub Universe in the GUP framework we refer to EH). Also in this case the 
singularity is probabilistically suppressed, ii) In the GUP Taub Universe, the wave pack- 
ets favour the establishment of a quantum isotropic Universe with respect to those in the 
WDW theory. 

5.16 Quantum Chaos 

In the previous section we have widely discussed the chaotic features of the Mixmaster 
dynamics; as usual in chaotic systems, we expect that this classical behavior is maintained 
also at a quantum level. Several attempts were made over the years to characterize the 
"quantum chaos" of the Bianchi IX model, analyzing the evolution of the shape of the 
wave functions I205ll206ll63l or the distribution of the energy levels fl58lf220l[2T9l. 

The numerical simulation performed in [205 , 206J is focused on a wave packet initially 
peaked in the center of the (/3+, /3_) plane and is a solution of the WDW equation. The 
evolution is numerically computed as a multiple-reflection problem of the wave-packet 
within moving potential walls. The quantized Mixmaster model possessing a quantum 
mixing property can also be seen as a quantum chaotic system. The mixing operates on 
the expectation values of the scale factors, which tend to an equilibrium value without 
any help of the matter field or of the cosmological term. A self-similar structure of the 
distribution of probability density emerges, strongly suggesting the existence of a fractal 
structure. 

From [63J one finds another feature of the classical chaos emerging at a quantum level. 
Adopting a path-integral representation of the dynamics, a Monte Carlo simulation of the 
dynamics shows how the anisotropy potential at the origin of the classical chaos yields 
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a quantum mechanical correlation between large universe volume and high anisotropy, 
although the peak at zero anisotropy is significant as well. This result suggests a possible 
failure of the quantized Mixmaster cosmology to represent the very Early Universe, or it 
might be a consequence of neglecting non gravitational fields. 

A different point of view is addressed in ||158l I220| , where the quantum chaos is out- 
lined when studying the distribution of the energy levels of the quantized dynamics, as- 
suming Dirichlet boundary conditions. In the Poincare upper half -plane representation, 
the authors evaluate approximately the first 10 3 eigenvalues, theie distribution and how 
the energy levels are influenced by the deformation of the circular side of the bounding 
triangle. 



124 



6 Inhomogeneous Mixmaster Model 



This Section is devoted to the generalization of the Mixmaster dynamics to the inhomo- 
geneous case and to its multidimensional extensions. We firstly provide the construction 
of the generic cosmological solution of the Einstein equations in the neighborhood of the 
Big Bang and then we discuss the so-called fragmentation process. The Hamiltonian for- 
mulation of the inhomogeneous Mixmaster is reviewed with particular attention to the 
role played by the Ricci scalar and the quantum picture is then briefly analyzed. We close 
this review showing how the chaos is a dimensional phenomenon summarizing the most 
important results obtained in the multidimensional extensions of GR. 

6.1 Formulation of the generic cosmological problem 

We will focus our attention on the generic behavior of the Einstein equations near the Big 
Bang, i.e., we will discuss the properties of their general cosmological solution, meaning 
with it a solution that possesses the right number of physical arbitrary functions allowing 
to specify arbitrary initial conditions on a generic non-singular space-like hypersurface 
at a fixed instant f, thus four arbitrary functions in vacuum, or eight if a perfect fluid is a 
source |319](for a discussion in presence of inflationary matter see [489]). 

Belinskii, Khalatnikov ans Lifshitz in the 70's derived such a solution and how its dy- 
namics resembles the one of the homogeneous indicess of types VIII and IX |5ll l56l |57| 
(see also |73l |70| ). The construction can be achieved firstly considering the generic so- 
lution for the individual Kasner epoch, and then providing a general description of the 
alternation of two successive epochs. The answer to the first question is given by the so- 
called generalized Kasner solution, while to the latter is found to be in close analogy to 
the replacement rule in the homogeneous indices 

For a general view on inhomogeneous cosmologies, see |108t 14591 l3l 11091 124TI 14381 14391 
14401 144TI 13001 l72l l405ll . 

6.1.1 The Generalized Kasner Solution 

The homogeneous indices of types VIII and IX provide the prototypes for the construction 
of the generic (non-homogeneous) solution of the Einstein equations in the neighbour- 
hood of the singularity. In [294], it is shown that the Kasner solution can be generalized 
to the inhomogeneous case and near the singularity, as 




(6.1.1) 



where 



b ~ 



m 



c ~ t Pn , 



(6.1.2) 
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and pi f , p m , p n are functions of spatial coordinates subjected to the conditions 

Vl {xf) + p m (x y ) + Pn(x y ) = pjV) + Pi(x y ) + P«(* 7 ) = 1 ■ (6-1-3) 

Differently from the homogeneous indicess, the frame vectors I, m, n are now arbitrary 
functions of the coordinates (subject to the conditions imposed by the Einstein equations 
(k|29U). 

The behaviour (|6.1.1|) cannot last up to the singularity, unless a further condition is im- 
posed on the vector 1 ( i.e., the one corresponding to the negative index p\) 

I . V A / = . (6.1.4) 

This condition reduces to three the number of arbitrary functions, i.e. one less than the 
generic solution: in fact, (|6.1.1[) possesses twelve arbitrary functions of the coordinates 
(nine components of the Kasner axes and three indexes pj(x 7 )), and satisfies the two 
Kasner relations (|6.1.3|) , the three Ooc Einstein equations, the three conditions from the 
invariance under three-dimensional coordinate transformations, and (|6.1.4[) . 

The generalized Kasner solution (|6.1.1|) does not require any additional condition, since 
the individual Kasner epoch lasts only for a finite interval of time, and contains, for ex- 
ample, four arbitrary functions in vacuum. 



6.1.2 Inhomogeneous BKL indices 

Let us now generalize our scheme by investigating the implications of the remaining con- 
ditions (|6.1.4[) . This analysis leads to the inhomogeneous BKL indices. 

We assume that the factors that determine the order of magnitude of the components 
of the spatial metric tensor ((6.1.1)> can be included in the functions a, b, c, which change 
with time according to (|6.1.2|) , i.e. the vectors /, m, n define the directions of the Kasner 
axes. For non homogeneous spaces there is no reason to introduce a fixed set of frame 
vectors, which would be independent of Kasner axes. 

The time interval of applicability of solution ( |6.1.1| ) is determined by conditions which 
follow from the Einstein equations. Near the singularity, the matter energy-momentum 
tensor in the 00- and a/5-components may be neglected 



^=3 t (^)+P^ = 0. (6.1.5b) 



1K « 



Solution (|6.1.1|) - (|6.1.3[) is obtained neglecting the three-dimensional Ricci tensor in 
d6.1.5bl ) and its validity is easily formulated in terms of the projections of the tensors along 
the directions Z, m, n [294j, which must satisfy the conditions 

Pi K. K « f ~ 2 , A » K. Pn. (6-1-6) 

The off -diagonal projections of (|6.1.5b|) determine the off-diagonal projections of the met- 
ric tensor (t]i m/ r]\ n , r\ mn ), which should only be small corrections to the leading terms of 
the metric, as given by (|6.1.1|) . In the latter, the only non- vanishing projections are the 
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diagonal (rjn, tj mm , t] nn ), and satisfy |5Z| 

Vim < yJVllVmm > K]\n < y/VllVnn , Vmn < y/tjmmVnn , (6.1.7) 

and consequently 

Pl m < ab/t 2 , Pi„ < ac/t 2 , P mn < be ft 2 . (6.1.8) 

In this case, the off-diagonal components of (|6.1.5b|) can be disregarded to leading order. 
The Ricci tensor P% for the metric (|6.1.1|) is given in Appendix D of |294fl . 
The diagonal projections P\, P%, P" contain the terms 

If alV A(al) \ k 2 a 2 k 2 a 4 

- - (6.1.9) 



2 \abc(l ■ [m x n]) J b 2 c 2 A 2 P ' 

and analogous terms with a I replaced by b m and c n, where 1/k denotes the order of 
magnitude of spatial distances over which the metric significantly changes, and A is the 
same as in (|3.6.5D . In d6.1.9|) , all the vectorial operations are performed as in the Euclidean 
case. 

According to (|6.1.6|) , we get the inequalities 



flv / fc7A<l, K/fcM<l, cv^7A<l, (6.1.10) 

which are not only necessary, but also sufficient conditions for the existence of the gener- 
alized Kasner solution. After conditions (|6.1.10|) are satisfied, all other terms in Pj, P™, P%, 
as well as P\ m , P\ ni P mn , automatically satisfy (|6.1.6|) and (|6.1.8|l as well. 
An estimate of these terms leads to the conditions 

k 2 

[a 2 b 2 a 3 b,..., a 2 bc,...) < 1 , (6.1.11) 

containing on the left-hand side the products of powers of two or three of the quantities 
which enter in (|6.1.10|) , and therefore are a fortiori true if the latter are satisfied. Inequalities 
(|6.1.11|) moreover represent a natural generalization of those which already appearing in 
the oscillatory regime in the homogeneous case. 

As t decreases, an instant t tr may eventually occur when one of the conditions (|6.1.10|) 
is violated (the case when two of these are simultaneously violated can happen when the 
exponents p\ and p2 are close to zero, corresponding to the case of small oscillations). 
Thus, if during a given Kasner epoch the negative exponent refers to the function a(t), 
i.e., pj = pi, then at t tr , we have 

atr = ^~l. (6.1.12) 

Since during that epoch the functions b(t) and c(t) decrease with t, the other two inequal- 
ities in (|6.1.10|) remain valid and at t ~ ttr we shall have 

btr < «ir / Ctr < «fr • (6.1.13) 

At the same time all conditions (|6.1.11|l continue to hold, and all off -diagonal projections 
of d6.1.5bl) may be disregarded. In the diagonal projections (|6.1.9|) , only terms containing 



127 



6 Inhomogeneous Mixmaster Model 



fl 4 / t 2 become relevant. In such remaining terms we have 

(a I • V A (a I)) = a (I • [Vfl x /]) + a 2 (I • V A /) = a 2 (I • V A /) , (6.1.14) 
i.e., the spatial derivatives of a drop out. 

As a result, we obtain the following equations for the replacement of two Kasner epochs 

(6.1.15a) 
(6.1.15b) 
(6.1.15c) 
(6.1.15d) 

which differ from the corresponding ones of the homogeneous indices d3-6.2H3.6T3l) only 
for the quantity 

A(x) = /; VA? 1 , (6.1.16) 
I ■ [m x n\ 

no longer being a constant, but a function of the space coordinates. Since d6.1.15cl) is a 
system of ordinary differential equations with respect to time where space coordinates 
enter parametrically only, such difference does not affect at all the solution of the equa- 
tions and the following map. Similarly, the law of alternation of exponents derived 
for homogeneous indices remains valid in the general inhomogeneous case. Numeri- 
cal evidences 1 506 1 support the point-like dynamics of the generic cosmological solution 
at the basis of the BKL conjecture. 





(abc)- 
abc 


i i2 « 2 

+ A 2b 2 c 2 


= 0, 


_ R m = 

m 


(abc)- 
abc 


2 a 2 
2b 2 c 2 


= 0, 


-R\ = 


(abc)- 
abc 


2 a 2 
2b 2 c 2 


= 0, 


-R°o = 


a b 
a b 


+ ^ = 0, 
c 





6.1.3 Rotation of the Kasner axes 

Even if the dynamics is quite similar to that of the homogeneous indicess in vacuum, the 
new feature of the rotation of the Kasner axes emerges. If in the initial epoch the spatial 
metric is given by d6.1.1|) , then in the final one we have 

Kp = aXl'p + fc 2 m>; + c 2 n>J , (6.1.17) 

with a, b, c given by a new set of Kasner indexes, and some vectors V , w! , n'. If we project 
all tensors (including h a a) in both epochs onto the same directions /, m, n, the turning of 
the Kasner axes can be described as the appearance, in the final epoch, of off-diagonal 
projections Y\\ m , f]j n , r] mn , which behave in time as linear combinations of the functions 
a 2 , b 2 , c 2 . Such projections do indeed appear |57]| , and induce the rotation of the Kasner 
axes. 

The main effects can be reduced to a rotation of the m— and n— axis by a large angle, 
and a rotation of the I— axis by a small one, thus neglecting the small changes of I the 
new Kasner axes are related to the old ones as 

l' = l, m' = m + a m l, n' = n + a n l, (6.1.18) 
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where the o~ m , a n are of order unity, and are given by 

<r m = ^-([Ixjnl-V^ + ^wVAlU T , 

Vi + 3pi I A A ) I ■ [mxn] 

<r n = V- I [" x • - ^-n ■ V A A — ^ T , (6.1.19a) 

P2 + 3pi \ 1 J A A J I -[mxn]' y ' 

and can be inferred also from the — a Einstein equations which play the role of con- 
straints to the space functions. 

The rotation of the Kasner axes (that appears for a matter-filled homogeneous space 
only) is inherent in the inhomogeneous solution already in the vacuum case. The role 
played by the matter energy-momentum tensor can be mimicked by the terms due to in- 
homogeneity of the spatial metric in the Einstein equations. As in the generalized Kasner 
solution, in the general inhomogeneous approach to the singularity the presence of mat- 
ter is exhibited in the relations between the arbitrary spatial functions which appear in 
the solution only. 

6.2 The fragmentation process 

We will now qualitatively discuss a further mechanism that takes place in the inhomoge- 
neous Mixmaster indices in the limit towards the singular point: the so-called fragmenta- 
tion process l30TTl299l 13741 . 

The extension of the BKL mechanism to the general inhomogenous case contains the 
physical restriction of the "local homogeneity": in fact, the general derivation is based on 
the assumption that the spatial variation of all spatial metric components possesses the 
same characteristic length, described by a unique parameter k, which can be regarded as 
an average wavenumber. Nevetheless, such local homogeneity could cease to be valid, as 
a consequence of the asymptotic evolution towards the singularity. 

The conditions (|6.1.3|) do not require that the functions p fl (x 7 ) have the same ordering 
in all points of space. Indeed, they can vary their ordering throughout space an infinite 
number of times without violating conditions (|6.1.3|) , in agreement with the oscillatory- 
like behaviour of their spatial dependence. Furthermore, the most important property 
of the BKL map evolution is the strong dependence on initial values, which produces an 
exponential divergence of the trajectories resulting from its iteration. 

Given a generic initial condition p^(x 7 ), the continuity of the three-manifold requires 
that, at two nearby space points, the Kasner index functions assume correspondingly 
close values. However, for the mentioned property, the trajectories emerging from these 
two values exponentially diverge and, since the p a {x Y ) vary within the interval [—1/3; 1] 
only, indeed the spatial dependence acquires an increasingly oscillatory-like behaviour. 
In the simplest case, let us assume that, at a fixed instant of time to, all the points of the 
manifold are described by a generalized Kasner metric, the Kasner index functions have 
the same ordering point by point, and pi(x 7 ), pi{x y ), P3(x 7 ) are described throughout 
the whole space, by a narrow interval of w-values, i.e. «G [K, K + 1] for a generic integer 
K. We refer to this situation asa manifold composed by one "island". We introduce the 
remainder part of u(x 7 ) as 

X(xT) = u(x^) - \u{x^)] , X G [0, 1) Vx G E, (6.2.1) 
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where the square brackets indicate the integer part. Thus the values of the narrow interval 
can be written as u °(xT) = K° + X°(je?). As the evolution goes by, the BKL mechanism 
induces a transition from an epoch to another; the n-th epoch is characterized by an in- 
terval [K — n, K — n + 1], until K — n = 0; then the era comes to an end and a new one 
begins. The new w^x 7 ) starts from u 1 = 1/X°, i.e., takes value in the interval [1, oo). 
Only very close points can still be in the same "island" of u values; distant ones in space 
will be described by very different integer K 1 and will experience eras of different length. 
As the singular point is reached, more and more eras take place, causing the formation of 
a greater and greater number of smaller and smaller "islands", providing the "fragmenta- 
tion" process. Our interest is focused to the value of the parameter k, which describes the 
characteristic length and increases as the islands get smaller. This implies the progressive 
increase of the spatial gradients and in principle could deform the BKL mechanism. 
We argued that this is not the case with a qualitative explanation: the progressive increase 
of the spatial gradients produces the same qualitative effects on all the terms present in 
the three-dimensional Ricci tensor, including the dominant ones. In other words, for each 
single value of k, in every island, a condition of the form 

inhomogeneous term k 2 t 2Ki f(t) 

~ <C 1 , 

dominant term k 2 ^ 1 

Ki = l- pi > 0, /(f) = 0(lnf, In 2 t), t < 1 , 

is still valid, where we call inhomogeneous the terms containing spatial gradients of the 
scale factors, which are evidently absent from the dynamics of the homogeneous cosmo- 
logical indicess. 

From this point of view, the fragmentation process does not produce any behaviour 
capable of stopping the iterative scheme of the oscillatory regime. 



6.3 Hamiltonian formulation and dry turbulence 



As mentioned above, a generic cosmological solution is represented by a gravitational 
field with all its degrees of freedom and, therefore, allowing to specify a generic Cauchy 
problem. 

In the ADM formalism, the three-metric tensor corresponding to such a generic indices 
reads as 

h Kfi = ei°5 ad O a b O d c d a y%y c , (6.3.1) 

where q a = q a (x,t) and y h = y b (x,t) are six scalar functions and 0\ = 0%(x) is a SO(3) 
matrix. Thus the action for the gravitational field is 



S = [ dtd 3 x (p a d t q a + U d d t y d -NH- N a H a ) , 



H = 



/ZxR 
1 

7^ 



E(p«) 2 4(l>) - h(3)R 



H a = n a d a y a + p a d a q a + IpaiO-^Ol , 
where p a and TI^ are the conjugate momenta of the variables q a and y h , respectively. 



(6.3.2a) 

(6.3.2b) 
(6.3.2c) 
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The ten independent components of a generic metric tensor are represented by the three 
scale factors q", the three degrees of freedom y", the lapse N and the shift-vector N a ; by 
the variation of the action (|6.3.2a|) with respect to p a ,Tl a , the relations 

d t y d = N a 3 a y rf (6.3.3) 

L,aFa \ b b J 

hold. 

We remind that a wide class of cosmological models resembling the behavior of the 
inhomogeneous Mixmaster is the so-called Gowdy cosmology 112281 [279 1. 



6.3.1 Solution of the super-momentum constraint 

By the usual Hamiltonian constraints (|6.3.2b|) and (|6.3.2c[) , the super-momentum can be 
solved by choosing the function y a as space coordinates, taking fj = t, and y a = y a (t,x), 
getting 

l dyf ZMU )a dy* 

and furthermore 



n* = -p.A-2p.(0- 1 )5- s 4. (6.3.5) 



q a (t,x) - 


-> <f(v*y) 




(6.3.6a) 


Va{t,x) - 


-» p' a (y>y) = 


p«(^y)/|/l 


(6.3.6b) 


d 
dt 


dy b d 


3 

dt] 


(6.3.6c) 


d 

dx* 


9y fc 9 
^ 9x<* 3y fc ' 




(6.3.6d) 



where |/| denotes the Jacobian of the transformation, relation (|6.3.6a[) in general holds 
for all the scalar quantities, while (|6.3.6b|) for all the scalar densities and action (|6.3.2a|) 
rewrites as 

S = [ d n d?y (p a d v q a + IpaiO-^Ol - NH) . (6.3.7) 
6.3.2 The Ricci scalar 

In this scheme, the potential term appearing in the super-Hamiltonian reads as 

U = rm (3)R = E A «° 2Q " + E D Q * +a o (d q , (d q ) 2 ,y /V ) (6.3.8) 

\J\ a b^c V 7 

where D = h\}\, h = exp yj H q", the Q fl are the anisotropy parameters (|4.5.4|) and A are the 
functions 

A* = E (o^VO« ( Vy c A Vy") . (6.3.9) 

k,j >■ ' 
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Assuming y a (t,x) smooth enough (which implies smoothness of the coordinates system 
as well), then all the gradients appearing in the potential U are regular, in the sense their 
behavior is not as strongly divergent as to destroy the billiard representation. It was 
shown [299J that the spatial gradients logarithmically increase in the proper time along 
the billiard's geodesies and result of higher order. Thus, as D — > the spatial curvature 
( 3 )R diverges and the cosmological singularity appears; in this limit, the first term of U 
dominates all the remaining ones and can be indicesed by the potential 

li = £®(Q 8 ), (6.3.10) 

a 

resembling the behavior of the Bianchi IX indices. By (|6.3.10|) the Universe dynamics inde- 
pendently evolves in each space point; the point-Universe moves within the dynamically- 
closed domain IIq and near the singularity we have dp a /dt] = 0. Thus, the term2p fl (0 _1 )^9 ; 
in H6.3.7[) behaves as an exact time-derivative and can be ruled out of the variational prin- 
ciple. 

The same analysis developed for the homogeneous Mixmaster model in Section l4~8l can 
be straightforwardly implemented in a covariant way (i.e. without any gauge fixing for 
the lapse function or for the shift vector). 

The super-Hamiltonian constraint is solved in the domain IIq as 



- Pr = -Kadm = W (£ 2 - i)p| + ^rrr < 6 - 3 - n ) 

and the reduced action reads as 

5Sn Q = 5 J dr]d 3 y (pfag + p g d n 6 - Hadm^t) = . (6.3.12) 

By the asymptotic limit (|6.3.10|) and the Hamilton equations associated with (|6.3.12[) we 
get de/dtj = de/drj = and therefore 'KADM^iy") is a constant of motion even in the 
non-homogeneous case. 



6.4 The Iwasawa decomposition 

We will now briefly discuss the generalization of the technique due to Chitre and Misner 
introduced in |164| . This technique was adopted by the authors to study the dynam- 
ics of the multi-dimensional Einstein equations coupled to dilatons and p-forms in the 
neighborhood of a generic singularity (see last section of this chapter). We will limit our 
discussion to the four dimensional case only. 

The Iwasawa decomposition of the spatial metric h a p is based on the choice of a unique 
oriented orthonormal spatial coframe {co a }, co a = e^dx" 1 defined by the relations 

e " a = £D£M* (6.4.1) 

|8 
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resulting in the following expression for 
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m 3 
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1/ 



^ = £exp(-2&")MSMS 



(6.4.2) 



(6.4.3) 



The existence and the uniqueness of this frame is a consequence of the uniqueness of the 
QR decomposition in linear algebra; furthermore the matrix M can be viewed as repre- 
senting the Gram-Schmidt orthogonalization of the spatial coordinate coframe. 

In this variables, the Lagrangian of the gravitational field reads 




^exp[2(b d -¥)]M c K m 



+ \ZhN 3 R, (6.4.4) 



where M denote the inverse matrix of M. 



As soon as the Legendre transformation is taken and the momenta n a and intro- 
duced, the gravitational Hamiltonian in vacuum is obtained 



H = N 



(6.4.5) 



In M164I it is argued that, as soon as the existence of a space-like singularity in the past 
is hypothesized, b a is expected to be timelike in the vicinity of this singularity (in the 
sense of the associated Minkowskian metric that endows the b fl -space). This allows to 
introduce a new set of variables, A and the '"orthogonal angular variables'" 7 H in which 
the Hamiltonian l6.4.5l reads (the lapse function has been rescaled and set equal to exp(2A)) 



-7H 



+ e 



2A 



J^c A exp(- 

A 



(6.4.6) 



here c A denotes some functions of spatial derivatives of the metric, off-diagonal metric 
variables, and momenta; iva(i) denote linear forms of the variables j a , i.e. Wa(j) = 
YLb{ w A)blb- It is expected the singularity to appear in the limit A —> oo; in this limit, 
each term of the summation becomes an infinite wall, and can be described by the usual 
generalized © function. Like in Chite-Misner picture, the are three dominant terms (in 
the vacuum case) that describe the billiard table. The coefficients c A are non negative 
functions of the variables, but "'generically'" they are positive [164]. The three dominant 
terms c^e 2A exp[— 2w^(7)e A ] are assumed to be the only relevant contributions for the 
asymptotic dynamics, while the other are subdominant and dropped away. The limiting 
Hamiltonian then reads 



H, 



asymptotic 



E©[ 

A=\ 



-2w A (y)] 



(6.4.7) 
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As in the previous sections, this Hamiltonian is independent of M, P* and A, suggesting 
the existence of asymptotic constants of the motion. The remaining variables describe free 
motion on the hyperbolic space, constrained by the three walls. 

This Hamiltonian formulation of Einstenian theory was adopted in the cosmological 
framework by Uggla and collaborators in a series of paper [250 . 493. 8J. In particu- 
lar, they split the field equations into decoupled equations for the conformal factor (by 
means of the conformal Hubble-normalization) and a coupled system of dimensionless 
equations for quantities associated with the dimensionless conformal metric. This re- 
duced dimensionless system carries the essential information about the problem, while 
the dimensional one allows to recover the physical metric. One of the main advantages of 
Hubble-normalization lies in the behaviour of the dynamical variables as the cosmologi- 
cal singularity is reached: the Hubble-normalized variables remain bounded despite the 
standard field ones do not (for further details see the recent review [250 J. 
When this formalism is applied to the generic field equations in the vicinity of the initial 
singularity, it can be demonstrated the existence of the so-called billiard attr actor [493 J. 
The main feature of this attractor is to give a precise meaning to the notion of piece-wise 
BKL approximation. This way the space-like cosmological singularity is silent, denoting 
the decoupling of space points on the causal horizon [8J. Furthermore it was established 
the existence of a duality of such framework with the cosmological billiard approach of 
Damour et al. Ifl64| . 

There are some other interesting features, which we mention here as a reference al- 
though demanding to the literature for details. For the Bianchi type IX there is a very 
strong limit on global vorticity, caused by the difficulty of "fitting in" a vortex in the type 
IX compact geometry B246H145ll39ll34| . Moreover, in the inhomogeneous case arises a lin- 
earisation instability [43j, i.e. linearisations around the killing vectors of compact spaces 
are not stable and series expansions around them in general do not form the leading order 
terms of a series converging to a true solution of the Einstein equations, unless staisfying 
further constraints 

6.5 Inhomogeneous quantum Mixmaster 

Since the spatial gradients of the configurational variables don't play a relevant dynami- 
cal role in the asymptotic limit T — > oo (indeed the spatial curvature behaves as a potential 
well), then the quantum evolution takes independently place in each space point and the 
total wave function of the Universe can be represented as 

E(t, u, v) = n. x £ Xl (T, u, v) (6.5.1) 

where the product is (heuristically) taken over all the points of the spatial hypersurface. 
However, in the spirit of the long-wavelength approximation, the physical meaning of 
a space point must be recovered on the notion of a cosmological horizon; in fact, we 
are dealing with regions where the inhomogeneity effects are negligible and this corre- 
sponds to super-horizon sized spatial gradients. Even if this request is statistically well 
confirmed to classical level[299j, in the quantum sector it acquires a precise meaning if 
we refer the dynamics to a lattice indicesing the space-time in which the spatial gradients 
of the configurational variables become real potential terms. In this respect, the geometry 
of the space-time is expected to acquire a discrete structure on the Planck scale and we 
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believe that a regularization of our approach could arrive from a loop quantum gravity 
treatment (HI. 

Despite this local homogeneous framework of investigation, the appearance near the sin- 
gularity of high spatial gradients and of a space-time foam (like as outlined in classical 
dynamics|299 . 374J) can be recognized in the above quantum picture too. In fact, the 
probability that in n different space points (horizons) the variables u and v take values 
within the same narrow interval decreases with n as p n , p being the probability in a sin- 
gle point; all probabilities are identical to each other and no interference phenomenon 
takes place. From a physical point of view, this consideration indicates that a smooth pic- 
ture of the large scale Universe is forbidden on a probabilistic level and different causal 
regions are expected to be completely disconnected from each other during their quan- 
tum evolution. Therefore, if we start with a strongly correlated initial wave function 
So(m, v) = E(to, u, v), its evolution toward the singularity induces increasingly irregular 
distributions, approaching (|6.5.1|) in the asymptotic limit T — > oo. 



6.6 Multidimensional oscillatory regime 

Let us consider a (d + 1) -dimensional space-time (d > 3), whose associated metric tensor 
obeys a dynamics described by a generalized vacuum Einstein equations 

^R lk = 0, (t,* = 0,1,. ..,<*), (6-6.1) 

where the (d + 1) -dimensional Ricci tensor takes its natural form in terms of the metric 
components gik{x l ) . 

In II173| , it is shown that the inhomogeneous Mixmaster behaviour finds a direct general- 
ization in correspondence to any value of d. Moreover, in correspondence to d > 9, the 
generalized Kasner solution acquires a generality character, in the sense of the number of 
arbitrary functions, i.e. without a condition analogous to (\ lambda gen). 
In a synchronous reference (described by the usual coordinates (t, x 7 )), the time-evolution 
of the d-dimensional spatial metric h a p(t,x 7 ) singles out an iterative structure near the 
cosmological singularity (t = 0). Each single stage consists of intervals of time (Kasner 
epochs) during which h a & takes the generalized Kasner form 

M^ 7 ) = X> 2P "'^' ( 6 - 6 - 2 ) 

8=1 

where the Kasner indeces p a (x 7 ) satisfy 

£>(x 7 ) = Zvl{x y ) = 1 ' (6-6-3) 

and Z 1 (x 7 ), ^(x 1 ) denote d linear independent 1-forms l a = m^dx 01 , whose components 
are arbitrary functions of the spatial coordinates. 

In each point of space, the conditions (|6.6.3|) define a set of ordered indexes {p a } (p\ < 
Pi < ■•■ < pi) which, from a geometrical point of view, fixes one point in R d , lying 
on a connected portion of a (d — 2) -dimensional sphere. We note that conditions (|6.6.3|) 
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require pi < and Pd-i > 0/ where the equality takes place for the values p\ = ... = 
Pd-i = and pd = 1 only. The following d-dimensional BKL map, linking the old Kasner 
exponents p a to the new ones q a , holds jl87l ft881 

~Vl - P P2 Pd-2 lc c 

?i=l + 2pi + P' K=i + 2 Vl + P C,d - 2 =l + 2p 1 + P {6 - 6A) 

p A -i +2 Pl + P p d + 2 Pl +P 

qi ~ 1= 1 + 2 P1 +P ' qd= l + 2p 1 + P (6A5) 

where P = Ef=2p«- 

As shown in 1 1741 1261 1 (see also 1299U302| ), each single step of the iterative solution is 
stable, in a given point of the space, if 

limt 2 ^R b a = 0. (6.6.6) 

f->0 

The limit H6.6.6D is a sufficient condition to disregard the dynamical effects of the spatial 
gradients in the Einstein equations. An elementary computation shows how the only 
terms capable to perturb the Kasner behavior in t 2 contain the powers t 2a " bc , where 
ot a }, c are related to the Kasner exponents as 

Kabc = 2p a + Va , (a ^b,a ^ c,b ^ d) , (6.6.7) 

d^a,b,c 

and for generic l a , all possible powers t 2DC,lbc appear in t 2 This leaves two possibilities 

for the vanishing of t 2 ^R a as t — > 0. Either the Kasner exponents can be chosen in an 
open region of the Kasner sphere defined in (|6.6.3[) , so as to make a„f, c positive for all 
triples a, b, c, or the conditions 

* abc (x 7 )>0 V(x\...,x d ) (6.6.8) 

are in contradiction with (|6.6.3[) , and one must impose extra conditions on the functions 
I and their derivatives. The second possibility occurs, for instance, for d = 3, since a abc 
is given by 2p a , and one Kasner exponent is always negative, i.e. Thus (|6.6.3|) is 

a solution of the vacuum Einstein equations to leading order if and only if the vector I 1 
obeys the additional condition 

l x ■ V A h = , (6.6.9) 

and this kills one arbitrary function. 

It can be shown |T74ll26T| that, for 3 < d < 9, at least the smallest of the quantities {6.6.7) . 
i.e. «i,d-i,d results to be always negative (excluding isolated points {p,} in which it van- 
ishes); for d > 10 an open region exists of the (d — 2) -dimensional Kasner sphere where 
a -\,d-i,d takes positive values, the so-called Kasner Stability Region (KSR). 
For 3 < d < 9, the evolution of the system to the singularity consists of an infinite number 
of Kasner epochs, while for d > 10, the existence of the KSR, implies a profound mod- 
ification in the asymptotic dynamics. In fact, the indications presented in [173. 302J in 
favor of the "attractivity" of the KSR, imply that in each space point (excluding sets of 
zero measure) a final stable Kasner-like regime appears. 

In correspondence to any value of d, the considered iterative scheme contains the right 
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number of (d + l)(d — 2) physically arbitrary functions of the spatial coordinates, re- 
quired to specify generic initial conditions (on a non-singular space-like hyper surf ace). 
In fact, we have d 2 functions from the d vectors I and d — 2 Kasner exponents; imparl- 
ance under spatial reparametrizations allows to eliminate d of these functions, and other 
d because of the Oa Einstein equations. This piecewise solution describes the asymptotic 
evolution of a generic inhomogeneous multidimensional cosmological indices. For a cos- 
mological application of the eleven space-times models see [383 J. A wide analysis of the 
dynamics concerning homogeneous and inhomogeneous multidimensional cosmologies 
see 1133611421 13371 l3l2ll478lllMl49Tl l30l 14121 IT23ll3#ll3Ml3^rai3^l3Ml33u1 . 

We summarize some properties about the insertion of p-forms and dilatons in the grav- 
itational dynamics. This is a wide topic, and for review we refer to ||164|| . 

The inclusion of massless p-forms, in a generic multi-dimensional indices |I1601 11611 
I163|) , can restore chaos when it is otherwise suppressed. In particular, even though pure 
gravity is non-chaotic in 11 space-time dimensions, the 3-forms of d + 1 = 11 supergravity lll72| 
make the system chaotic (those p-forms are part of the low-energy bosonic sector of 
super string /M- theory indices ||161H160|| ). 

The billiard description in the four dimensional case is quite general and can be extended 
to higher space-time dimensions, with p-forms and dilatons I28ll ll62ll340| . If there are n 
dilatons, the billiard is a region of the hyperbolic space TI ( j +)J _i, and in the Hamiltonian 
each dilaton is equivalent to the logarithm of a new scale factor. 

The other ingredients that enter the billiard definition are the different types of the walls 
bounding it: symmetry walls related to the off-diagonal components of the spatial metric, 
gravitational walls related to the spatial curvature, and p-form walls (electric and mag- 
netic). All of them are hyper-planar, and the billiard is a convex polyhedron with finitely 
many vertices, some of which are at infinity. 

For the cosmological application of the brane framework[344J, see 12431 14981 1141) , in 
which the homogeneous models and the Mixmaster chaos are discussed. 
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By this work, we resumed many of the key efforts made over the last four decades to 
characterize, as accurate as possible, the nature of the asymptotic behavior of the Uni- 
verse when it approaches the singularity. In particular, the description of the oscillatory 
regime for the Bianchi VIII and IX models outlines that the singular character of the so- 
lution survives also when the isotropy condition is relaxed, but the detailed features of 
the dynamics can be deeply altered, driving up the smooth behaviour of the FLRW scale 
factor to a chaotic evolution of the no longer equivalent directions. 

However, the astonishing result consists of the possibility to extend the oscillatory 
regime to the generic inhomogeneous case, as far as a sub-horizon sized geometry is 
concerned, especially in view of the stability this picture acquires with respect to the 
presence of a matter fluid. We think that this very general dynamics does not provide 
only the important, but somehow academical, proof about the existence of the singularity 
(as a general feature of the Einstein equations under cosmological hypotheses), but it also 
represents the real physical arena to implement any reliable theory of the Universe birth. 
Indeed, both from a classical and from a quantum point of view, the inhomogeneous 
Mixmaster offers a scenario of full generality to investigate the viability of a theoreti- 
cal conjecture, without the serious shortcoming of dealing with specific symmetries (in 
principle even not appropriate to the aim of the implementation, i.e. consider the use of 
the minisuperspace when the quantization is performed, disregarding causality). In this 
sense, the results we presented here are only the basic statements for future developments 
on the nature of the generic singularity. An important issue will be to fix the chaotic fea- 
tures as expectable properties of the Universe origin, when a convincing proposal for the 
quantization of gravity will acquire the proper characteristics of a Theory. Of course, in 
a quantum picture, the chaotic nature of the Universe would be translated into certain 
indicators of the state arrangements. The transition to the classical limit of an expanding 
Universe must be addressed as a crucial stage for the early cosmology investigation. 

On the basis of the material included in this work and in view of the present knowledge 
of fundamental interactions, we can formulate the following proposal for the origin of the 
actual Universe from the generic singularity. Far from being a definitive picture of the 
Universe birth, it better represents the course grain that would be useful, in our opinion, 
as a guideline for future investigations in this field. 

We summarize our point of view into the following individual, but correlated, steps: 

• The Universe was born in a quantum configuration, in which the gravitational field 
and the kinetic term of the scalar inflaton field are the only important degrees of 
freedom. In this scenario, the scalar field would play the role of a time variable (in 
a relational sense) and a Big-Bounce could be inferred. 

• The presence of the scalar field, as well as other pure quantum effects, would imply 
the absence of real quantum chaos features. Thus the Universe could reach configu- 
ration of semiclassical nature, as far as regions of enough large volume are explored. 
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However, we could infer that, without an additional contribution, the classical limit 
could not arise. Two main terms can be suggested as good candidates for this scope, 
the quadratic potential of small oscillations and a cosmological constant term; we 
propose that these two terms could act as synergic ones to each other. 

• After the classical limit is established, the achievement of an isotropic and homoge- 
neous Universe, up to a large scale, is got by the action of the classical slow-rolling 
phase of the scalar field. In this respect it is crucial that the cosmological term can act 
before the three-scalar of curvature becomes important enough, for instance during 
a long era of small oscillations or in the damped oscillation of small anisotropies. 

However, we have to stress that the reliability of this picture is concerned with a strong 
cosmological term, to some extents, significantly greater than the one expected by the 
current theory of inflation. Thus the profile of the inflationary Universe has, in this per- 
spective, to be modified even to include aspects eventually related to a quantum, or semi- 
classical de-Sitter phase. 

Finally, we emphasize that the traced ideas rely on the assumption, made in this review 
article, to deal only with the Einstein-Hilbert action. More general proposals for the quan- 
tum interaction of gravity and matter fields could fix different character for the Universe 
birth, but eventually enforce the presented point of view. 

We would like to thank all the colleagues who gave valuable feedback on our work, 
and in particular, John Barrow, Dieter Lorenz-Petzold and Henk van Elst 
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